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EXECUTIVE  SUMMARY 


Six  DECADE  related  problem  areas  were  the  focus  of  Code  6720  efforts  in  FY  1992  in  the 

Plasma  Radiation  Source  (PRS)  Theory  Program: 

(1)  An  important  need  exists  to  better  understand  and  to  improve  the  performance  of  large  aspect 
ratio  wire  array  implosions,  i.e.,  those  implosions  for  which  the  initial  array  diameter  is  greater 
than  2  cm.  The  ability  to  extend  the  kilovolt  emission  spectrum  in  DECADE  toward  warm 
x-rays  depends  on  the  effectiveness  of  these  implosions.  While  kinetic  energy  arguments 
can  provide  sufficient  guidance  for  scaling  pulsed  power  machine  design,  DNA  is  ultimately 
concerned  with  achieving  the  best  performance  from  the  x-ray  simulators  that  utilize  these 
machines.  In  order  to  begin  an  evaluation  of  the  effectiveness  of  large  aspect  ratio  implosions 
to  convert  kinetic  energy  into  kilovolt  x-rays,  a  series  of  PRS  experiments  was  planned  for 
the  SATURN  machine  at  Sandia  National  Laboratory.  These  experiments  were  carried  out  in 
August,  1992.  They  were  also  carefully  designed  to  increase  the  maximum  yield  of  kilovolt  x- 
rays  from  aluminum  array  implosions  on  SATURN  from  its  previous  high  value,  and  to  allow 
meaningful  comparisons  with  similar  experiments  conducted  on  the  Double  EAGLE  machine 
at  Physics  International,  Inc.  to  be  made.  While  the  experiments  did  achieve  a  significant 
increase  in  kilovolt  yield,  they  did  not  perform  entirely  up  to  expectations,  especially  for  the 
large  aspect  ratio  implosions.  The  design  of  these  and  of  follow-up  experiments  is  discussed 
in  Section  I  of  this  report.  Serious  analysis  of  the  SATURN  experiments  will  take  place  in 
FY93. 

(2)  Because  experimental  implosions  are  much  softer  (i.e.,  much  lower  in  density  and  much 
fatter  in  extent)  than  classically  calculated  1-D  MHD  implosions,  it  is  important  to 
phenomenologically  determine  enhancements  to  classically  calculated  plasma  viscosities, 
heat  conductivities,  and  electrical  resistivities  in  order  to  explain  the  present-day  experiments. 
This  soft  implosion  modeling  is  also  needed  to  provide  a  more  realistic  scaling  of  load 
designs  from  present-day  machines  to  DECADE.  It  was  found  this  year,  for  example,  that  soft 
implosion  modeling  of  argon  implosions  significantly  modifies  the  scaling  of  K-shell  yield 
with  load  mass  and  implosion  velocity.  In  particular,  the  breakpoint  masses  for  achieving 
efficient  K-shell  emission  are  shifted  upwards  from  their  predicted  hard  implosion  values.  It 
was  also  found  that  soft  implosion  modeling  can  achieve  reasonably  good  agreement  with 
experimental  density,  temperature,  and  K-shell  yield  measurements,  and  that  it  can  also 
predict  the  experimental  trends  of  this  data.  This  work  is  described  in  Section  11. 

(3)  As  machine  currents  increase,  load  current  densities  will  increase,  and  the  need  to  determine  the 
current  flow  dynamics  in  the  presence  of  possible  plasma  microturbulence,  runaway  electron, 
and  non-Maxwellian  electron  distribution  effects  could  become  increasingly  important. 
Current  profiles  determine  the  force  dynamics  on  the  PRS  plasma,  the  coupling  of  the  PRS 
load  to  the  generator,  and  the  transfer  of  energy  from  the  electrons  to  the  ions  in  the  load  and 
subsequently  into  x-rays.  Progress  on  the  analysis  of  this  electron  dynamics  in  DECADE-class 
machines  is  discussed  in  Section  133. 

(4)  Because  of  the  finite  number  of  return  current  posts  used  in  PRS  diode?1  and  because  of  the 
three  dimensional  flow  patterns  cf  gas  jets,  exploding  wire  arrays,  ar.d  n  etal  vapor  jets,  2-D 
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MHD  modeling  of  PRS  load  implosions  is  essential  in  order  to  optimize  load  symmetry  and 
kinetic  energy  generation  and  conversion  on  axis.  This  problem  may  be  especially  critical  for 
large  aspect  ratio  implosions.  For  wire  arrays,  r-0  calculations  are  needed  to  determine  the 
optimal  positioning  of  the  wires.  For  gas  puff  implosions,  r-z  calculations  are  needed  to  help 
determine  optimal  nozzle  designs  and  tilt  angles.  This  latter  work  on  Double  EAGLE  and 
SATURN  is  discussed  in  Section  IV.  The  Double  EAGLE  work  has  provided  an  important 
benchmark  for  extending  these  calculations  to  DECADE. 

(5)  Because  DECADE  is  designed  to  be  an  Inductive  Energy  Store  (IES)  generator,  the  effect  of 
Plasma  Opening  Switch  (POS)  operation  on  PRS  performance  must  be  investigated.  There 
are  many  uncertainties  as  to  how  POSs  operate,  and  a  number  of  POS  models  have  been 
developed  to  describe  this  operation.  One  model,  recently  developed  and  suitable  for  use  in 
a  transmission  line  model  description  of  power  flow  in  DECADE,  is  discussed  in  Section 
V.  During  FY92,  work  was  carried  out  to  benchmark  this  model  against  experimental  data 
obtained  from  Code  6770’s  Hawk  POS.  Work  was  also  carried  out  to  demonstrate  a  procedure 
for  optimizing  PRS  performance  within  a  power  flow  calculation  using  any  POS  model.  This 
work  will  be  published  in  the  DoD’s  Journal  of  Radiation  Effects. 

(6)  Finally,  a  small  effort  continued  in  FY92  to  develop  scalable  L-shell  models  for  use  in 
modeling  the  performance  of  PRS  loads  on  DECADE.  It  is  discussed  in  Section  VI.  This 
effort  now  makes  use  of  atomic  code  capabilities  that  were  developed  at  the  Los  Alamos 
National  Laboratory.  There  are  two  problems  of  critical  importance  to  DECADE  that  this 
work  addresses.  One  is  the  problem  of  scaling  L-shell  emission  behavior  with  atomic  number. 
The  other  is  the  problem  of  studying  L-shell  bum-through  in  order  to  ignite  the  K-shell. 


ADVANCED  CONCEPTS  THEORY  ANNUAL  REPORT  1992 

FINAL  REPORT 

I.  DESIGNING  WIRE  ARRAY  PRS  EXPERIMENTS  FOR  PRESENT-DAY  PULSED  POWER 
GENERATORS 

In  March  of  1992,  a  series  of  aluminum  PRS  experiments  were  designed  for  the  SATURN 
machine  at  Sandia  National  Laboratory.  The  purpose  of  these  experiments  was  to  test  the  K-shell 
yield  scaling  that  had  been  predicted  for  aluminum  PRS  loads,1  that  had  been  partially  confirmed  in 
experiments  at  Physics  International  Inc.,2  and  that  were  being  used  to  predict  PRS  performance  on 
DECADE.  The  basic  concept  being  tested  by  the  PI  experiments  was  that,  for  efficient  production 
of  K-shell  emission,  a  z-pinch  array  load  had  to  be  imploded  to  a  final  kinetic  energy  per  ion,  A,, 
in  excess  of  (and  usually  a  small  multiple  of;  a  minimum  energy,  J£m»n,  which  scales  with  the 
atomic  number,  Z,  of  the  load: 


Emin  =  1.012Z*-*"  eV/ian. 

For  aluminum,  this  minimum  energy  is  approximately  12  keV/ion.  The  dimensionless  parameter, 
q,  defined  by  q  =  Ki/Emint  thus  should  satisfy  the  criterion,  q  >  1,  i.e.,  K{  >  12  keV/ion  for 
aluminum. 

In  Ref.  1,  MHD  calculations  had  been  carried  out  that  showed  how  the  K-shell  yield  from 
aluminum  would  scale  with  m,  the  imploded  array  mass  per  unit  length,  when  q  was  held  constant 
or  with  q  when  m  was  held  constant.  These  calculations,  which  were  carried  out  using  classical 
plasma  conductivities,  produced  hard  implosions,  i.e.,  implosions  in  which  the  plasma  assembled 
on  axis  into  a  tight  pinch  of  radius  on  the  order  of  0.1  mm.  The  experiments  that  were  conducted 
at  PI  to  test  the  validity  of  the  calculated  yields  were  not  designed  to  hold  either  m  or  q  constant. 
In  them,  rnrij  was  held  roughly  constant,  where  r0  is  the  initial  radius  of  the  imploded  array. 
This  condition  corresponds  to  holding  either  Kt>  the  total  kinetic  energy  of  the  implosion,  or  the 
product  mq  roughly  constant.  Scaling  relations  had  been  derived  from  the  MHD  calculations  which 
predicted  that  the  yields  in  the  PI  experiments  would  be  roughly  constant  as  well,  corresponding 
to  a  fixed  fractional  (~  30%)  conversion  of  kinetic  energy  into  K-shell  x-rays  for  all  implosions  of 
different  radii,  r0,  for  which  mrjj  =  constant. 

The  PI  experiments  confirmed  that,  when  the  condition,  q  >  1,  was  satisfied,  the  pinch 
behaved  as  an  efficient  bulk  K-shell  emitter  with,  at  minimum,  the  kinetic  energy  conversion 
efficiencies  predicted  by  the  calculations.  However,  two  important  differences  between  theory 
and  experiment  were  observed.  First,  the  K-shell  yields  were  not  flat,  but  had  a  peak  for  small 
q  values  (q  ~  2  or  3),  corresponding  to  a  kinetic  energy  conversion  efficiency  of  >100%,  which 
is  a  theoretical  impossibility.  Second,  the  experimental  implosions  were  much  softer  than  the 
calculated  implosions,  i.e.,  they  assembled  on  axis  with  radii  on  the  order  of  1  mm.  The  first 
difference  suggests  the  conjecture  that  an  anomalously  high  amount  of  Ohmic  heating  occurs  under 
certain  plasma  conditions  on  axis  as  the  current  continues  to  cook  the  pinch.  The  second  points  to 
the  conjecture  that  the  implosions  have  anomalously  high  (MHD  turbulence  induced)  viscosities 
and  heat  conductivities  as  well.  Calculations  utilizing  enhanced  electrical  resistivities,  viscosities, 
and  heat  conductivities  appear  to  confirm  these  conjectures;  however,  more  experimental  guidance 
is  needed  to  confirm  the  use  of  these  effects  in  larger  current  pulsed  power  machines. 

Titus,  the  first  rationale  for  carrying  out  experiments  on  SATURN  was  to  determine  how 
the  experimental  results  of  Double  EAGLE  would  scale  to  higher  currents  and  larger  imploded 
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aluminum  masses.  Because  SATURN  is  a  much  softer  machine  than  Double  EAGLE,  (since 
the  PRS  load  implosion  dynamics  reacts  back  on  the  current  flowing  through  the  load  much 
more  strongly  in  SATURN  than  in  Double  EAGLE),  it  was  important  to  design  the  SATURN 
experiments  using  the  lumped  circuit  model  for  SATURN.  The  same  procedure  would  be  needed 
for  Double  EAGLE  if  greater  control  of  the  experimental  parameters  is  required  than  was  achieved 
in  the  Ref.  2  experiments. 

The  circuit  model  that  was  used  to  design  the  SATURN  experiments  is  shown  in  Figure  1 
of  Ref.  3.  The  pulsed  power  generator  is  described  as  a  time  dependent  voltage  source,  V'(t), 
driving  a  line  resistance,  Zq,  a  line  inductance,  Lpt  and  a  dynamic  imploding  z-pinch  load  with 
a  time  dependent  inductance,  L(r(t)).  By  using  a  slug  model  description  of  the  dynamic  load 
to  determine  r(f),  the  location  of  die  shell  of  imploding  plasma,  the  PRS  load  dynamics  can  be 
initially  simplified  to  the  problem  of  solving  the  following  two  equations  of  motion:3 


dl  (  l  dr  \ 


(1) 


d?r 

mdt*  = 


2£  r’ 


(2) 


where  I(t)  is  the  current  flowing  through  the  load,  Lt  =  Lp  —  L0in(r/r0),  Lo  = 
(1  —  l/JV)(/i0f)/(2ir),  N  is  the  number  of  wires  in  the  array,  t  is  the  length  of  the  array, 
and  no  —  47rl0-7  henry/m.  The  inductance  of  the  load  is  given  by  L  =  -L0ln(T/r0).  A  set 
of  four  benchmark  calculations  were  carried  out  to  test  the  use  of  this  circuit  model  to  explain 
previous  PRS  implosion  data  that  had  been  taken  at  the  SATURN  facility.4  The  benchmark  worked 
well.  The  predicted  implosion  times,  total  implosion  kinetic  energies,  and  kinetic-energies-per-ion 
were  in  accord  with  experimental  observations. 

Three  sets  of  experiments  were  planned.  They  were  designed  using  commercially  available 
aluminum  wire  sizes,  and,  in  an  attempt  to  achieve  maximum  implosion  symmetry,  the  arrays 
were  taken  to  contain  24  wires.  The  first  set  of  experiments  were  designed  to  closely  parallel  the 
Double  EAGLE  experiments  so  that  direct  inferences  as  to  how  the  K-shell  yields  were  scaling  to 
higher  mass  loads  could  be  made.  Therefore,  as  in  the  Double  EAGLE  experiments,  the  initial 
array  diameters  were  taken  to  be  <  2  cm.  In  the  first  set  of  experiments,  the  total  kinetic  energy 
generated  in  each  implosion  was  170  kJ  and  tj’s  ranged  in  value  from  less  than  1  to  6.  If  30% 
of  the  kinetic  energy  could  be  converted  to  K-shell  emission,  the  minimum  expected  output  of 
K-shell  x-rays  in  these  experiments  would  be  51  kJ.  Eqs.  (1)  and  (2)  describe  the  dynamics  of 
hard  MHD  implosions  only  to  the  point  where  the  plasma  is  close  enough  to  the  axis  to  generate  a 
back  pressure  comparable  to  the  j  x  B  forces,  which  are  described  by  the  I2  term  on  the  right  side 
of  Eq.  (2).  Thus,  the  solutions  to  Eqs.  (1)  and  (2)  must  be  terminated  before  this  condition  ensues. 
The  parameters  of  the  first  set  of  SATURN  experiments  that  are  given  in  the  following  table  were 
calculated  from  Eqs.  (1)  and  (2)  by  terminating  the  solution  in  each  case  at  r(Lmp)  =  1.5  mm, 
where  Ump  is  by  definition  the  implosion  time. 


The  calculations  in  Ref.  1  also  suggested  that  K-shell  yields  would  continue  to  rise  as  a 
function  of  ij  for  a  fixed  m  until  fairly  large  t)  values  were  reached.  For  example,  peaks  in  the 
calculated  yields  ocurred  at  rj  =  7.5  and  rj  =  27.5  for  m  —  60  and  200  /ig/cm  respectively. 
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SATURN  PRS  Experiments 

K.E.=  170  kJ.  Return  current  radius  =1.4  cm,  array  length  =  2.0  cm 


Wire  Diam 

(mil) 

Wire  Mass 
(pg/cm) 

Array  Rad 

(cm) 

H 

#  Wires 

2.5 

2050 

0.638 

0.95 

24 

2.0 

1310 

0.67 

1.6 

24 

1.7 

950 

0.7 

2.2 

24 

1.5 

740 

0.73 

wm 

24 

1.2 

470 

0.79 

4.2 

24 

1.0 

330 

0.86 

6.0 

24 

Such  large  t/  values  could  be  reached  on  the  SATURN  generator  by  going  to  large  aspect  ratio 
implosions,  i.e.,  those  implosions  for  which  the  initial  array  diameter  was  larger  than  2  cm.  In 
argon  and  krypton  gas  puff  experiments,  large  aspect  implosions  had  been  observed  to  produce  the 
highest  yields.  This  observation  provided  a  second  rationale  for  proposing  a  second  set  of  large 
aspect  ratio  SATURN  experiments  in  which  the  total  kinetic  energy  generated  in  the  implosions 
was  held  fixed  at  230  kJ.  A  30%  conversion  of  this  energy  to  kilovolt  x-rays  would  produce  a 
K-shell  yield  of  69  kJ.  The  following  table  lists  the  calculated  parameters  for  these  large  aspect 
ratio  aluminum  array  implosions: 

SATURN  PRS  Experiments 

K.E.=  230  kJ.  Return  current  radius  =  2.4  cm,  array  length  =  2.0  cm 


Wire  Diam 
(mil) 

Wire  Mass 
(pg/cm) 

Array  Rad 
(cm) 

H 

#  Wires 

1.5 

740 

1.1 

3.5 

24 

1.2 

470 

1.145 

5.6 

24 

1.0 

330 

1.21 

8.3 

24 

0.75 

185 

1.365 

15 

24 

0.7 

160 

1.405 

17 

24 

0.6 

120 

1.52 

22 

24 

A  third  set  of  special  case  experiments  was  also  proposed.  For  example,  by  going  out  to 
a  radius  of  ro  =1.7  cm,  a  large  implosion  velocity  and  a  total  kinetic  energy  of  280  kJ  could 
be  generated  using  either  aluminum  or  titanium  wires.  The  final  velocity  achievable  in  such  an 
implosion  would  be  sufficient  to  generate  an  t)  larger  than  1  for  titanium  (for  which  Emin  =  83.4 
keV/ion)  provided  the  finest  0.8  mil  titanium  wire  is  used.  The  masses  of  these  two  implosions, 
listed  in  the  following  table,  were  large  enough  to  anticipate  that  good  conversion  of  kinetic  energy 
into  kilovolt  x-rays  was  possible  following  assembly  on  axis  of  the  pinch: 
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K.E.s=  280  kJ. 


SATURN  PRS  Experiments 
Return  current  radius  =  3.0  cm,  array  length  =  2.0  cm 


Element 

Wire  Diam 
(mil) 

Wire  Mass 
(fig/cm) 

Array  Rad 

(cm) 

B 

#  Wires 

A1 

1.2 

473 

1.7 

n 

24 

Ti 

0.8 

350 

1.7 

ESI 

24 

Two  other  special  array  experiments  were  designed  to  simulate  an  exploding  wire  by 
positioning  the  wires  so  close  to  the  axis  that  they  would,  immediately  upon  exploding,  form  a 
plasma  on  axis  with  sufficient  back  pressure  to  negate  the  generation  of  kinetic  energy.  The  wire’s 
mass  loading  was  chosen  so  that  the  wires  would  quickly  assemble  on  axis  into  an  anticipated 
Bennett  equilibrium.  Under  these  conditions,  the  ability  of  Ohmic  heating  to  ionize  the  bulk 
plasma  into  the  K-shell  could  be  observed  as  a  test  of  the  conjecture  that  a  strong  anomalous 
heating  could  be  generated  in  a  z-pinch  on  axis.  The  parameters  that  were  chosen  for  these  two 
experiments  are : 


SATURN  PRS  Experiments 

K.E.=  0  kJ.  Return  current  radius  =1.4  cm,  array  length  =  2.0  cm 


Element 

Wire  Diam 

(mil) 

Wire  Mass 
(fig/ cm) 

Array  Rad 
(cm) 

B 

A1 

1.7 

950 

0.2 

0 

24 

K.E.=  0  kJ.  Return  current  radius  =  0.5  cm,  array  length  =  2.0  cm 


Element 

Wire  Diam 
(mil) 

Wire  Mass 

(fig! cm) 

Array  Rad 

(cm) 

V 

#  Wires 

A1 

2.0 

1310 

0.2 

0 

24 

Finally,  the  array  parameters  for  one  last  (large  aspect  ratio)  special  case  were  calculated.  The 
goal  in  this  experiment  was  to  implode  a  3  cm  long  aluminum  array  of  sufficient  mass  and  77  in 
order  to  achieve  the  largest  possible  yield  in  these  experiments.  The  total  kinetic  energy  generated 
in  this  case  is  330  kJ;  a  30%  conversion  efficiency  would  yield  99  kJ  of  kilovolt  x-rays.  The 
parameters  for  this  case  are: 
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SATURN  PRS  Experiments 
Return  current  radius  =  2.4  cm,  array  length  =  3.0  cm 


Element 

Wire  Diam 
(mil) 

Wire  Mt,ss 
(/ig/cm) 

Array  Rad 

(cm) 

n 

#  Wires 

A1 

1.2 

473 

1.5 

5.8 

24 

Shortly  after  the  above  experiments  had  been  completed,  R.  Spielman  at  Sandia  indicated 
that  a  significant  improvement  in  aluminum  K-sheli  yield  had  been  attained  on  Saturn  in  the  0.86 
cm  radius,  330  /tg/cm  shot  (an  increase  of  30  to  40%  over  previous  yields).  An  improvement  of 
this  order  was  indicated  by  1-D  MHD  calculations  that  focused  on  the  thermalization  of  kinetic 
energy  as  the  primary  mechanism  for  producing  K-shell  emission.  These  same  calculations  predict 
that  the  yield  should  continue  to  increase  for  the  larger  aspect  ratio  implosions  (array  diameters 
>  2  cm);  however,  R.  Spielman  also  indicated  that  these  large  aspect  ratio  A1  experiments,  just 
completed  in  August,  1992,  did  not  work  as  well  as  hoped.  When  these  experiments  were  planned 
in  March  the  goal  was  to  tradeoff  load  mass  for  implosion  velocity,  using  the  early  hard  implosion 
calculations  as  a  guide.  They  had  showed  K-shell  yield  increasing  as  a  function  of  77  up  to  fairly 
high  values.  Implosion  symmetry  is  very  important  in  implosions  that  start  from  large  radii, 
however,  and  it  was  hoped  that,  by  using  24  wires,  the  SATURN  experiments  had  a  better  chance 
of  success  than  earlier  large  radius  experiments  on  Double  EAGLE  or  BJS.  However,  the  pinhole 
pictures  of  these  high  aspect  ratio  implosions  apparently  showed  poor  implosion  quality. 

Based  on  these  very  preliminary  and  undocumented  findings,  a  new  set  of  aluminum  PRS 
experiments  was  planned  for  the  Phoenix  facility  at  the  Naval  Surface  Weapons  Center.  A  circuit 
model  for  Phoenix  was  obtained  from  P.  Corcoran  of  PSI.  When  run  in  short  circuit  mode,  it 
produced  about  5.3  MA  peak  current  with  a  risetime  of  approximately  100  ns.  Efforts  to  design 
some  experiments  for  Phoenix  were  then  begun.  They  were  influenced  by  two  other  important 
recent  findings: 

(1 )  Recent  large  aspect  ratio  experiments  with  krypton  at  PI  were  very  successful  (in  contrast 
to  A1  on  SATURN),  i.e.,  the  peak  in  the  kilovolt  yield  curve  (L-shell  radiation)  occurred  at  a  4  cm 
nozzle  diameter.  This  is  a  larger  diameter  than  was  used  on  any  of  the  Sandia  wire  experiments. 
Therefore,  all  loads  being  equal,  the  Sandia  experiments  should  have  mimicked  the  success  of  the 
PI  krypton  experiments.  Apparently,  even  24  wires  with  6  or  8  return  current  posts  do  not  perform 
as  well  as  one  cylindrical  flow  of  gas  with  the  same  return  current  posts.  The  important  question 
is  why  not? 

(2)  Calculations  on  argon  that  were  carried  out  for  the  BEAMS  92  conference  in  May  showed 
that,  for  soft  implosions,  load  mass  was  more  important  than  77,  i.e.,  the  K-shell  yield  curves 
saturated  rapidly  with  t?  (77  <  6)  for  the  kinds  of  masses  used  in  the  SATURN  and  PI  experiments. 
(For  the  PI  krypton  experiments,  we  found  that  the  L-shell  yield  peaked  for  77  —  1.)  These  soft 
implosion  (also  called  turbulent)  results  are  shown  in  Figure  1  in  comparison  to  the  200  /ig/cm  hard 
implosion  (non-turbulent)  result.  In  a  soft  implosion  calculation  77*  represents  the  (normalized) 
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total  energy  coupled  to  the  load:  kinetic  +  pdV  work  +  Ohmic  heating.  As  Figure  2  shows,  in  a 
soft  implosion,  only  55  to  60%  of  the  coupled  energy  is  kinetic,  which  accounts  for  the  softness  of 
the  implosion. 

We  conclude  from  these  findings  that  large  aspect  ratio  wire  experiments  should  be  carried  out 
in  the  future  at  SATURN,  Phoenix,  and  Double  EAGLE  to  test  a  different  machine/load  tradeoff: 
that  between  mass  and  implosion  radius.  Experiments  in  which  tj  is  fixed  at  small  values  from 
2  to  6  should  be  done  with  an  increasing  array  mass  (or  increasing  number  of  wires)  in  order 
to  generate  increasing  kinetic  energies.  A  series  of  such  Phoenix  experiments  were  designed  in 
September.  Like  the  SATURN  experiments,  they  utilize  standard,  commercially  available,  wire 
sizes,  and  they  are  tailored  to  the  fact  that  Phoenix  has  eight  return  current  posts.  The  proposed 
aluminum  wire  array  parameters  are  listed  in  the  following  three  tables: 


n  =  2, 


Phoenix  Aluminum  Experiments 
Return  current  radius  =  3.0  cm,  array  length  =  2.0  cm 


Wire  Diam 

Wire  Mass 

Array  Rad 

K.E. 

#  Wires 

(mil) 

(Mg/cm) 

(cm) 

(kJ) 

1.0 

219 

0.75 

38 

16 

1.5 

246 

0.79 

42 

8 

1.2 

315 

0.29 

54 

16 

1.7 

316 

0.9 

54 

8 

2.0 

438 

1.15 

74 

8 

1.2 

473 

1.325 

24 

Return  current  radius  =  3.0  cm,  array  length  =  2.0  cm 

Wire  Diam 

Wire  Mass 

Array  Rad 

K.E. 

#  Wires 

(mil) 

(/xg/cm) 

(cm) 

(kJ) 

1.2 

158 

0.87 

40 

8 

0.75 

185 

0.93 

47 

24 

1.0 

219 

1.02 

57 

16 

1.5 

246 

1.1 

63 

8 

1.2 

315 

1.32 

81 

16 

These  experiments  differ  from  the  SATURN  experiments  because,  as  the  array  diameter  is 
increased  to  generate  more  total  kinetic  energy,  the  r\  values  (i.e.  the  final  implosion  velocities) 
are  kept  moderately  low  while  the  array  mass  is  increased.  The  large  array  diameter  SATURN 
experiments  were  designed  for  roughly  the  same  mass  values  that  are  shown  in  the  above  tables, 
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Phoenix  Aluminum  Experiments 
r)  —  A,  Return  current  radius  =  3.0  cm,  array  length  =  2.0  cm 


Wire  Diam 

(mil) 

Wire  Mass 

0ig/cm) 

Array  Rad 

(cm) 

K.E. 

(kJ> 

#  Wires 

1.2 

158 

1.08 

54 

8 

0.7 

161 

1.08 

55 

24 

0.75 

185 

1.17 

63 

24 

1.0 

219 

1.32 

75 

16 

1.5 

246 

1.55 

85 
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but  for  much  larger  rjs.  Again,  earlier  aluminum  experiments  at  Double  EAGLE  had  shown  peak 
K-shell  yields  for  17  values  of  2  to  3,  and  they  converted  close  to  100%  of  the  implosion  kinetic 
energy  to  K-shell  x-rays.  If  similar  results  could  be  obtained  in  these  Phoenix  experiments,  they 
might  produce  maximum  K-shell  yields  of  between  70  to  80  kilojoules. 

Because  it  is  important  to  discover  why  K-shell  yields  tend  to  fall  off  as  the  array  diameter  is 
increased,  it  was  also  proposed  that  one  or  more  of  the  large  aspect  ratio  implosion  cases  in  the 
above  list  of  experiments  be  carried  out  in  two  geometries.  These  geometries  are  displayed  in  the 
following  three  pictures. 


16  CYLINDRICAL  WIRES 

.  • 


It  is  standard  practice  to  mount  the  array  wires  on  the  same  circle  as  shown  in  the  first  picture. 
However,  in  this  case,  the  wires  nearest  the  return  current  posts  will  experience  the  strongest 
j  x  B  forces  and  will  be  accelerated  to  the  axis  sooner  than  the  other  wires.  Thus,  a  symmetric 
beginning  for  the  wires  will  not  lead  to  a  symmetric  ending  on  axis,  a  condition  that  is  essentia1  for 


18  CONCENTRIC  WIRES,  8  WIRE  OFFSET 


24  CONCENTRIC  WIRES,  16  WIRE  OFFSET 


the  efficient  production  of  K-shell  radiation.  For  this  reason,  we  propose  that  a  second  geometry 
be  employed,  illustrated  above  in  the  lower  two  pictures  depending  on  the  number  of  wires.  In 
this  case,  the  wires  are  mounted  on  two  concentric  circles.  The  wires  nearest  the  return  current 
posts  are  mounted  on  the  outer  circle,  and  the  other  wires  are  mounted  nearer  to  the  axis  on  the 
inner  circle.  The  goal  of  these  experiments,  like  the  goal  of  the  tilted  nozzle  puff  gas  experiments 
recently  conducted  at  PI,  is  to  determine  hether  the  K-shell  yield  will  increase  as  the  inner  circle 
radius  is  reduced.  One  might  begin  with  a  1/2  mm  radial  separation  between  the  inner  and  outer 
circles,  and  then  double  it  to  1  mm,  and  then  double  it  again  to  2  mm.  for  example. 

A  circuit  model  for  Double  EAGLE  was  also  used  to  design  a  set  of  experiments  for  Double 
EAGLE  that  are  the  equivalent  of  the  above  Phoenix  experiments.  They  take  into  account  the  fact 
that  Double  EAGLE  has  only  6  return  current  posts,  while  Phoenix  has  8,  and  they  again  make 
use  of  commercially  available  wire  sizes.  Since  the  idea  for  these  experiments  is  to  get  large 
aspect  ratio  aluminum  experiments  to  work  more  efficiently,  the  comparison  of  well  diagnosed 
Double  EAGLE  with  comparable  Phoenix  experiments  should  help  us  to  understand  what  the 
problems  are  in  achieving  this  performance.  As  in  the  Phoenix  experiments,  it  will  be  important  to 
understand  the  role  of  mass  /  kinetic  energy  tradeoffs  in  achieving  these  efficiency  improvements. 
If  these  experiments  can  be  made  to  work  as  well  as  earlier  PI  experiments,  then  100%  of  the 
kinetic  energies  listed  below  would  be  converted  to  K-shell  x-rays,  and  Double  EAGLE  would 
nearly  double  its  previous  maximum  A1  K-shell  yield.  The  proposed  experiments  are  listed  in  the 
following  three  tables: 

Because  Phoenix  is  a  newly  designed  and  operated  machine,  its  lumped  circuit  model  has  yet 
to  be  tested.  Another  reason,  therefore,  for  comparing  Double  EAGLE  and  Phoenix  experiments 
is  to  obtain  an  important  benchmark  for  the  Phoenix  circuit  by  way  of  Double  EAGLE.  The  earlier 
experiments  that  were  carried  out  at  Double  EAGLE  were  both  well  diagnosed  and,  consequently, 
well  benchmarked.  Hence,  they  have  provided  a  valuable  validation  of  the  Double  EAGLE  lumped 
circuit  model.  We  demonstrate  this  fact  with  the  following  example.  Additionally,  it  illustrates 
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Double  EAGLE  Aluminum  Experiments 


*7  =  2,  Return  current  radius  =  2.0  cm,  array  length  =  2.0  cm 


Wire  Diam 
(mil) 

Wire  Mass 

Ozg/cm) 

Array  Rad 

(cm) 

K.E. 

(kJ) 

#  Wires 

1.0 

164 

0.8 

28 

12 

1.2 

236 

0.91 

41 

12 

2.0 

328 

1.07 

56 

6 

1.2 

355 

1.14 

61 

18 

*7  =  3,  Return  current  radius  =  2.0  cm,  array  length  =  2.0  cm 


Wire  Diam 
(mil) 

Wire  Mass 
(pg/cm) 

Array  Rad 

(cm) 

K.E. 

(kJ) 

#  Wires 

1.2 

118 

0.92 

30 

6 

1.0 

164 

1.04 

43 

12 

1.5 

185 

1.09 

47 

6 

1.2 

236 

1.26 

61 

12 

1.0 

246 

1.3 

63 

18 

*7  =  4,  Return  current  radius  =  2.0  cm,  array  length  =  2.0  cm 


Wire  Diam 
(mil) 

Wire  Mass 
Oxg/crn) 

Array  Rad 

(cm) 

K.E. 

(kJ) 

#  Wires 

0.7 

121 

1.115 

41 

18 

0,75 

139 

1.18 

47 

18 

1.0 

164 

1.3 

52 

12 

1.5 

185 

1.4 

64 

6 

that  die  full  capacity  of  Double  EAGLE  to  produce  K-shell  emission  from  aluminum  loads  has 
yet  to  be  realized.  The  aluminum  experiments  proposed  above  should,  hopefully,  demonstrate  this 
fact.  If  not,  then,  they  ought  to  provide  us  with  some  reasons  why  not. 

In  the  following  example,  the  circuit  equations,  (1)  and  (2)  above,  are  used  with  the  circuit 
elements  appropriate  to  Double  EAGLE.  A  systematic  examination  of  the  implosion  dynamics  up 
to,  but  just  prior  to,  assembly  on  axis,  is  then  carried  out  by  varying  the  mass  of  the  aluminum  array. 
The  experimental  case  in  Ref.  (2)  being  modeled  is  the  one  where  164  pg/cm  was  imploded  from 
a  radius  of  0.75  cm.  From  Ref.  (2),  the  implosion  time  was  observed  to  be  95  ns  and  the  implosion 
produced  28  kJ  of  K-shell  yield  with  an  *?  value  of  1.7.  The  results  of  the  calculations  that  were 
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performed  with  the  Double  EAGLE  slug  model  circuit  equations  are  shown  in  Figs.  (3)-(6). 

As  more  mass  is  imploded  from  the  same  initial  array  radius  (of  0.75  cm),  one  sees  from 
Fig.  (3)  that  the  implosion  time  is  stretched  from  1 10  to  260  ns.  Small  mass  implosions  reach  the 
axis  early  during  the  rise  of  the  current  pulse;  whereas,  the  large  mass  implosions  reach  the  axis 
late  in  the  current  pulse  after  it  has  peaked  (see  Fig.  (4)).  The  labels,  R,  L,  and  G,  in  Figs.  (3)-(6) 
place  into  correspondence  the  mass  loadings  with  the  three  current  traces  and  kinetic  energies  that 
are  generated  when  the  masses  implode.  Fig.  (5)  shows  that  the  implosion,  L,  that  produces  the 
most  kinetic  energy  (~  50  kJ)  corresponds  to  an  implosion  that  occurs  just  after  peak  current  at 
roughly  167  ns.  However,  in  Fig.  (6),  one  sees  that  this  implosion  has  an  t?  value  of  only  0.4 
and,  therefore,  it  should  be  a  poor  radiator  of  kilovolt  x-rays.  Thus,  even  though  this  particular 
choice  of  array  parameters  produces  good  coupling  to  the  Double  EAGLE  generator  and  a  near 
optimal  amount  of  kinetic  energy  generation,  it  also  produces  a  poor  kilovolt  x-ray  radiator.  The 
Double  EAGLE  load  designs  presented  in  the  above  tables  were  clearly  designed  to  overcome  this 
limitation.  They  were  determined  to  optimize  kinetic  energy  as  well  as  kilovolt  x-ray  production. 
Unfortunately,  the  present  impediment  to  their  successful  realization  appears  to  be  the  requirement 
that  they  are  initiated  from  array  diameters  larger  than  2  cm.  An  important  near  term  objective  of 
the  DECADE  PRS  program  should  be  to  remove  this  limitation. 

Fig.  (3)  also  shows  an  extrapolation  of  the  mass  versus  implosion  time  curve  to  approximately 
the  mass  value  of  164  /ig/cm.  These  slug  model  calculations  then  predict  that  this  mass  should 
implode  to  the  axis  in  roughly  105  ns.  As  seen  in  Fig.  (4),  this  time  can  easily  be  adjusted  down  to 
95  ns.  given  the  experimental  uncertainty  of  identifying  the  beginning  of  the  current  pulse.  When 
the  kinetic  energy  curve  is  extrapolated  to  105  ns  in  Fig.  (5),  it  shows  that  25  kJ  of  kinetic  energy  is 
generated  in  the  experiment  corresponding  to  a  greater  than  100%  efficiency  of  converting  kinetic 
energy  to  kilovolt  x-rays.  Fig.  (6)  confirms  that  a  164  /ig/cm  aluminum  implosion  from  a  radius 
of  0.75  cm  has  enough  implosion  velocity  to  be  a  good  K-shell  radiator,  i.e.,  in  this  case,  r)  >  1 .7. 

Enhanced  transport  MHD  calculations  were  performed  to  test  the  quality  of  the  energy 
conversion  to  be  expected  in  both  the  proposed  Phoenix  and  Double  Eagle  aluminum  experiments. 
These  calculations  were  done  for  soft  implosions  in  accord  with  the  HEART  Conference  work  to 
be  published  in  the  DoD  Journal  of  Radiation  Effects,  i.e.,  the  thermal  conductivity  was  enhanced 
by  a  factor  of  30,  electrical  resistivity  was  enhanced  by  20,  and  the  viscosity  was  enhanced  by 
25.  The  K-shell  emission  comes  predominantly  from  thermalization  of  kinetic  energy  because  the 
current  is  terminated  when  the  voltage  across  the  plasma  diode  reaches  3  MV.  The  Double  Eagle 
calculations  show  that  K-shell  yields  of  20  kJ/cm  can  be  achieved.  This  is  specifically  true  for 
the  proposed  experiment  in  which  the  mass  loading  is  246  fig/ cm  and  the  array  radius  is  1.3  cm. 
Similar  calculations  made  for  previous  Double  Eagle  experiments  underpredicted  the  measured 
K-shell  yield  (by  as  much  as  7  kJ/cm),  especially  for  loads  in  excess  of  200  fig/cm.  It  is  thought 
that  current-on  effects  such  as  Ohmic  heating  are  responsible  for  this  under-prediction.  If  these 
same  effects  are  present  in  the  above  proposed  experiments,  then  one  would  predict  that  K-shell 
yields  of  25  to  30  kJ/cm  ought  to  be  achieved.  Similar  calculations  for  Phoenix  show  that,  without 
taking  into  consideration  any  current-on  effects,  K-shell  yields  of  30  kJ/cm  are  possible.  Hus 
predicted  yield  is  larger  than  the  Double  Eagle  yield  because  of  the  larger  current  that  Phoenix 
produces  at  current  tum-off  time.  If  these  predictions  are  not  observed,  then  they  may  provide  an 
indicator  for  determining  the  propensities  for  the  two  machines  to  short-circuit  in  the  diode  region. 
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The  above  yield  results  represent  what  is  believed  to  be  achievable  on  Double  EAGLE  and 
Phoenix  given  that  the  load  design  promotes  a  symmetric  implosion.  Although  detailed  data 
analysis  still  needs  to  be  performed  on  the  SATURN  aluminum  experiments,  it  is  believed  that 
asymmetric  behavior  may  be  responsible  for  the  fact  that  the  larger  aspect  ratio  wire  loads  did  not 
perform  as  well  as  expected.  Although  large  aspect  ratio  argon  and  krypton  gas  puff  experiments 
at  PI  did  exhibit  symmetric  implosion  behavior  and  produce  K-shell  and  L-shell  yields  as  large  or 
larger  than  predicted,  these  encouraging  results  are  not  yet  reproducible  with  wire  arrays.  Thus, 
much  experimental  and  theoretical  load  design  work  needs  to  be  done  if  DECADE  is  to  reach  its 
full  potential  as  an  x-ray  simulator,  especially  using  wire  arrays. 
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Figure  1. 


•  Loads  are  driven  by  linearly  ramped  current  profiles 

from  an  initial  radius  of  1.0  cm.  The  current  is  terminated 
at  the  time  a  0-D  slug  calculation  predicts  the  outer  radius 
of  the  load  is  within  0.14  cm  of  the  axis. 

•  r\*  is  the  coup  led  energy-per-ion  normalized  by  the  minimum 
energy  needed  to  ionize  to  the  K-shell  (38  keV  for  argon). 

rjm  =  (JxB  work  +  ohmic  heat  i ng)/(number  of  ions  x  38  keV) 

•  These  K-shell  yields  also  include  >  IkeV  radiation,  however, 
the  >  3  keV  radiation  represents  at  least  85%  of  the  K-shell 
yield. 
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Ratio  of  Peak  Kinetic  Energy  to 
Coupled  Energy  Versus  Mass  and 
Coupled  Energy  (77*) 


Non-Turbulent  200  jig/cm 


Turbulent  600  #ig/cm 


Turbulent  200  fig/cm 
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Figure  2. 

rj *  is  the  coupled  energy-per-ion  normalized  by  the  minimum 
energy  needed  to  ionize  to  the  K-shell  (38  keV  for  argon). 

rj*  =  (JxB  work  +  ohmic  heat ing)/(number  of  ions  x  38  keV) 


is  the  peak  kinet  ic-energy-per-ion  normalized 
by  the  minimum  energy  needed  to  ionize  to  the 
K-shell  (38  keV  for  argon). 


current  traces  of  the  different  mass  implosions  of  Fig.  (3)  are  shown  with  three 
G,  enlarged  to  allow  easy  correspondence  to  be  made  with  the  masses  in  Fig.  (3). 
radius  is  0.7S  cm,  the  array  length  is  2  cm,  and  the  return  current  radius  is  2  cm. 
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Figure  5.  The  total  kinetic  energies  generated  in  eac 
shown  as  a  function  of  implosion  time.  The  initial  am 
cm,  and  the  return  current  radius  is  2  cm. 


II.  COMPARISON  of  PI  ARGON  EXPERIMENTS  WITH  1-D 
ENHANCED  TRANSPORT  CALCULATIONS 

INTRODUCTION 


1  2 

In  work  that  was  presented  at  the  1992  HEART  and  BEAM  conferences,  *  a 
phenomenological  treatment  of  hydromagnetic  turbulence  and  microturbulence 
was  included  in  the  standard  non- turbulent  1-D  MHD  model  in  order  to  obtain 
better  agreement  with  experimental  measurements  than  had  been  provided  by 
the  standard  1-D  model.  Historically  the  non-turbulent  model  over-predicts 
the  experimental  stagnation  densities  and  radiative  powers  by  orders  of 
magnitude.  The  phenomenological  treatment  consisted  of  enhancing  transport 
by  using  multipliers  for  the  classical  transport  coefficients  of 
resistivity,  viscosity,  and  thermal  conduction.  A  choice  of  coefficients  was 
made  so  that  the  calculated  stagnation  ion  densities  and  electron 
temperatures  approximated  those  measured  in  a  series  of  Double  Eagle  argon 
experiments . 

The  results  of  the  HEART  and  BEAM’S  work  showed:  (1)  Because  of  the 
softer  implosions  (lower  stagnation  ion  density)  generated  by  enhancing  the 
transport  coefficients,  the  turbulent  calculations  require  more  mass  than 
non-turbulent  ones  to  reach  an  efficient  K-shell  emission  scaling  regime  (> 
202  conversion  of  kinetic  energy  into  K-shell  photons).  (2)  The  range  of 
values  for  If  (kinetic-energy-per-ion  normalized  to  the  amount  of  energy 
required  to  ionize  to  the  K-shell)  over  which  efficient  scaling  can  be 
obtained  for  a  given  mass  is  limited  to  low  values  in  the  enhanced 
transport  scaling  calculations.  (3)  When  turbulence  is  included,  the  onset 
of  K-shell  burn-through  and  lower  plasma  density  severely  limit  the  maximum 
K-shell  yield  that  is  attainable  from  a  given  mass  (as  opposed  to  the  non- 
turbulent  case).  However,  other  than  these  quantitative  differences  between 
turbulent  and  non-turbulent  transport  models,  the  results  of  earlier  non- 
turbulent  transport  scaling  work  are  still  largely  valid.  In  particular,  the 
existence  of  an  efficient  scaling  regime  for  producing  K-shell  radiation, 
which  is  a  function  of  mass  loading,  machine  energy,  and  atomic  number  of 
the  material,  is  still  supported  by  the  enhanced  transport  calculations. 
This  is  important  because  it  provides  strong  theoretical  support  for  the 
belief  that  high  Z  experiments  performed  on  powerful  PRS  machines  such  as 
DECADE  will  exhibit  efficient  K-shell  x-ray  production.  Thus,  enhanced 
transport  calculations  can  provide  valuable  guidance  for  designing  future 
experiments  to  be  performed  on  these  high  power  machines. 
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The  turbulence  coefficients  used  in  the  HEART  and  BEAM’S  studies  were 
obtained  from  stagnation  density  and  temperature  comparisons  with  PI 
experiments  all  having  the  same  52  pg/cm  argon  load.  We  have  since  extended 
the  HEART  and  BEAM’S  work  by  examining  how  this  determination  of  turbulence 
coefficients  is  altered  as  the  mass  of  the  argon  load  used  in  the 
experiments  is  increased  substantially  above  52  pg/cm.  In  this  case  1-D 
simulations  are  performed  of  Double  Eagle  argon  experiments  having  mass 
loads  of  42,  85,  125  and  165  pg/cm.  In  each  simulation  the  transport 
coefficients  are  varied  until  a  reasonable  match  with  experiment  is  found. 
The  technique  (HEART  paper)  of  multiplying  the  classical  coefficients  by 
constant  multipliers  is  utilized;  however,  in  the  interest  of  limiting  the 
parameter  space  of  possible  coefficients  it  was  assumed,  somewhat 
arbitrarily,  that  the  mechanism  responsible  for  enhanced  viscosity  is  also 
responsible  for  an  equal  enhancement  of  ion  thermal  conduction.  The 
calculated  quantities  that  are  matched  to  the  experimental  data  are  the  mass 
averaged  temperature  and  density  of  the  K-shell  emitting  region.  The  latter 
quantities  were  inferred  from  the  measured  K-shell  spectrum,  pulsewidth,  and 
pinhole  picture.  The  MHD  model  is  the  same  as  that  described  in  the  HEART 
paper.  All  of  the  simulations  employ  a  circuit  model  of  Double  Eagle  so  that 
the  driving  current  and  the  plasma  load  dynamics  in  the  calculation  is  the 
same  as  the  machine  dynamics.  The  calculations  are  begun  at  the  same  radius 
as  the  experiments  with  93  percent  of  the  mass  concentrated  within  an 
annular  shell  having  the  identical  thickness  as  the  0.4  cm  nozzle  exit.  The 
remaining  mass  is  distributed  with  a  mass  density  that  exponentially  decays 
towards  the  axis.  This  initialization  procedure  provides  for  numerical 
stability  and  it  also  simulates  the  backpressure  effects  of  plasma  that  is 

blown  toward  the  axis  by  the  gas  jets. 

1-4 

A  major  premise  of  past  work  and  also  of  this  present  analysis  is 
that  thermalization  of  kinetic  energy  determines  the  state  of  the  stagnation 
plasma.  This  is  insured  in  the  calculations  by  terminating  the  current 
(current-off  calculation)  when  the  voltage  across  the  diode  reaches  3  MV. 
Because  there  is  no  current  beyond  a  certain  early  time  in  the  on-axis 
dynamics,  no  attempt  is  yet  made  to  model  ohmic  heating,  anomalous  heating, 
and  plasma  confinement  by  magnetic  fields  during  these  latter  stages  of  the 
implosion.  All  of  these  phenomena  will  tend  to  increase  the  K-shell  yield. 

The  HEART  study  showed  that  when  late  time  current  effects  are 
neglected,  enhanced  resistivity  does  not  play  as  strong  a  role  as  enhanced 
viscosity  and  enhanced  ion  and  electron  thermal  conduction  in  lowering 
stagnation  ion  densities.  In  this  regard,  we  are  going  to  take  the  point  of 
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view  that  prior  to  stagnation,  most  of  the  major  differences  between 
experiment  and  non-enhanced  transport  calculations  are  attributable  to 
macroscopic  turbulence.  Thus,  only  the  artificial  viscosity  v,  which  heats 
directly  the  ions,  and  ion  thermal  conductivity  are  enhanced  in  these 
calculations.  Our  reason  for  not  enhancing  the  electron  thermal  conduction 
is  discussed  in  more  detail  at  the  end  of  the  results  section. 

Note,  since  late  time  current  effects  are  neglected  and  also  because  it 
is  the  nature  of  a  1-D  model  to  neglect  zippering  and  other  large  scale 
multi-dimensional  macroscopic  motion  (all  of  which  can  significantly 
lengthen  the  time  scale  over  which  kinetic  energy  is  thermalized  and 
consequently  can  broaden  the  radiation  pulse)  it  is  unreasonable  to  expect 
exact  agreement  between  calculations  and  experiment.  As  mentioned  at  the 
beginning  of  the  introduction  we  are  primarily  seeking  a  reduction  in  the  2 
-  3  order  of  magnitude  differences  in  density  and  radiative  powers  that 
exists  between  1-D  non- turbulent  calculations  ard  experiments. 

RESULTS 

The  results  of  this  theoretical  and  experimental  comparison  are 
summarized  in  Tables  1  through  4  for  the  42,  85,  125,  and  165  pg/cm  argon 
loads,  respectively.  The  mass  of  each  plasma  load  was  determined  from  a  0-D 
slug  model  analysis  of  the  experimental  implosion  time.  The  K-shell  yields 
were  measured  using  kapton  filtered  tantalum  calorimeters  and  also  by 
numerical  integration  of  filtered  x-ray  diode  (XRD)  emission  power  signals. 
The  spatial  and  temporal  behavior  of  K-shell  emission  was  determined  using 
the  arrangement  of  Deeney  et.  al.^  Applying  the  analysis  technique  of 
Coulter  et.  al.6  to  the  time  integrated  and  resolved  K-shell  data  provided 
the  density,  mass  fraction,  and  electron  temperature  of  the  K-shell  emission 
region.  This  required  knowledge  of  the  outer  radius  of  the  K-shell,  the  K- 
shell  power,  and  the  He  to  Ly  line  ratio. 

As  in  the  HEART  and  BEAM’S  work,  it  was  possible  to  find  a  broad  range 
of  values  in  (x^,u)  parameter  space  that  produced  stagnation  temperatures 
and  ion  densities  in  acceptable  agreement  with  experiment.  Over  an 
appreciable  time  span  during  peak  K-shell  emission,  electron  temperatures 
could  be  matched  within  252  and  ion  densities  could  be  matched  within  a 
factor  of  3.  In  particular,  multiples  of  >10  to  25  times  the  classical 
conduction  and  viscosity  coefficients  for  the  42  pg/cm  load,  multiples  of  25 
to  40  for  the  85  pg/cm  and  125  pg/cm  load  and  multiples  of  >25  to  50  for  the 
165  pg/cm  load  all  produced  acceptable  agreement  with  experimental 
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stagnation  densities  and  electron  temperatures  in  the  K-shell  emission 
region.  Based  on  this  limited  data  set  it  appears  that  the  amount  of 
enhanced  transport  required  to  give  experimental  agreement  could  be 
increasing  by  a  factor  of  2  over  the  range  of  masses  from  42  to  165  pg/cm. 
However,  this  upward  variation  is  expected  to  saturate  with  increasing  mass 
because  enhancements  of  40  or  greater  over  classical  ion  conduction  are 
close  to  the  limit  of  physically  allowable  conduction.  To  insure  that 
calculations  do  not  go  beyond  the  bounds  of  physical  possibility,  flux 
correctors  will  be  utilized  in  future  models. 

The  most  representative  calculation  for  each  mass  load  produced  a  full- 
width-half-maximum  radiation  pulse  width  within  a  factor  of  two  or  three  of 
the  experiments.  This  is  good  agreement  considering  that  late  time  current 
behavior  and  axial  motion  such  as  zippering  are  not  modeled.  The  mass 
calculated  to  be  emitting  in  the  K-shell  as  well  as  the  outer  radius  of  the 

K-shell  emission  at  the  time  of  stagnation  are  in  excellent  agreement  with 

experiment.  These  results  indicate,  both  experimentally  and  theoretically, 
that  the  K-shell  emission  of  the  lower  mass  loads  are  dominated  by  bulk 

plasma  behavior,  whereas,  the  inner  core  plasma  does  the  K-shell  radiating 

for  the  larger  mass  loads.  Since  if  is  less  than  unity  for  the  more  massive 
loads,  they  could  not  be  bulk  K-shell  emitters. 

The  calculated  peak  K-shell  powers  agree  within  a  factor  of  1.5  with  the 

experiments  for  the  lower  mass  loads  and  within  a  factor  of  5  for  the  higher 

mass  loads.  Again,  this  is  good  agreement  considering  that  the  calculated 

and  experimental  ion  densities  only  agree  to  within  factors  of  2  -  3  and 

that  optically  thin  radiative  powers  are  very  sensitive  to  the  ion  density, 

2 

i.e.  radiative  power  a  , 

Since  K-shell  yields  are  dependent  upon  the  amount  of  magnetic  work  done 
on  the  plasma,  the  comparison  between  calculated  and  experimental  K-shell 
yields  has  large  uncertainties  at  this  stage  because  we  do  not  have 
sufficient  knowledge  of  the  magnetic  work  that  was  done  on  the  experimental 
plasmas.  As  soon  as  more  information  on  the  time  resolved  and  integrated 
minimum  outer  radius  of  L-shell  emission  is  made  available  to  us,  better 
determinations  of  this  magnetic  work  can  be  made.  Thus,  for  the  moment,  we 
find  that  the  calculated  K-shell  yields  for  the  more  massive  roads  are  in 
excellent  agreement  with  experiment,  whereas,  the  calculated  yields  for  the 
lower  mass  loads  are  substantially  smaller  than  the  corresponding 
experimental  yields.  It  is  likely  that  late  time  current  effects  are 
responsible  for  this  behavior.  One  reason  for  this  likelihood  is  provided  by 
the  observation  that  the  calculated  kinetic  energy  after  stagnation  is  •»  70 


21 


percent  of  the  original  magnetic  work  energy  coupled  to  the  plasma  for  the 
42  jug/cm  plasma  while  it  is  only  ~  10  percent  for  the  165  pg/cm  load.  Thus, 
the  lower  mass  loads  undergo  essentially  an  elastic  stagnation  in  which  very 
little  of  their  kinetic  energy  is  radiated  away  as  it  is  thennalized  in  the 
current-off  calculations.  Conversely,  most  of  this  energy  is  radiated  upon 
thermalization  for  the  more  massive  loads  even  though  only  a  small  fraction 
of  the  total  mass  is  needed  to  do  the  K-shell  radiating.  If  the  current 
continues  to  confine  the  plasma  after  stagnation  in  the  experiments,  then 
all  of  the  pre-stagnation  kinetic  energy  will  be  eventually  radiated  from 
the  low  mass  loads.  Since  the  value  of  tj  is  >  1.0  for  these  loads,  the  bulk 
of  the  plasma  is  hot  enough  that  a  significant  portion  of  this  radiation 
will  come  from  the  K-shell. 

Another  possible  late  time  current  effect  that  could  be  responsible  for 
the  experimental  K-shell  -ields  bei~$  larger  than  the  calculated  values  for 
the  low  mass  plasmas  is  ohmic  or  anomalous  heating.  If  the  plasma  is 
confined  on  axis  then  it  most  likely  will  be  in  Bennett  equilibrium  and  then 


B2/8ir  «  N.(l+Z)Te/(jrr2) 


(1) 


where  B  is  the  magnetic  field  strength,  r  is  the  stagnation  radius,  and  Z  is 
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the  effective  charge.  Substituting  B  «  2I/rc,  and  fT  -  m/(2Z*1.67  x  10  ) 

in  the  above  expression,  we  find 


I2  =  96  x  m  x  T 


(2) 


where  I  is  the  current  in  MA,  T  is  in  eV,  and  m  is  the  mass-per-unit- length 

e 

(g/cm).  In  order  to  produce  K-shell  emission  it  is  best  to  operate  near  the 

optimal  electron  temperature  required  to  effectively  radiate  in  the  K-shell, 
2  9  4 

i.e.  T  ~  0.3  Z  . 
opt 

Substituting  TQpt  into  expression  (2)  gives 


_2  _  _2.9 

I  ~  28  x  m  x  Z 


(3) 


For  Double  Eagle,  with  a  stagnation  current  of  3  MA,  and  Z  =  18  for  argon, 
the  optimal  mass  loading  to  achieve  a  K-shell  Bennett  equilibrium 
confinement  is  ~  73  p g/cm,  which  lies  between  the  42  and  85  Mg/cm  loads.  The 
heating  mechanism  responsible  for  maintaining  the  equilibrium  has  not  yet 
been  determined  and  is  under  investigation. 


The  mass  of  each  load  was  determined  using  a  0-D  slug  model  and 
knowledge  of  the  experimental  implosion  times  T,  which  is  the  time  to  the 
onset  of  the  K-shell  emission.  The  fact  that  a  0-D  model  does  not  include 
plasma  back  pressure  effects,  which  are  especially  important  for  those  loads 
which  start  close  to  the  axis,  could  explain  the  differences  between  the 
calculated  T  and  experimental  T  for  the  most  massive  loads.  If  this  is  true, 
then  125  jig/cm  and  165  pg/cm  are  larger  than  the  actual  respective  load 
masses. 

In  the  HEART  work,  viscosity,  electron  thermal  conduction,  and  ion 
thermal  conduction  were  enhanced  until  agreement  with  a  52  pg/cm  argon  load 
was  obtained;  this  occurred  for  multipliers  of  20  -  40.  Preliminary 
calculations  that  included  electron  thermal  conduction  also  produced 
experimental  agreement  for  the  42  pg/cm  and  85  pg/cm  loads.  However,  for  the 
more  massive  loads,  with  enhanced  electron  thermal  conductions  even  as  large 
as  a  factor  of  40,  the  stagnation  ion  densities  of  the  K-shell  emission 
region  were  still  nearly  two  orders  of  magnitude  larger  than  the 
experimental  densities.  This  is  the  primary  reason  enhanced  electron  thermal 
conduction  was  not  considered  in  this  work. 

SUMMARY  AND  CONCLUSIONS 
3  4 

In  two  earlier  papers,  *  a  series  of  1-D  MHD  calculations  were 
described  that  showed  how  the  imploded  mass  and  the  kinetic  energy  generated 
per-ion  influence  the  K-shell  emission  characteristics  of  a  z-pinch. 
However,  these  calculations  predicted  an  ideal  implosion  dynamics,  which  was 
not  seen  in  the  experiments.  Based  on  the  experimental  evidence  and  on  2-D 
MHD  z-pinch  calculations,  it  seems  likely  that  some  or  all  of  the  missing 
physics  in  the  1-D  calculations  is  caused  by  plasma  turbulence.  In  the  HEART 
and  BEAM’S  work  this  hypothesis  was  tested  by  employing  a  phenomenological 
approach  of  incorporating  turbulence  effects  into  the  1-D  model.  Multipliers 
were  introduced  for  the  plasma  (artificial)  viscosity,  the  ion  and  electron 
heat  conductivity,  and  the  electrical  resistivity.  It  was  then  possible  to 
make  an  appropriate  choice  for  these  multipliers  so  that  the  1-D  MHD 
calculations  produced  approximately  the  same  plasma  conditions  as  a  52  pg/cm 
PI  argon  gas  puff  experiment.  Again  this  provided  supporting  evidence  for 
the  presence  of  plasma  turbulence  in  z-pinch  experiments. 

In  this  paper,  a  similar  approach  as  outlined  above  was  used  in  order  to 
examine  phenomenologically  the  scaling  of  turbulence  effects  with  load  mass. 
This  was  accomplished  by  comparing  enhanced  transport  calculations  with  a 
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series  of  PI  argon  experiments  in  which  the  load  mass  was  varied.  The  masses 
in  these  experiments  ranged  from  42  to  165  pg/cm.  Preliminary  calculations 
had  showed  that  agreement  could  not  be  attained  with  the  experimental  plasma 
conditions  of  the  K-shell  emission  region  when  enhanced  electron  thermal 
conduction  was  included  in  the  calculations  for  the  larger  mass  loads.  For 
this  reason,  enhanced  electron  thermal  conduction  was  not  included  in  this 
work,  i.e.  macroscopic  turbulence  effects  were  modeled  only  by  enhanced 
viscosity  (artificial),  and  ion  thermal  conduction. 

We  have  found  that  it  is  possible  to  obtain  reasonable  agreement  with 
experimental  conditions  during  stagnation  of  the  z-pinch  over  a  range  of 
mass  loads  by  enhancing  the  viscosity  and  ion  thermal  conduction 
coefficients.  A  factor  of  >  10  to  25  enhancement  of  the  transport 
coefficients  was  required  in  order  to  get  the  42  pg/cm  load  to  approximate 
the  experimental  conditions  while  it  took  twice  this  factor  for  the  largest 
(165  fig/ cm)  load.  It  is  expected  that  the  amount  of  enhancement  required  to 
approximate  the  experimental  conditions  will  saturate  with  increased  mass, 
because  eventually,  a  flux  limit  for  physically  allowable  transport  will  be 
reached.  This  hypothesis  will  be  tested  by  including  ion  thermal  flux 
correctors  in  the  calculations  and  comparing  the  results  with  larger  mass 
experiments,  e.g.  Saturn  argon  experiments. 

Both  calculations  and  experiments  support  the  idea  that  the  bulk  of  the 
plasma  emits  in  the  K-shell  when  1 ]  >  1  and  that  it  is  only  the  plasma  core 
which  radiates  in  the  K-shell  when  JJ  <  1.  Previous  current-off  calculations 
which  included  enhanced  electron  thermal  conduction  did  not  show 
proportionately  as  much  K-shell  emission  from  the  core  of  f]  <  1  plasmas, 
most  likely  because  the  energy  of  the  core  is  transferred  too  rapidly  to  the 
bulk  of  the  plasma. 

The  stagnation  of  the  42  pg/cm  load  was  very  elastic,  i.e.  only  30 
percent  of  the  pre-stagnation  kinetic  energy  was  radiated  away  in  these 
current-off  calculations.  Given  that  the  current  is  still  present  during 
stagnation  in  the  experiment,  most  of  the  remaining  energy  (70  percent) 
would  be  confined  and  radiated.  This  is  a  possible  explanation  for  why  the 
K-shell  yield  from  the  lower  mass  plasmas  was  higher  in  the  experiments  than 
predicted  by  the  current-off  calculations.  We  showed  that  ohmic  heating  or 
anomalous  heating  during  Bennett  confinement  can  "cook"  the  lower  mass 
plasmas  while  they  exist  in  the  K-shell  ionization  stages. 

In  the  future  we  hope  to  investigate  the  scaling  of  turbulence 
coefficients  to  even  larger  masses  than  examined  here  and  also  to  different 
Z  elements.  If  a  phenomenological  determination  of  the  scaling  of  turbulence 


coefficients  with  Z  and  mass  can  be  made,  then  this  knowledge  will  be  used 
to  reexamine  the  BEAM'S  work  in  order  to  obtain  a  more  accurate 
determination  of  the  mass  boundaries  for  scaling  of  K-shell  emission 
efficiently  with  Z,  mass  loading,  and  kinetic  energy  (machine  energy).  As 
illustrated  in  this  study,  effects  such  as  plasma  confinement  and  heating 
during  stagnation  may  play  a  significant  role  in  influencing  the  K-shell 
yield.  By  continuing  our  strong  collaboration  with  the  experimentalists  the 
plasma  conditions  under  which  these  effects  are  important  can  be  determined 
by  comparing  the  experimental  results  with  calculations  in  which  the  full 
current  or  a  short-circuit  current  are  present  during  stagnation.  All  of 
this  work  has  substantial  implications  for  the  performance,  design,  and 
optimization  of  DECADE  PRS  loads. 
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Table  1.  Experimental  and  calculated  parameters  for  the  42  ug/cm  argon 

gas  puff  load,  is  the  time  duration  of  the  full-vidth-half-maximum 

of  the  K-shell  radiation  pulse.  t+  and  T-  are  the  respective  times  of  the 

right  and  left  hand  edges  of  H  is  the  magnetic  vork  done  on  the 

plasma  divided  by  the  product  of  the  ion  density  and  the  minimum  energy 

required  to  reach  the  K-shell  (E  .  ~  40  keV/ion). 

im 
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Argon  mass  is 

85  pg /cm  Experim 

Initial  radius  is  Results 
X.25  cm 


Experimental  IK.  and  v  x  1  IK,  and  v  x  25  K.  and  n  M 
l“i  —  Results 


Results 


K-shell  radiation 
pulse  width 

rFVHM  <nS) 


<Ni>  of  K-shell 
emission  region 
(10  cm*  ) 


<Te>  of  K-shell 
emission  region 
(keV) 


Hass  emitting  in 
the  K-shell 
(pg/cm) 


Outer  Radius  of 
K-shell  emission 
region  (cm)  at  T 


Minimum  outer 
radius  of  L- she 11 
emission  region 
(cm) 


K-shell  emission 
peak  power  (GV/cm) 


K-shell  yield 
(kJ/cm) 


Magnetic  work  done 
on  plasma  (kJ/cm) 


Ohmically  Coupled 
Energy  (kJ/cm) 
(kJ/cm) 


Implosion  time 
r  in  (ns) 


Current  at  t  (MA) 


Kinetic  energy 
after  stagnation 
t  «  •  (kJ/cm) 


Not 

Available 


Table  2.  Experimental  and  calculated  parameters  for  the  85  Mg/ cm  argon 
gas  puff  load.  is  the  time  duration  of  the  full-vidth-half-maximum 

of  the  K-shell  radiation  pulse.  T+  and  t~  are  the  respective  times  of  the 
right  and  left  hand  edges  of  h*  is  the  magnetic  work  done  on  the 

plasma  divided  by  the  product  of  the  ion  density  and  the  minimum  energy 
required  to  reach  the  K-shell  (Emjn  ~  40  keV/ion). 
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Argon  mass  is 

— 

125  rig/cm 

Experimental 

and  v  x  1 

K  and  V  x  25 
Results 

Kj  and  v  x  40 

Inital  radius  is 

Results 

Results 

Results 

1.0  cm 

K-shell  radiation 
pulse  width 

Not 

TFVHM  <ns) 

Available 

<  1 

~  3 

**  4 

<Ni>  of  K-shell 

T+ 

14. 

2.0 

1.2 

emission  region 
<10A*  cm"J ) 

Not 

T 

30. 

6.0 

4.0 

Available 

f- 

100. 

5.0 

3.6 

<Te>  of  K-shell 

T+ 

2.0 

2.0 

2.0 

emission  region 

1.6 

T 

2.0 

2.3 

2.0 

(keV) 

T- 

1.3 

1.4 

1.4 

Hass  emitting  in 
the  K-shell 

Not 

9. 

24. 

29. 

(jjg/cm) 

Available 

Outer  Radius  of 
K-shell  emission 

not 

0.02 

0.04 

0.07 

region  (cm)  at  r 

available 

Minimum  outer 

time  -  NA 

radius  of  L-shell 

integrated 

0.13 

0.14 

0.14 

emission  region 

time  -  NA 

(cm) 

resolved 

K-shell  emission 
peak  power  (GW/ cm) 

194 

2000 

800 

600 

K-shell  yield 

(kJ/cm) 

2.6 

1.4 

2.8 

2.33 

Magnetic  work  done 

Not 

on  plasma  (kJ/cm) 

Available 

10.5 

10.8 

11.2 

Ohmically 

Coupled  Energy 

Not 

1.1 

1.1 

1.0 

(kJ/cm) 

Available 

Not 

V 

Available 

0.88 

0.91 

0.93 

Implosion  time 
r  in  (ns) 

100 

107 

108 

108 

Current  at  T  (MA) 

3.2 

2.8 

2.8 

2.8 

Kinetic  energy 
after  stagnation 

Not 

0.6 

1.7 

2.3 

t  *  *  (kJ/cm) 

Available 

Table  3.  Experimental  and  calculated  parameters  for  the  125  ug/cm  argon 
gas  puff  load.  is  the  time  duration  of  the  full-vidth-half -maximum 

of  the  K-shell  radiation  pulse,  x  and  r~  are  the  respective  times  of  the 
right  and  left  hand  edges  of  h  is  the  magnetic  vork  done  on  the 

plasma  divided  by  the  product  of  the  ion  density  and  the  minimum  energy 
required  to  reach  the  K-shell  (Effljn  “  40  keV/ion). 
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Argon  mass  is 

165  fig/cm 

Experimental 

K.  and 

V1x  1 

K.  and 
t/Sc  25 

K.  and 
v  x  40 

K.  and 
v1x  50 

Initial  radius  is 

Results 

Results 

Results 

Results 

Results 

0.75  cm 

K-shell  radiation 
pulse  vidth 

Mot 

mm 

TFVHM  (nS) 

Available 

<  1 

-  2 

B 

~  4 

<Ni>  of  K-shell 

r+  io. 

6. 

3. 

2. 

emission  region 

2. 

r  50. 

12. 

6. 

5. 

(10  cm  ) 

T-  200. 

20. 

7. 

5. 

<Te>  of  K-shell 

QS1S| 

1.8 

1.8 

emission  region 

2.0 

1.8 

(keV) 

1 

3afl 

1.2 

1.3 

B 

Mass  emitting  in 
the  K-shell 

10. 

5. 

18. 

16. 

18. 

(f/g /cm) 

Outer  Radius  of 
K-shell  emission 
region  (cm)  at  T 

0.05 

.007 

0.02 

0.03 

0.04 

Minimum  outer 

time  -  MA 

radius  of  L- she 11 

integrated 

0.12 

0.13 

0.12 

0.13 

emission  region 

time  -  HA 

(cm) 

resolved 

K-shell  emission 
peak  power  (GW/ cm) 

140 

500 

850 

700 

630 

K-shell  yield 

(kJ/cm) 

1.7 

0.33 

2.2 

2.3 

2.3 

Magnetic  work  done 

Hot 

BB 

on  plasma  (kJ/cm) 

Available 

a. 7 

8.8 

BB 1 

9.2 

Ohmic ally 

Coupled  Energy 

Hot 

■ 

■f 

Ea 

B 

(kJ/cm) 

Available 

b 

B 

Hfl 

B 

Hot 

V 

Available 

0.55 

0.55 

0.57 

0.58 

Implosion  time 
r  in  (ns) 

95 

100 

100 

100 

100 

Current  at  T  (MA) 

2.9 

2.7 

2.7 

2.6 

Kinetic  energy 
after  stagnation 

Not 

0.8 

1.0 

t  »  •  (kJ/cm) 

Available 

B 

■1 

Table  4.  Experimental  and  calculated  parameters  for  the  165  Vg/cm  argon 
gas  puff  load,  Is  the  tiae  duration  of  the  full-vidth-half-maxiaum 

of  the  K-shell  radiation  pulse.  t+  and  t~  are  the  respective  tines  of  the 
right  and  left  hand  edges  of  h  Is  the  oagnetic  vork  done  on  the 

plasaa  divided  by  the  product  of  the  ion  density  and  the  minimum  energy 
required  to  reach  the  K-shell  (Enin  -  40  keV/ion). 


ill.  NON-MAX  WELLIAN  DISTRIBUTION  FUNCTIONS  IN  DECADE-CLASS  Z-PINCHES 


Introduction 

Non-Maxwellian  electron  distribution  functions  could  become  a  critical  factor  determining 
how  an  optimal  PRS  implosion  will  perform  at  DECADE  current  levels  and  beyond.  As  the 
mobile  component  of  the  plasma,  electrons  provide  most  of  the  thermal  pressure  and  carry  most 
of  the  current.  Varying  the  largely  electronic  thermal  and  electrical  conductivity  will  change  a 
“hard”  implosion  into  a  “soft”  implosion,  or  vice-versa.  Ohmic  plasma  heating  is  mediated  by 
the  electrons,  which  absorb  electrical  energy  from  the  generator  and  transfer  it  to  ions  by  both 
elastic  and  inelastic  collisions.  The  atomic  line  radiation  consequent  from  inelastic  electron-ion 
collisions  is  an  important  plasma  diagnostic,  and  could  play  a  crucial  role  in  the  energy  balance  that 
determines  whether  there  is  radiative  collapse.  Also  of  concern,  observations  of  beams  of  runaway 
electrons  under  certain  pinch  conditions  show  that  there  is  sometimes  energy  transfer  that  does  not 
result  in  plasma  heating;  these  runaway  electrons  may  also  be  associated  with  hot  spots,  poorly 
understood  regions  of  enhanced  radiation  and  probable  hydrodynamic  instability.  For  all  these 
reasons,  a  good  understanding  of  the  factors  influencing  the  electron  distribution  function  is  an 
important  part  of  understanding  present  PRS  experiments  and  of  designing  future  more  powerful 
PRS  devices. 

In  DECADE,  the  peak  current  will  be  many  times  that  of  earlier  machines.  This  increased 
current  means  increased  electron  flow  velocities  and  greater  electric-field  heating  of  the  electrons, 
which  means  more  favorable  conditions  for  ion-acoustic  microturbulence.  The  higher  plasma 
temperatures  in  DECADE  could  mean  increased  importance  for  electron-ion  excitation  collisions, 
a  factor  that  in  lower  power  machines  made  a  difference  of  only  a  few  percent.  (The  importance 
of  these  inelastic  collisions  is  proportional  to  temperature  for  fixed  ion  populations,  although  the 
population  levels  also  depend  on  temperature).  Higher  inductive  current  flows  in  DECADE  also 
mean  larger  electric  fields,  both  parallel  to  and  transverse  to  the  current,  which  could  change  the 
way  ohmic  heating  affects  the  distribution  function  from  lower-power  machines. 

Generally,  it  is  assumed  that  the  electron  distribution  has  a  Maxwellian  shape,  because  electron- 
electron  collisions,  which  are  always  present,  cause  the  distribution  to  relax  towards  a  Maxwellian 
equilibrium.  PRS  timescales  are  nanoseconds  to  hundreds  of  nanoseconds.  Electron-electron 
collisions  equilibrate  the  system  on  the  collisional  timescale: 

Tee  =  0.2588-^— ps,  (1) 

mg  log  A 

where  kT  is  the  temperature  in  eV,  nig  is  the  density  in  units  of  1018  cm-3,  and 
log  A  =  log  [5.1596(Kr)3/2/mg2]  is  the  Coulomb  logarithm.  At  low  energies,  below  the  thermal 
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energy,  electron-electron  collisions  are  increasingly  dominant,  but  equilibration  is  slower  for  the 
high-energy  electrons  in  the  tail  of  the  distribution.  In  the  distribution  tail,  the  effects  of  ohmic 
heating  and  collisional  cooling  to  produce  non-Maxwellian  distributions  are  quite  pronounced, 
especially  at  high  temperatures,  and  are  strongly  influenced  by  the  strength  of  the  magnetic  field. 
In  addition,  the  effects  of  these  processes  can  be  magnified  by  plasma  microturbulence,  which  if 
present  can  also  greatly  increase  the  resistivity  and  viscosity. 

In  this  report,  we  lay  out  a  general  analysis  of  how  the  electron  distribution  function  develops 
in  a  PRS  implosion.  These  results  form  a  comprehensive  basis  for  finding  plasma  conditions 
(and  implosion  parameters)  that  create  non-Maxwellian  distributions.  In  the  first  section,  we 
describe  the  plasma  conditions  and  the  methods  of  study.  The  second  section  sketches  the  relevent 
kinetic  theory,  which  includes  electron-electron  collisions,  a  model  of  ohmic  heating,  and  inelastic 
collisions  with  ions.1  We  then  discuss  ion-acoustic  turbulence,  whose  properties  and  likely 
presence  in  implosions  have  been  presented  earlier.2  We  show  how  ion-acoustic  turbulence  can 
be  simply  included  in  our  model,  and  introduce  a  parameter  to  estimate  its  strength. 

The  consequences  of  this  kinetic  theory  are  explored  in  section  4.  The  actual  shape  of 
the  distribution  function  depends  on  the  interplay  between  the  Maxwellian-producing  electron- 
electron  collisions  and  the  various  other  processes,  whose  effects  depend  primarily  on  temperature, 
current  density  and  magnetic  field  strength.  Ohmic  heating  in  a  strong  magnetic  field  produces 
depleted-tail  distributions  of  a  well-established  form,*  but  for  moderate  magnetic  field  strengths 
or  in  the  presence  of  ion-acoustic  turbulence  more  complicated  distribution  function  shapes  are 
produced.  Inelastic  electron-ion  collisions,  surprisingly,  tend  to  produce  relatively  enhanced 
distribution  function  tails.  We  estimate  the  relative  importance  of  these  processes,  depending  on 
plasma  conditions  and  on  the  strength  of  turbulence. 

To  obtain  practical  results  from  the  model,  we  need  a  way  to  determine  the  local  electric  field 
in  the  plasma.  In  fact,  all  that  is  generally  available  is  an  estimate  of  the  total  current  through 
the  pinch,  from  which  we  can  estimate  the  local  current  density,  based  on  assumptions  about  the 
radial  current  profile.  The  current  density  is  directly  related  to  the  current-aligned  component  of 
the  electric  field  by  the  transverse  Spitzer  resistivity.4  In  addition  to  this,  a  large  radial  component 
of  the  electric  field  is  needed  to  drive  the  current  across  the  pinching  magnetic  field.  It  is  the  total 
electric  field,  not  just  the  current-aligned  component,  that  is  used  in  our  model.  Therefore,  we 
show  in  section  5  how  to  calculate  the  resistivity  in  terms  of  the  actual  distribution  functions;  this 
allows  us  to  find  the  electric  field  in  terms  of  the  current  density.  Simplified  analytic  expressions 
for  the  conductivity  are  introduced  in  this  section,  which  are  useful  for  showing  the  basic  scaling 
properties  of  the  resistivity.  Numerical  calculations  of  the  resistivity  of  self-similar  distributions8 
are  also  presented,  showing  non-Maxwellian  effects  for  various  magnetic  fields.  The  increase  of 
resistivity  with  turbulence  is  shown. 
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Runaway  phenomena  are  briefly  discussed  in  section  6.  In  contrast  to  tokamacs,  where  current 
is  mostly  along  magnetic  field  lines  and  runaways  are  commonly  observed,  PRS  devices  have 
transverse  current  and  magnetic  field,  and  freely-accelerated  runaways  tend  to  be  suppressed. 
Runaways  are  possible  where  the  electric  field  exceeds  the  magnetic  field  strength;  where  and  how 
this  situation  exists  depends  sensitively  on  how  current  penetrates  into  the  conducting  plasma,  a 
subject  that  will  be  studied  in  detail  in  the  near  future.  Here,  we  discuss  conditions  for  runaway 
generation  and  show  how  runaways  would  be  manifested  in  our  model.  In  particular,  the  increasing 
anisotropy  of  the  electron  distribution  with  increasing  current  is  pointed  out. 

Finally,  we  illustrate  the  model  developed  here  by  calculating  in  section  7  the  time-evolution 
of  the  electron  distribution  function  in  a  selenium  z-pinch  at  different  temperatures  and  magnetic 
field  strengths.  An  energy  equilibrium  is  assumed,  where  ohmic  heating  is  balanced  by  inelastic 
collisions  and  radiative  dissipation.  The  K-  and  L-shell  inelastic  collision  strengths  for  selenium 
were  derived  from  radiative  rates  determined  by  the  Radiation  Hydrodynamics  Branch,5  which  are 
also  used  in  Branch  MHD  calculations.  Ion  populations  were  determined  by  collisional-radiative 
equilibrium  (CRE)  calculations  for  the  temperature  and  density. 

These  results  represent  only  the  initial  step  in  applying  the  kinetic  model  developed  here  to 
study  PRS  electron  distributions.  The  model  includes  the  important  influences  on  the  distribution, 
even  nonlinear  turbulence  phenomena.  When  this  is  coupled  with  the  comprehensive  atomic 
database  of  the  Radiation  Hydrodynamics  Branch,  as  we  have  done,  model  predictions  can  be 
used  to  study  the  effect  of  non-Maxwellian  distributions  on  energy  balance  within  the  device. 
The  transport  coefficients  (e.g.,  resisistivity)  from  the  model  can  be  used  to  improve  MHD 
modeling.  The  next  step  will  be  to  determine  for  what  parameters  non-Maxwellian  corrections 
are  significant;  for  this,  a  self-consistent  scheme  for  determining  the  current  and  field  distribution 
must  be  developed.  Plans  for  this  step  are  discussed  in  the  conclusion. 

1.  Theory 

We  are  looking  at  the  evolution  of  the  electron  distribution  function  in  a  high  atomic  number 
plasma  in  an  imploding  z-pinch.  The  ions  are  much  heavier  than  the  electrons,  and  reach 
equilibrium  much  faster,  and  so  we  assume  the  ion  distribution  to  be  a  Maxwellian.  Electron- 
electron  collisions  are  very  important  in  energy  exchange  and  in  determining  the  isotropic  part  of 
the  distribution,  but  are  not  so  important  in  exchange  of  momentum;  therefore  we  ignore  the  direct 
contribution  of  electron-electron  collisions  to  the  current  (that  is,  we  assume  a  Lorentzian  plasma). 
Electron-ion  collisions  are  the  dominant  factor  determining  the  current;  these  are  much  more 
frequent  than  electron-electron  collisions,  so  that  the  directional  part  of  the  electron  distribution  is 
in  quasi-equilibrium.  Inelastic  collisions  affect  only  the  isotropic  part  of  the  distribution;  they  are 
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not  frequent  enough  to  contribute  on  the  electron-ion  timescale. 

The  system  is  composed  of  electrons  and  various  ion  species  a,  each  with  charge  Zae ,  mass 
ma  and  density  na.  The  mean  ion  charge  number  is  Z  =  ne/rii,  where  ne  =  ]T)a  naZa  and 
n<  =  na\  we  assume  Z  »  1.  The  plasma  frequency  for  species  a  is  u)a  =  V^zl  e2nQ/mQ. 
The  electric  field  E  and  magnetic  field  B  are  constant  and  uniform;  the  electrostatic  force  on 
electrons  is  ma  —  eE,  and  the  vector  cyclotron  frequency  is  uc  =  ( e/mc)B .  Spatial  variation  of 
all  quantities  is  neglected  in  the  model. 

The  ions,  with  mass  ma  much  greater  than  the  electron  mass  m,  are  assumed  to  maintain  a 
Max  well -Boltzmann  distribution  fa(v).  The  electrons  are  described  by  a  distribution  function 
f{v,t),  normalized  to  ne.  Here,  a  Cartesian  tensor  expansion  is  used®  with  only  the  first-order 
term  retained,  so  that 

f(v,t)tz  f0(v)  +  v-fi{v).  (2) 

Thus,  for  example,  the  average  energy  (E)  is  determined  by  fo(v): 

(E)  =  —  f  ~mv2  f0(v)v2dv  (3) 

J o  2 

and  the  current  density  J  is  determined  by  fi(v): 

3  =  ^TeJo  {enev)fi(v)v2dv.  (4) 

In  general,  the  Cartesian-tensor  expansion  has  an  infinite  number  of  terms,  with  term  fn(y)  being 
an  nth  order  tensor.  The  magnitude  of  /i(t>)  determines  the  directionality  of  the  distribution:  The 
distribution  is  isotropic  (though  not  necessarily  Maxwellian)  if  fi(v)  ~  0,  and  the  greatest  possible 
anisotropy  is  if  f\{y)  =  fo(v).  Except  in  special  circumstances,  neglect  of  the  higher-order  terms 
implies  that  fi(v)  -C  fo(v). 

The  “temperature”  of  the  nonequilibrium  system  is  defined  to  be  proportional  to  the  average 
energy  (in  units  of  eV,  when  in  numerical  formulas): 

kT  =  |(B),  (5) 

and  the  thermal  velocity  is  defined  as: 


The  (Rutherford)  collision  frequency  of  a  particle  of  species  a  as  it  moves  with  velocity  va 
through  a  number  density  np  of  species  0  particles  is: 

=  Ya0~Y  (V 

vl 
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where 


4tt  ZlZle* 

Y<*0  —  2  ^°6  Aap .  (8) 

771  a 

and  log  A  is  the  Coulomb  logarithm.7  The  characteristic  frequency  in  this  system  is  vr  =  uee(vth)- 
The  inelastic  collisions  considered  here  transfer  energy  from  the  electrons  to  the  ions,  exciting 
the  ion  to  a  higher  energy  level.  De-excitation  collisions,  which  transfer  energy  from  ions  to 
electrons,  are  neglected.  The  transition  energy  from  level  a  to  level  6  is  e0f,  below  this  electron 
energy  there  is  no  energy  transfer.  The  actual  inelastic  collision  frequencies  are  complicated 
functions  of  atomic  parameters,  but  they  can  be  conveniently  expressed  in  terms  of  slowly- varying 
semi-empirical  collision  strengths,  of  order  one.  The  collision  frequency  for  an  electron  of  velocity 
v  to  excite  an  ion  from  level  a  to  level  6  is 

Vab  =  ««  |*ag  — 8(e  -  Cat)  V,  (9> 

where  e  =  \mv2,  Ry  =  13.6  eV  is  the  Rydberg  energy,  oq  is  the  Bohr  radius,  ga  is  the  multiplicity 
of  level  a,  and  £2ab(e)  is  the  collision  strength  mentioned  above. 


~£afc)  v, 


2.  Basic  description:  without  turbulence 

The  electron  distribution  function  evolves  according  to  the  Boltzmann  equation: 

df 

—  +  (a  +  V  X  U>e)  •  Vvf  =  Cel  +  Cinel  (10) 

with  small-angle  collisions  Cei  described  by  a  Landau  collision  term,8  and  inelastic  collisions 

Cinei  included  as  discussed  below.  Using  the  first-order  Cartesian-tensor  expansion,  assuming 

the  ions  to  be  in  a  spherically-symmetric  Maxwellian  distribution  and  neglecting  the  effect  of 

— • 

electron-electron  collisions  on  /i(w)  (i.e.,  assuming  a  Lorentzian  plasma),  we  obtain  two  coupled 
kinetic  equations: 

-J-  =  Ceb  +  Cfp  +  C*  (11) 

^  +  «' X  fi  =  -VciMfi  •  (12) 

The  isotropic  distribution  fo(v)  is  determined  by  the  electric-field  heating  term  Ceb ,  a  Fokker- 
Planck  elastic  collision  term  Cfp  and  an  inelastic  collision  term  C*.  The  electric-field  heating 
term  is: 

=  (13) 
This  is  evidently  proportional  toj-E  (cf.  eq.  (4)).  The  elastic  collision  term  describes  the  effects 
of  multiple  small-angle  electron-electron  collisions: 

[/» l  /•(•>“  dv'+rt{hl  Mv>*  dv'+v[  Mv>' iv')  ■ 

(14) 
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The  inelastic  term  is  of  the  form 


C* 


^[-^(e)/o(e)  +  \j  j  v*{e  +  £j)fo(e  +  e,)j 


(15) 


where  fj(e)  is  the  energy-deoendent  collision  frequency  (given  in  eq.  (9))  and  ej  is  the  threshold 

energy,  for  the  jth  inelastic  collision  process.  This  term  cont"'~,'es  the  number  of  ions  and 

free  electrons,  as  it  should,  without  ionization  or  recombination:  An  inelastic  collision  simply 

transfers  electrons  within  the  distribution  from  energy  t  to  energy  e  -  £j.  The  factor  y'(e  +  £j)/e 

makes  allowance  for  the  greater  phase-space  volume  available  at  energy  e  +  e,:  the  number  of 

particles  scattered  between  energy  e  and  e  +  tj  is  always  u*(e)f(e)4ir%/e  de.  The  contribution 

of  inelastic  collisions  to  the  fi  equation  is  neglected  because  v*  <  (u* /vri  <x  Z~2),  and  so 

■  ■< 

inelastic-collision  changes  to  f\  are  relatively  slow. 

The  inelastic  collision  term  C*  is  related  to  the  excitation  rates  commonly  used  in  MHD  codes. 
The  excitation  rate  from  state  a  to  state  b  is  defined  such  that  dnb/dt  =  noWtLb\  it  is  the 
inelastic  collision  frequency  ulb  times  the  number  of  electrons  available  in  state  a  to  excite  the 
transition.  Thus,  the  rate  depends  on  the  electron  distribution  function  assumed  for  its  computation 
(usually  a  Maxwellian): 

w.t  =  —  r  vlh(v)f,ivWdv.  (16) 

n*  Jo 

All  information  about  the  collision  is  found  in  the  collision  strength  £20fc(e).  Although  the  rates 
are  computed  as  shown  from  the  collision  strengths,  it  is  not  feasible  to  deduce  collision  strengths 
from  rates.*  We  have,  however,  successfully  reconstructed  collision  strengths  from  rates  by  fitting 
to  the  assumed  functional  form  Clab  —  Sloe4**-  This  self-consistent  procedure  entails  only  a  small 
error  in  computing  radiative  output  (determined  by  recalculating  rates  from  the  deduced  collision 
strengths),  and  will  give  good  estimates  of  C*  for  most  collision  strengths. 

As  a  step  towards  solving  the  coupled  eqs.  (11-12),  we  now  find  an  approximate  solution 
to  eq.  (12),  thus  determining  Ceb  and  leaving  only  one  equation  to  solve.  We  note  that  /i 
equilibrates  at  the  high  frequency  of  */„•  =  Zuee,  and  so  assume  it  to  be  in  quasi-static  equilibrium 
with  /o.  The  equilibrium  solution  to  eq.  (12)  is  in  general 


wc  x  a - a Jc  —  ueid 


dfo 

dv 


(17) 


Note  that  when  a  ||  wc,  this  is  independent  of  wc,  as  it  should  be.  Here,  we  consider  the 

*  It  is  in  principle  possible  to  invert  the  numerical  integration  sums,  which  form  a  product  of  a 
Vandermonde  matrix  with  a  vector  of  collision  strength  values  on  a  grid.  In  practice,  however,  the 
matrix  is  much  too  ill-conditioned. 
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complementary  situation,  a  _L  wc,  where 

r  _  JLA 

°EB  3u2  &v 

In  summary,  the  kinetic  equation  without  turbulence  is  just  eq.  (11),  with  eq.  (18)  for  Ceb . 
eq.  (14)  for  Cfp  and  eq.  (15)  for  C*. 


dfo] 


^(l  +w2/vM  dv 


(18) 


3.  Inclusion  of  turbulence 


The  non-turbulent  plasma  is  well  defined  in  terms  of  slowly-varying  average  quantities  like 
the  density.  Although  there  are  always  microscopic  fluctuations  from  the  average  values,  these  are 
small  in  magnitude.  But  in  turbulent  plasmas,  the  fluctuations  are  not  small,  and  can  greatly  affect 
plasma  dynamics.  Turbulent  fluctuations  generally  involve  collective  modes  of  the  plasma  and  are 
generated  by  a  plasma  instability;  at  a  minimum,  the  turbulent  mode  cannot  be  strongly  damped. 
We  have  conducted  a  survey2  of  important  instabilities  for  typical  z-pinch  parameters  and  found 
that  the  most  likely  important  sources  of  turbulence  are  the  ion-acoustic  and  lower-hybrid  drift 
instabilities.  Here  we  look  exclusively  at  the  ion-acoustic  instability,  which  we  found  to  be  favored 
at  certain  times  in  a  variety  of  high-Z  z-pinches.  Although  the  lower-hybrid  drift  instability  might 
be  important,  it  may  not  have  sufficient  time  to  become  established  in  these  short  implosions. 

The  ion-acoustic  dispersion  relation  (real  part)  is 


“(*)  = 


kc, 

y/i+k*\*D* 


(19) 


where  c,  =  y/ZKT€jmi  is  the  ion-sound  velocity  and  Xp  =  y/KTe/4imee2  is  the  electron  Debye 
length.  The  growth  or  damping  rate  is  given  by  the  imaginary  part  of  w,  and  for  a  mode  to 
be  important  its  growth  time  r  =  $rm(t i;)-1  must  be  positive,  and  significantly  shorter  than  the 
implosion  timescale.  In  a  stationary  plasma  and  when  electron  and  ion  temperatures  are  equal,  the 
ion-acoustic  mode  is  Landau-damped.  Instability  occurs  only  when 


{ 


ZTe/Ti  >  1 

Uo  >  c , 


(20) 


where  the  electron  fluid  velocity  is 


u0  =  i/eri'. 


(21) 


Thus,  for  ion-acoustic  turbulence  to  be  significant  in  an  implosion,  there  must  be  both  a 
significant  electron-ion  thermal  disequilibrium  and  a  large  flow  velocity  (high  current  or  low 
density).  Electron-ion  equilibrium  depends  on  details  of  the  plasma  heating:  Electric-field  heating 
gives  energy  primarily  to  the  electrons,  which  then  collisionally  transfer  this  energy  to  ions.  The 
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equilibration  timescale  is  v"-1,  which  we  have  already  argued  is  very  short.  When  there  is  strong 
heating,  however,  such  as  early  in  runin  and  at  stagnation,  the  electron  temperature  can  exceed 
ion  temperature.  Even  a  small  such  disequilibrium  can  trigger  instability  if  the  ions  are  strongly 
ionized  ( Z  »  1).  As  the  electrons  heat  up  relative  to  the  ions,  they  interact  with  the  current  flow 
to  trigger  the  instability.  Without  ion  damping  (an  extreme  limit),  the  maximum  growth  rate  is 
such  that  u?jTmaz  «  3.19vt/i/«o.  or 


«  Oft  ,  ft_5  y/ fli7lK.T / Z 

rmat  ~  2.30  x  10  — - : - ns, 

J 


(22) 


where  m  =  mi/mpr^ton  is  the  mass  number  of  the  ions.  This  could  easily  translate  to  a  growth 
timescale  of  nanoseconds,  the  characteristic  implosion  time. 

After  triggering,  the  instability  develops  nonlinearly,  and  will  persist  even  at  saturation,  when 
the  conditions  drop  below  triggering  levels  and  the  growth  rate  is  zero.  Energy  from  the  current 
flow  continues  to  go  into  fluctuations,  as  reflected  in  the  anomalously  high  resistivity  and  heating 
levels  observed.9'10 

We  now  show  how  to  use  quasiiinear  theory  to  incorporate  turbulence  into  our  kinetic  model. 
We  will  derive  a  parameter  for  the  turbulence  strength  and  use  it  to  estimate  the  impact  of 
fully-developed  ion-acoustic  microturbulence.  In  quasiiinear  theory,  the  distribution  function  is 
written  as  the  sum  of  a  slowly-varying  average  (f(v))  (identical  to  f(v,t)  of  eq.  (2))  and  a 
rapidly-varying  fluctuation  6f(v,t)  (introduced  here).  Other  quantities  (e.g.,  fields)  are  similarly 
analyzed.  If  this  expansion  is  used  in  eq.  (10),  the  average  and  fluctuating  distribution  satisfy  the 
“linearized”  equations:3 


£  +  (a  +  vx  u?e).V„ 
^  +  (e  +  vxwc)-V, 


f(v,t)  =  -{6a*Vv6f) 
6f(v,t)  =  -ffi.  V*/(v) 


(23) 

(24) 


Eq.  (23)  is  just  the  Boltzmann  equation,  eq.  (10),  omitting  the  collision  terms  but  adding  a 
“quasiiinear”  source  term  on  the  right-hand  side.  The  difference  of  two  similar  second-order  terms 
has  been  neglected  in  eq.  (24),  on  the  argument  that  their  influence  on  the  average  distribution  is 
effectively  third  order.  The  RHS  of  eq.  (23)  provides  a  collision-like  turbulence  term  for  eq.  (10). 
Evaluating  it  requires  solving  for  both  the  fluctuating  distribution  and  the  consequent  fluctuating 
fields  from  eq.  (24)  and  Maxwell’s  equations.  Making  the  fluctuating  densities  and  fluctuating 
fields  self-consistent  results  in  a  dispersion  equation,  and  the  undamped  solutions  of  the  dispersion 
equation  define  the  turbulent  modes. 

Eq.  (24)  is  in  the  form  of  the  linearized  Vlasov  equation  and  can  be  solved  in  the  same 
way.  Fortunately,  the  magnetic  field  can  be  ignored  in  the  solution,  as  long  as  it  doesn’t 
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significantly  affect  the  predicted  dispersion  relation.  Therefore,  we  can  consistently  study  plasma 
microturbulence  from  electrostatic  waves  while  retaining  the  dynamical  effects  of  the  magnetic 
field  in  eq.  (10).  Thus,  Fourier  transforming  to  solve  eq.  (24)  and  using  Poisson’s  equation  to 
obtain  the  fluctuating  fields  from  the  fluctuating  densities,  the  new,  quasilinear  term  has  a  diffusive 
form: 

—{S3.  •  Vv6f)  =  V„  •  [l?  •  Y„/j  (25) 

where  the  diffusion  tensor  is 


T>a  =  d’k^w(t)g(^(k)  -k-v) 


(26) 


and  W(k)  is  the  electrostatic  wave  energy  density,  defined  such  that  Jd3k  W(k)  =  ( 6E2)/A-n . 

The  specific  form  of  the  quasilinear  diffusion  tensor  V  depends  on  the  turbulent  mode  w(fc) 
and  energy  density  W(k).  If  we  use  the  real  part  of  the  ion-acoustic  dispersion  relation,  eq.  (19), 
to  solve  for  V,  and  then  perform  a  Cartesian-tensor  expansion,  we  obtain  quasilinear  terms  in  both 
the  scalar  and  vector  equations  (eq.  (11)  and  eq.  (12)): 


■[v-V.f 

(0) 

4ir2e2£l 

m2 

■[v.v.f 

(1) 

47r2e2£2 
m2v 2 

(v  >  c.) 


£1  AW 

v2  dv  \t>  dv  ) 
e?A 


dv 


(27) 

(28) 


Both  terms  are  zero  for  very  low  velocities,  when  v  <  c,.  The  parameter  £2,  a  measure  of  the 
turbulence  energy,  is  the  average  fluctuation  energy  density  per  wavelength: 


Q  =  4tt 


J  kW(k)dk. 


(29) 


With  the  turbulence  term  of  eq.  (28)  added,  eq.  (12)  becomes  a  vector  differential  equation 
in  t>  which  is  difficult  to  solve  even  with  our  assumption  of  equilibrium  /a.  Fortunately,  it 
is  a  good  assumption  to  neglect  the  differential  terms  in  both  eq.  (28)  and  eq.  (27),  since 
cl/vth  —  Zme/mi  <  1.  Thus: 


where  the  turbulence  strength  parameter  is 

a  —  47t'2e2Q 

-  Zm2vlhvR' 


(32) 
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Thus,  ion-acoustic  turbulence  just  acts  to  multiply  the  electron-ion  collision  frequency  by  a  factor 
/9.  The  magnitude  of  can  be  roughly  estimated  by  noting  that,  if  l  is  the  dominant  length  scale 
of  the  fluctuations,  then  ft  «  i  x  [Fluctuation  energy  density],  and  so 


& 


u. 


Zvr  A  d 


} E , 


(33) 


where  Je  is  the  fraction  of  thermal  energy  contained  in  the  fluctuations. 

The  fluctuation  energy  Je  is  generally  proportional  to  N  =  n\*D,  the  number  of  particles 
in  a  Debye  sphere;  for  most  plasmas,  N  »  1.  The  exact  relationship  depends  on  the  degree  of 
turbulence  in  the  plasma.  In  a  non-turbulent  plasma,  the  well-known  result  is  that  Je  ~  TV-1;  for 
moderate,  fully-developed  turbulence,  where  the  instability  growth  has  saturated,  Je  ~ 
and  for  strong  turbulence,  Je  ~  JV°,  and  the  fluctuation  energy  is  comparable  to  the  thermal 
energy.11  We  account  for  this  by  writing 


Je  -  N~*  (34) 

where  0  <  7  <  1;  for  saturated  microturbulence,  7  as  Using  the  relationship 

u}e/vR  =  2-rrN/log  A,  where  2a:/  log  A  is  of  order  one,  the  turbulent  enhancement  in  collision 
frequency  is 

?  ~  <3S) 

Table  1  shows  the  value  of  N  for  various  plasma  conditions.  For  ion-acoustic  turbulence,  it  is 
commonly  assumed  that  l  fa  Xd-  There  are  a  wide  range  of  possible  enhancements  to  the  collision 
frequency  due  to  ion-acoustic  microturbulence;  for  typical  PRS  conditions,  factors  of  10-1 00  are 
not  unreasonable. 


Table  1.  Debye  length  (nm),  N  =  particles  in  Debye  Sphere. 


^  0 

§ ! 
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10 

Temperature  (< 
100 

jV) 

500 

1000 

1016 

235 

N=130 

743 

N=4,100 

1,660 

N=45,900 

2,350 

N= 130,000 

101* 

23.5 

N=13 

74.3 

N=410 

166 

N=4,590 

235 

N=  13,000 

1020 

2.35 

N=1.30 

7.43 

N=41 

16.6 

N=459 

23.5 

N=  1,300 

1022 

0.24 

N=0.13 

0.74 

N=4.10 

1.66 

N=45.9 

2.35 

N=130 

The  increase  in  the  effective  collision  frequency  is  enough  in  itself  to  significantly  increase  the 
resistivity  of  the  plasma.  In  addition,  the  turbulence  changes  the  shape  of  the  distribution  function 
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in  a  way  that  also  increases  the  resistivity.  The  modification  to  the  kinetic  equation,  eq.  ( 1 1 ),  due 
to  turbulence  is  found  by  including  the  turbulence  term  of  eq.  (31)  in  the  solution  to  eq.  (12). 
The  solution  for  Ceb  remains  in  the  form  of  eq.  (18),  but  instead  of  the  effective  collision 
frequency 

vlif  =  (!  +  0)uei  (36) 

must  be  used.  If  the  magnetic  field  is  low,  this  increase  in  collision  frequency  changes  the 
magnitude  of  the  ohmic  heating  term;  if  the  magnetic  field  is  high,  the  higher  collision  frequency 
changes  primarily  the  nature  of  the  heating. 

4.  Characteristics  of  Non-Maxwellian  Terms 

For  computation  it  is  useful  to  introduce  dimensionless  units,  such  that 

T  =  URt 

e  =  mv2 /2kT. 

We  define  the  turbulence-modified  effective  collision  frequency  at  the  thermal  velocity  as 
veff  ~  The  fields  are  normalized  to  Q.c  =  wc(ur,  and  E  =  Ef  Ed,  where 

Ed  =  mvthVR/e  is  the  Dreicer  field,12  sufficient  to  cause  runaway  in  unmagnetized  systems. 
The  new  distribution  function,  normalized  to  1,  is  /(e)  =  (t>Jfc/2n)/(v).  It  obeys  the  normalized 
kinetic  equation,  which  is,  from  eq.  (1 1), 

~  =  CebW  +  CFP(t)  +  C*(e)  (39) 

where  the  terms  on  the  right-hand  side  are  given  below.  In  normalized  units,  the  characteristic 
time  in  which  any  term  will  significantly  modify  the  distribution  function  is  equal  to  the  magnitude 
of  that  term,  in  units  of  u^1.  The  Fokker-Planck  elastic  collision  term  is  derived  from  eq.  (14): 

c"  =  {*  jf  M***  +  If  [/AM***  +  ^  JTA(,)  *] } '  (40) 

The  small-angle  collisions  represented  by  Cfp  cause  the  electrons  to  approach  a  Maxwell- 
Boltzmann  distribution  from  any  initial  state.  Acting  against  the  tendency  to  approach  a 

Maxwellian  are  the  heating  term  represented  by  Ceb  and  the  radiative  cooling  term  represented 

by  C*(E).  The  importance  of  these  other  terms  in  driving  the  distribution  from  a  Maxwellian  can 
be  estimated  by  comparing  their  magnitude  to  Cfp- 
The  ohmic  heating  term  in  eq.  (39)  is,  from  eq.  (18), 

[9('>f  ]  <41> 


(37) 

(38) 
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where 


and  the  shape  function 


AEB  =  ^ 

Zt/R 


9(e)  = 


(0 


determines  where  the  heating  is  focused.  When  there  is  no  turbulence,  veff/vR  =  Z  (see  eq.  (36)). 
This  form  of  Ceb  is  most  convenient  when  wc  >  i/e//. 

The  inelastic-collision  term  in  eq.  (39)  is  just  as  in  eq.  (15),  but  with  the  collision  frequency 
normalized  by  dividing  by  vr: 

c' = E  hm*  -  •*) + «• + •  («> 

Here,  aj?s  «  1/137  is  the  fine  structure  constant,  /«  is  the  fraction  of  ions  in  state  a,  is  the 
multiplicity  of  state  a,  and  £lai(e)  is  the  collision  strength  discussed  in  connection  with  eq.  (9). 
£lab(e)  varies  slowly  with  e. 

The  effect  of  ohmic  heating  on  the  distribution  is  determined  by  both  the  relative  magnitude 
of  Ceb  and  die  shape  function  5(c),  which  contributes  to  the  magnitude  when  wc  <C  uK/f.  The 
relative  magnitude  is  roughly 

(45) 

Ohmic  heating  significantly  influences  the  distribution  function  only  where  \Ceb\  >  1.  If 
the  magnetic  field  is  small  (wc  <  *>«//),  \Pbb\  is  proportional  to  E7/veffVR,  and  electron- 
ion  collisions  and  turbulence  limit  the  influence  of  ohmic  heating  on  the  distribution.  The 
low-field  ohmic  heating  rate  is  inversely  proportional  to  (the  dimensional  heating  rate  is 
CsB(e)  multiplied  by  vr).  For  a  high  magnetic  field  (wc  >  ve//),  \Ceb\  is  proportional 
to  (* 'eff/vR)(E/B)2,  and  the  magnetic  field  limits  the  relative  influence  on  the  distribution. 
High-field  ohmic  heating  increases  linearly  with  1 /«//,  because  collisions  increase  the  current  flow 
by  interfering  with  the  magnetic  field  confinement. 

The  shape  function  g(t)  determines  the  form  of  the  ohmically  heated  distribution.  As  shown 
in  fig.  (1),  it  increases  rapidly  from  zero  and  asymptotically  approaches  its  maximum  value  of  1 .  In 
general,  low-energy  particles  (much  below  the  inflection  point)  experience  little  heating,  because 
the  large  low-energy  collision  rates  effectively  randomize  the  electron  motion  and  prevent  the 
electric  field  from  doing  net  work  on  electrons.  High-energy  particles  (much  above  the  inflection 
point)  are  easier  to  heat  up  because  more  phase  space  is  available  to  them. 

When  the  magnetic  field  is  low,  g(e)  rises  rapidly  (~  e8)  in  the  most  populated  energy  range. 
The  resulting  electron  distribution  is  greatly  enhanced  at  high  energies  compared  to  a  Maxwellian, 
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reflecting  the  runaway  phenomenon,  wherein  faster  electrons  are  more  efficiently  heated  than 
slower  electrons.  For  intermediate  magnetic  fields,  the  distribution  is  non-Maxwellian,  with  a 
local  enhancement  due  to  heating.  With  a  large  magnetic  field,  g(e)  is  nearly  constant  for  most 
electrons,  and  the  distribution  shape  is  self-similar  (as  discussed  below). 

In  a  moderate  magnetic  field,  ohmic  heating  results  in  a  distribution  function  with  an 
enhancement  at  intermediate  energies  but  with  a  depleted  high-energy  tail  relative  to  a  same- 
temperature  Maxwellian.  From  the  form  of  Ceb .  it  is  clear  that  this  enhancement  will  be 
somewhere  near  the  rapid  rise  of  g(e),  although  exactly  where  will  depend  on  the  shape  of  the 
distribution  function.  The  slope  of  g{e)  is  maximized  at  the  inflection  point 


eo  = 


(46) 


and  so  this  makes  a  convenient  estimate  of  the  focus  of  ohmic  heating.  Actually,  the  falloff  of 
the  distribution  function  with  energy  will  make  e0  an  underestimate.  Fig.  (2)  shows  numerically- 
determined  values  for  the  energy  of  greatest  ohmic  heating,  when  the  heated  distribution  is  a 
Maxwellian. 

For  very  weak  magnetic  fields,  the  ohmic  heating  focus,  with  its  distribution  function 
enhancement,  moves  out  to  higher  energies,  and  the  depletion  in  the  tail  above  eo  becomes 
less  important.  When  uc  -*  0,  the  enhanced-tail  runaway  distribution  is  produced.  For  very 
strong  magnetic  fields,  the  ohmic  heating  focus  moves  to  lower  energies,  and  when  wc  oo  the 
distribution  takes  on  a  depleted-tail,  or  “flat-topped”,  self-similar  form. 

Ohmic  heating  in  a  strong  magnetic  field  produces  a  “self-similar”  solution1*  for  the 
distribution,  in  the  same  form  as  the  distributions  produced  by  inverse-bremsstrahlung  laser 
heating.  When  uc  >  i^//,  then  g(e)  as  1,  and  if  inelastic  collisions  are  neglected  the  solution  to 
eq.  (39)  has  the  form: 

/-<*.*)  -  (47) 

Vm 

where  the  normalization  and  kinetic  energy  integrals  determine 


dm 

Vm 


m 


87rr(3/m) 
3  r(3/m) 


2  T(5/m) 


kT0 


(48) 

(49) 


and  kT0  =  1  for  the  normalized  distributions.  With  electron-electron  collisions  alone,  m  =  2 
(a  Maxwellian).  With  ohmic  heating  alone,  m  =  5.  With  ohmic  heating  plus  electron-electron 
collisions,  2  <  m  <  5:* 


m  =  2  + 


1  +  1.66/a0-724 


(50) 
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where 


The  ratio  of  the  self-similar  distribution  to  a  same-temperature  Maxwellian  is  plotted  in  fig.  (3). 

The  characteristics  of  ohmically-heated  distribution  functions  discussed  above  are  illustrated 
in  fig.  (4)  and  fig.  (5).  These  show  the  evolution  of  an  initially-Maxwellian  distribution  function 
under  the  influence  of  ohmic  heating  only.  In  fig.  (4),  the  magnetic  field  is  low:  u>c/i/e/ /  =  1  /99, 
and  the  enhancement  near  e  =  8  can  be  seen.  In  fig.  (5),  the  magnetic  field  is  moderately  high: 
Uc/veff  =  99,  and  the  distribution,  while  enhanced  near  e  =  1,  has  mostly  the  character  of  the 
self-similar  solution. 


The  inelastic  collision  term  has  a  different  (non-differential)  character  from  the  other  terms 
in  eq.  (39);  to  estimate  its  importance  relative  to  electron-electron  elastic  collisions,  we  first 
simply  take  the  ratio  of  the  inelastic  to  elastic  collision  frequencies,  given  by  eq.  (9)  and  eq.  (7), 


respectively 


(52) 


where 


A*  =  1.838  x  10“2 


ItT  faQib 
.ZlogAee  ga 


This  ratio  is  generally  small;  since  the  only  terms  not  of  order  one  are  Z  and  log  A,  which  are  each 
near  10,  an  effect  of  perhaps  a  few  percent  might  be  anticipated  from  inelastic  collisions.  Greater 
effects  should  appear  in  the  tail  of  the  distribution,  in  particular  at  energies  greater  than  Z  times  the 
thermal  energy.  The  effect  of  inelastic  collisions  will  be  significantly  greater,  however,  at  higher 
temperatures  (above  1  keV). 

The  effect  of  inelastic  collisions  is  to  enhance  the  distribution  function  tail,  relative  to  a 


same-temperature  Maxwellian.  This  is  paradoxical,  because  collisional  excitation  shifts  electrons 
out  of  distribution  tail,  into  the  bulk.  But  the  scattering  cools  the  system,  and  the  net  result  is  a 
relative  enhancement  in  the  tail,  compared  to  the  cooled  Maxwellian.  In  thermal  equilibrium,  or 
for  a  fixed  cooling  rate,  inelastics  tend  to  produce  enhancements  at  low  energy,  depletions  near 
the  threshold  energy,  and  enhancements  in  the  tail.  This  is  illustrated  in  fig.  (6),  which  plots  the 
time-development  of  an  initially  Maxwellian  distribution  under  the  influence  of  inelastic  collisions 
alone. 


Compared  to  ohmic  heating,  inelastic  collisions  are  much  more  important  at  low  electron 
energy,  but  much  less  important  at  high  energy.  For  a  Maxwellian  distribution,  the  approximate 
ratio  is 

Cob.rn.ic  y^gg  (c(y  +  e3)-3]e2 
CJnel  A*  (y  +  e8)2  ( 
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where 


y  =  (".///“c)2- 


(55) 


ihe  energy  dependence  of  this  ratio  is  shown  in  fig.  (7),  for  various  magnetic  field  strengths. 
The  figure  shows  the  ratio  of  ohmic  to  inelastic  terms,  divided  by  the  high-energy  limit  of  this 
expression,  which  is 


C ohmic  _  72.5  Z(i/eff/i/n)\og  A  /  £  \ 

@lnel  e—too  K'^'{fa^ab/ 5o)  \^c/ 


(56) 


5.  Conductivity  in  an  inductive  system 

The  kinetic  theory  presented  here  is  expressed  in  terms  of  the  electric  field  (E  or  £)  in  the 
plasma.  In  an  inductive  system  like  a  z*pinch,  however,  it  is  more  convenient  to  use  the  current 
density  j  as  a  parameter.  Current  and  electric  field  are  related  by  the  conductivity  cr,  which  in  a 
magnetized  system  is  a  tensor: 

J  —  a  ■  E  (57) 

The  resistivity  tensor  is  defined  as  17  =  <r-1.  We  will  now  show  how  to  calculate  “scalar” 
conductivities  (or  resistivities),  which  we  can  use  to  find  the  effective  electric  field  for  the  kinetic 
model.  This  can  also  be  used  to  correct  MHD  transport  codes  when  the  distribution  function  is 
non-Maxwellian. 

There  are  two  useful  ways  to  define  scalar  conductivities.  Both  ways  use  the  same  conductivity 
parallel  to  the  magnetic  field,  which  is  unaffected  by  the  field.  Consider  the  current  perpendicular 
to  the  magnetic  field.  In  the  first  approach,  the  “field  representation”,  the  electric  field  is  given 
and  the  scalar  conductivities  determine  the  components  of  the  current  density  vector  parallel  to 
and  perpendicular  to  it.  In  the  second  approach,  the  “current  representation”,  the  current  density 
is  given,  and  the  scalar  resistivity  determines  the  components  of  the  electric  field  parallel  and 
perpendicular  to  the  current. 

If  total  current  density  is  j  (transverse  to  the  magnetic  field)  and  total  electric  field  is  E,  the 

— * 

field  representation  Ohm’s  law  gives  the  current  components  j±  parallel  to  E  and  jA  perpendicular 
toE: 

j±  —  tr±E  (58) 

=  <taE  (59) 

The  components  of  the  current  distribution  function  fi(v)  in  these  directions  can  easily  be  found 
in  eq.  (17).  Similarly,  E  has  components  E±  parallel  to  and  EA  perpendicular  to  j,  which  are 
determined  the  current  representation  Ohm’s  law: 

E±  =  Tjxj  (60) 

Ea  =  t?Aj-  (61) 
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Note  that  it  is  not  true  that,  e.g.,  ij±  =  o-J1.  The  relationship  between  the  field- 
representation  conductivities  and  the  current-representation  resistivities  is  found  by  rotating 
from  one  representation  to  the  other: 


V±  = 


<r± 

+*A 


Va 


*I  +  *A 


(62) 

(63) 


The  current-representation  resistivities  are  those  most  often  used4’8  in  MHD  models,  and  rj±  is 
referred  to  as  the  “Spitzer”  resistivity. 

As  the  magnetic  field  increases  to  infinity,  the  field-representation  conductivities  go  to  zero: 
The  drift  velocity  becomes  vanishingly  small,  and  the  electron  gyroradius  also  vanishes,  eliminating 
the  effect  of  collisions.  This  is  called  magnetic  insulation.  Paradoxically,  the  current-representation 
resistivities  do  not  become  infinite  for  an  infinite  magnetic  field;  in  fact,  the  Spitzer  transverse 
resistivity  r)±  for  infinite  magnetic  field  is  well  known  to  be  only  about  twice  the  zero-field  value. 
To  explain  this  paradox,  we  take  as  an  example  a  simple  conductivity  model.14  By  following  the 
motion  of  a  single  electron  under  electric  field  acceleration  and  collisional  drag,  we  can  estimate 
the  field-representation  conductivity  as 


*1  = 


l  +  (u»c/l/)2’ 


_*  — 
°A  = 


(h>e/i/)cr| 

1  +  (wc/i/)2 


(64) 


where  v  is  a  collision  frequency  and  <rjj  =  ne2  /mu  is  the  zero-field  conductivity.  These 
conductivities  do  properly  go  to  zero  as  the  magnetic  field  becomes  infinite.  The  current- 
representation  resistivities  are: 


»li  =  (<rii)~1»  Va  =  1  (65) 

For  this  model,  the  “Spitzer”  resistivity  rj±  is  completely  unaffected  by  the  magnetic  field.  How 
can  this  be?  When  the  magnetic  field  is  increased,  the  cross-field  resistivity  rjA  also  increases, 
indicating  an  increase  in  the  electric  field  perpendicular  to  the  current.  If  the  magnetic  field 
increases  to  infinity,  so  does  this  cross-current  electric  field.  The  ratio  between  the  field  Ex  along 
the  current  and  the  field  EA  transverse  to  the  current  is 


E±/Ea  =  vju}e. 


(66) 


The  current  representation,  in  postulating  a  finite  current,  does  not  allow  a  magnetic  field  cutoff: 
the  resistivity  simply  determines  the  electric  fields  necessary  to  generate  the  given  current.  This 
resolves  the  paradox:  A  “Spitzer”  current  flowing  in  an  infinite  magnetic  field  is  driven  by  an 
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infinite  electric  field  transverse  to  the  current  flow.  In  the  field  representation,  by  contrast,  a  finite 
electric  field  is  postulated,  and  the  resulting  currents  are  calculated.  An  inductive  system  (like  a 
PRS)  is  best  described  in  the  current  representation,  since  the  driver  is  capable  of  producing  large 
electric  fields  to  maintain  the  current  flow. 

In  our  kinetic  model,  the  field  representation  conductivities  can  be  straighforwardly  calculated 
from  eq.  (17),  using  eq.  (4): 


r  ~  du  i. 

3  v*lf  Jo  1  +  (“e/"*//)2*3  Se 
=  lxL(xP\  P  £,/i  af‘ 

3u.jt  U/,7  Jo  l+WK.ffJV  9' 


(67) 

(68) 


This  form  is  most  convenient  for  iac  <  t /e//,  but  these  expressions  are  general.  The  Spitzer 
resistivity  can  now  be  found  from  eq.  (62).  In  the  highly-magnetized,  large-u)c  limit,  where  the 
self-similar  solution  of  eq.  (47)  is  valid,  the  Spitzer  resistivity  has  a  simple  limit: 


lim  Tj±  = 
UJe—+  oo 


32tt2  u% 

~  5  /o(0)‘ 


(69) 


In  the  unmagnetized  limit,  the  expression  involves  an  integral  which  must  be  evaluated  for  the 
particular  distribution  function.  In  any  case,  the  resistivity  is  of  the  form 


Vx  = 


Z  log  A 

(«T)3/2 


sec. 


(70) 


The  constant  Ax  is  shown  in  Table  2  for  integral  values  of  the  self-similar  parameter  m,  both  for 
unmagnetized  and  strongly  magnetized  plasma,  with  no  turbulence  included. 


Table  2.  Resistivity  coefficient  A±  for  different  distributions 


m  —  2 

771  =  3 

m  =  4 

m  =  5 

o 

II 

cq 

3.3798  x  10"15 

3.5612  x  10~15 

3.6521  x  lO"15 

3.7056  x  10"15 

B  —  oo 

1.1475  x  10~14 

7.1223  x  10~15 

5.7471  x  10_1S 

5.1133  x  10"15 

In  the  unmagnetized  plasma,  a  Maxwellian  distribution  has  the  highest  resistivity,  because 
higher-m  distributions  have  fewer  fast,  tail  electrons.  The  resistivity  for  an  m  =  5  distribution 
is  slightly  (0.91  times)  lower  than  a  Maxwellian.  In  the  strongly-magnetized  plasma,  conversely, 
the  Maxwellian  has  the  lowest  resistivity.  The  strong  magnetic  field  efficiently  confines  the  tail 
electrons  and  makes  the  low-energy,  bulk  electron  the  current  carriers.  The  magnetized  resistivity 
of  the  77i  =  5  distribution  is  less  than  that  of  a  Maxwellian  by  a  factor  of  2.24. 

Table  3  lists  the  resistivities  of  several  elements  (most  of  which  are  used  as  PRS  radiators),  for 
comparison  with  the  magnitudes  of  the  plasma  resistivities  presented  here. 
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Table  3.  Selected  electrical  resistivities15 


Given  in  units  of  10  18  sec;  1  microohm-cm  =  |  x  10  17  sec. 

ELEMENT 

77°  K 

273°  K 

373°  K 

Mg 

0.69 

4.33 

6.22 

Fe 

0.73 

9.89 

16.3 

Cu 

0.22 

1.73 

2.49 

Ag 

0.33 

1.68 

2.37 

Au 

0.56 

2.27 

3.16 

Zn 

1.22 

6.11 

8.67 

A1 

0.33 

2.72 

3.94 

Se 

15.0 

13.3 

- 

In  the  presence  of  turbulence,  the  resistivity  uniformly  increases  due  to  the  increase  in  i/c//, 
although  the  extent  of  the  increase  depends  most  importantly  on  the  magnetic  field  strength,  form 
of  the  distribution  function,  and  temperature.  In  fig.  (8),  the  increase  in  resistivity  with  turbulence 
strength  is  shown  for  the  m  —  2  and  m  —  5  distributions  in  a  large  magnetic  field.  With  a  weak 
magnetic  field,  the  curves  would  look  similar,  but  the  Maxwellian  (m  =  2)  resistivity  would  be 
less  than  that  of  the  m  =  5  distribution,  and  the  effect  of  turbulence  would  be  less  pronounced. 

By  symmetry,  the  current  density  j  must  be  along  the  pinch  ( z )  axis.  This  situation  is 
enforced  by  the  presence  of  radial  charge-separation  electric  fields  (E±,EA).  As  we  have  shown, 
the  radial  electric  fields  can  become  extremely  large.  The  total  electric  field  makes  an  angle 
Oj  =  tan~1(<r|j/<rj.)  with  the  pinch  axis.  The  magnitude  of  the,  total  electric  field,  used  in  our 
kinetic  model,  is  related  to  the  total  (axial)  current  density  by: 


For  a  given  current  density,  E  increases  with  magnetic  field,  roughly  as  B  for  high  fields. 

The  electron  flow,  as  given  by  /i(v),  is  not  generally  parallel  to  the  current,  but  rather  changes 
direction  with  velocity.  Thus,  there  are  always  electron  flows  intermixing  throughout  the  radial 
extent  of  the  pinch.  These  flows  have  the  potential  to  both  enhance  radial  transport  and  seed 
instabilities  not  discussed  here.  Such  anisotropies  in  the  flow  are  discussed  below. 

6.  Runaway  electrons 

Because  the  electron  collision  rates  decrease  rapidly  with  relative  energy  (see  eq.  (7)),  particles 
with  large  enough  kinetic  energy  are  not  affected  by  collisions  and  become  runaways.  With  no 
magnetic  field,  “runaway”  electrons  are  those  with  a  velocity  greater  than  the  critical  velocity12 
ve  «  vthy/2ED/E.  The  critical  velocity  vc  is  independent  of  vth,  but  as  the  unmagnetized  system 
is  heated  and  vth  increases,  an  increasing  number  of  electrons  have  velocities  above  vc  and  thus 


A  magnetic  field  inhibits  this  free  acceleraton  of  electrons.  In  a  collisionless  plasma  with 
B  >  E  (in  Gaussian  units),  free-streaming  electron  motion  across  field  lines  is  impossible,  being 
replaced  by  a  combination  of  cyclotron  and  E  x  B  drift  motions.  When  B  <  E,  collisionless 
electrons  can  be  accelerated  across  field  lines,  although  the  magnetic  field  reduces  the  magnitude 
of  the  acceleration.  This  high-electric  field  acceleration  might,  in  connection  with  neck  formation, 
cause  z-pinch  runaways.16  The  large  radial  fields  required  for  axial  current  propagation  in  the 
magnetized  system  could  also  promote  this  type  of  runaway,  as  could  the  presence  of  turbulence. 
Predictions  of  electric  and  magnetic  fields  inside  the  pinch  require  a  model  of  how  current  and 
fields  penetrate  inside  and  their  consequent  relationship. 

In  the  bulk  of  the  current-carrying  region,  B  >  E,  and  high-energy  electrons  are  confined  by 
the  magnetic  field.  To  identify  “runaways”,  we  look  for  combinations  of  parameters  that  generate 
directional  and/or  monoenergetic  beams.  Finite-size  effects  can  be  estimated  using  the  mean  free 
path  and  the  physical  dimensions  of  the  system.  Or,  the  rate  at  which  high-energy  electrons  are 
produced  can  be  found  by  defining  a  runaway  velocity  threshold  and  calculating  the  electron  flux 
past  this  threshold. 

Anisotropy  in  the  electron  velocity  distribution  is  given  by  fi(v):  Any  fi{v)  >  0  indicates 
a  directed  beam  at  velocity  v.  The  angle  between  fi(v)  and  the  electric  field  is  easily  computed. 
Both  lie  in  the  r-z  plane,  perpendicular  to  B\  it  follows  that  the  angle  6  between  fi{y)  and  the  z 
axis  is 

tan#  =  8,  -  — c3/2,  (72) 

Ueff 

where  Oj  is  the  axial  angle  of  E.  For  a  large  magnetic  field,  this  anisotropy  between  average 
current  and  individual  particle  flow  acts  to  isotropize  the  distribution;  electron-ion  collisions,  by 
interfering  with  the  magnetic  field,  act  to  preserve  the  electric-field  anisotropy.  Just  the  opposite 
happens  with  a  low  magnetic  field,  where  electron-ion  collisions  isotropize  the  distribution. 

Future  work  on  runaways  will  be  directed  towards  determining  how  current  and  fields  penetrate 
into  the  pinch.  Given  the  current  and  field  distribution,  the  kinetic  theory  developed  here  can  be 
utilized  to  investigate  where  and  how  high-energy  and  high-anisotropy  electron  flows  originate. 

7.  Model  calculations 

The  actual  electron  distribution  function  is  a  product  of  the  interplay  between  electron -electron 
collisions,  ohmic  heating  and  electron-ion  excitation  collisions,  as  described  above.  To  investigate 
the  behavior  of  the  distribution  function  in  high-power  PRS  machines  such  as  DECADE,  we  have 
numerically  solved  eq.  (39)  for  the  time-evolving  distribution.  In  particular,  we  are  looking  for 
an  equilibrium  distribution  shape,  which  will  serve  as  a  generalization  of  the  self-similar  form 
of  eq.  (47)  to  moderate  magnetic  fields  and  the  presence  of  inelastic  collisions.  Therefore,  we 
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set  the  ohmic  heating  strength  Aeb  such  that  the  ohmic  heating  rate  equals  the  inelastic  cooling 
rate,  which  is  dependent  only  on  the  plasma  parameters  (e.g.,  temperature,  ionization,  electron 
distribution  function).  We  neglected  turbulence;  any  turbulence  would  act  like  an  increased 
ionization  level  Z. 


Table  4.  Parameters  used  in  numerical  solutions 


Case  B  (MG) 

epocu.  Current  (A/cm2) 

Cj  mei 

_ 

kT  =  500  eV,  Z  = 

24.01,  Radiative  power  3  x  101S 

W 

1  0.93 

9.0  6  x  10® 

3.08 

2  29.26 

4.6  8  x  10® 

0.19 

kT  =  1500  eV,  Z  = 

=  28.96,  Radiative  power  8  x  1015  W 

3  0.335 

9.0  7  x  10® 

0.50 

4  10.54 

4.6  9  x  10® 

0.35 

The  numerical  solution  was  for  a  selenium  plasma,  of  ion  density  n,-  =  1020  cn«~s.  The 
parameters  are  summarized  in  table  4.  Hie  atomic  model  used  for  selenium  included  95  transitions 
for  different  ionization  stages  between  the  ground  state  and  first  excited  state;  of  these,  35  were 
found  to  make  some  contribution  under  the  given  plasma  conditions.  Ion  populations  were 
determined  by  a  collisional-radiative  equilibrium  code.  The  transitions  used  are  numbered  and 
listed  in  table  5.  The  relative  importance  of  each  transition  can  be  seen  in  the  bar  charts  of  fig.  (9) 
and  fig.  (10),  which  show  the  contribution  of  each  transition  to  the  total  radiative  power  for  each 
of  the  two  temperatures  considered. 

The  model  results  are  shown  in  fig.  (1 1),  fig.  (12),  fig.  (13)  and  fig.  (14),  respectively,  for  the 
four  cases  described  in  Table  4.  For  a  given  temperature,  increasing  the  magnetic  field  changes 
the  distribution  from  the  runaway  distribution  of  fig.  (11)  to  the  self-similar  type  distril  Jtion 
of  fig.  (12).  The  low-energy  enhancement  of  the  ohmically-heated  distribution  seen  in  fig.  (5) 
is  missing  from  fig.  (12)  because  of  the  rapid  isotropization  from  electron-electron  collisions  at 
low  energies.  For  a  given  magnetic  field,  increasing  the  temperature  significantly  reduces  the 
non-Maxwellian  effects  for  these  parameters.  This  is  both  because  of  the  increasing  influence 
at  higher  temperatures  of  inelastic  collisions  relative  to  ohmic  heating  and  because  the  higher 
conductivities  at  high  temperatures  lead  to  lower  electric  fields  in  the  plasma.  Still,  distribution 
function  enhancements  at  higher  temperatures  do  generate  significant  numbers  of  electrons  at  high 
energies  (e.g.,  at  1500  eV,  e  =  10  corresponds  to  15  keV). 
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Table  5.  Transitions  used  in  numerical  solutions 


Index 

Threshold  (eV) 

Ionization 

Transition 

1. 

62.1 

Li-like 

0-1 

2. 

1956.7 

Li-like 

1  -  2 

3. 

68.1 

Be-like 

0-1 

4. 

1954.5 

Be-like 

0-2 

5. 

1886.5 

Be-like 

1  -  2 

6. 

87.7 

B-like 

0-1 

7. 

1712.3 

B-like 

0-2 

8. 

1624.5 

B-like 

1  -2 

9. 

117.3 

C-like 

0-1 

10. 

1765.4 

C-like 

0-2 

11. 

1648.2 

C-Iike 

1  -  2 

12. 

156.4 

N-like 

0-1 

13. 

1688.3 

N-like 

0-2 

14. 

1531.8 

N-like 

1  -2 

15. 

163.6 

O-like 

0-1 

16. 

1597.1 

O-like 

0-2 

17. 

1433.5 

O-like 

1-2 

18. 

1921.1 

O-like 

1-3 

19. 

180.0 

F-like 

0-1 

20. 

1642.2 

F-like 

0-2 

21. 

1462.2 

F-like 

1-2 

22. 

1890.2 

F-like 

1-3 

23. 

1897.4 

F-like 

1-4 

24. 

1530.3 

Ne-like 

0-1 

25. 

1996.4 

Ne-like 

0-2 

26. 

466.1 

Ne-like 

1-2 

27. 

663.8 

Ne-like 

1-3 

28. 

765.0 

Ne-like 

1-4 

31. 

468.7 

Na-like 

1-2 

33. 

656.3 

Na-like 

1-3 

35. 

750.5 

Na-like 

1-4 

8.  Conclusions 

We  have  presented  a  model  that  consistently  includes  the  effects  of  electron-electron  collisions, 
ohmic  or  electric-field  heating,  distribution  function  cooling  by  inelastic  collisions  with  ions,  and 
ion-acoustic  turbulence.  The  actual  distribution  function  shape  depends  on  the  interaction  between 
these  factors. 

At  moderate  temperatures  (under  1  keV),  the  ohmic  heating  term  is  the  dominant  influence  on 
the  distribution.  This  is  true  even  when  the  ohmic  heating  and  inelastic  cooling  rates  are  equal,  as 
shown  in  fig.  (1 1)  and  fig.  (12).  The  magnitude  of  the  ohmic  term  is  proportional  to  the  square  of 
the  electric  field  (cf.  eq.  (41)),  which  is  in  turn  proportional  to  the  current  density  and  the  magnetic 
field.  The  relationship  is  given  in  eq.  (71);  a  rough  estimate  comes  from  using  the  simple  model 
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of  eq.  (64): 


(73) 


where  the  dimensionless  quantity  £  is  the  ratio  of  the  electric  field  to  the  Dreicer  field.  The  exact 
electric  field  can  always  be  calculated  from  eqs.  (67-68).  Since  £  is  determined  by  the  current,  we 
can  define  a  minimum  current  density  for  which  ohmic  heating  is  important  (see  the  discussion  of 
eq.  (45)).  The  condition  will  depend  on  the  magnetic  field,  density  and  temperature  at  each  point, 
and  so  assessing  the  importance  of  ohmic  heating  in  any  PRS  will  require  a  magnetohydrodynamic 
model  that  predicts  both  current  penetration  and  heat  flow.  As  we  have  shown,  ohmic  heating 
effects  on  the  distribution  are  directly  multiplied  by  microturbulence. 

Inelastic  collisions  become  more  important  at  higher  temperatures  (cf.  eq.  (53))  and  could 
modify  the  distribution  by  up  to  10%  at  keV  temperatures.  Inelastic  collisions  oppose  high- 
magnetic  field  ohmic  heating,  cooling  the  distribution  so  as  to  produce  a  relatively  enhanced  tail. 
At  high  temperatures  and  currents,  ohmic  heating  and  inelastics  could  balance  each  other  and 
produce  a  distribution  that  was  close  to  Maxwellian.  Turbulence  would  disrupt  this  balance. 

When  strong-magnetic  field  ohmic  heating  dominates,  the  self-similar  distribution  /«(  c)  is  a 
good  analytic  approximation  to  the  actual  distribution  function.  In  hotter  plasmas,  or  for  lower 
current  density,  or  when  the  current  is  concentrated  on  the  plasma  edge  so  that  25/2?  >  1,  another 
analytical  approximation  must  be  found  for  the  distribution  function  form.  It  is  always  true  that 
the  distributrion  function  can  be  found  numerically,  by  integrating  eq.  (39). 

Establishing  the  scaling  of  non-Maxwellian  effects  for  DECADE  shot  parameters  is  the  next 
important  step.  For  this,  we  will  need  an  analytic  model  of  the  current  and  field  distribution  in 
the  device.  We  will  first  assume  a  uniform  axial  current  density,  investigating  non-Maxwellian 
distributions  produced  in  an  energy  equilibrium.  This  will  establish  a  baseline  for  non-Maxwellian 
production,  and  allow  us  to  estimate  the  change  in  PRS  radiative  output  from  that  predicted  using 
a  Maxwellian  distribution.  We  can  then  more  realistically  use  a  Bennett  pressure  equilibrium  to 
relate  current  and  temperature  profiles,  and  determine  the  effect  of  non-Maxwellian  distributions 
and  turbulence17  on  the  Pease-Braginskii  current  and  radiative  collapse.  Radiative  effects  in 
DECADE-class  z-pinches  are  likely  to  be  strongly  influenced  by  these  phenomena. 
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Figure  Captions 

Figure  1.  Plot  of  the  shape  function  g(t)  for  different  values  of  u >c/veff-  From  left  to  right,  the 
curves  are  for  =  10, 1  and  0.1,  respectively. 

Figure  2.  Numerically  determined  ohmic  heating  focus  energy  as  a  function  of  u)c/vef /.  The  focus 
energy  is  the  point  of  greatest  enhancement  of  an  initially  Maxwellian  distribution  function  when 
affected  by  ohmic  heating  only.  Each  curve  is  for  a  different  distribution  function  temperature 
kTo'.  The  middle  curve  is  for  the  initially  normalized  distribution,  kTq  =  1,  the  upper  curve  is  for 
kTo  =  20  and  the  bottom  curve  is  for  kT0  =  0.1. 

Figure  3.  Ratio  of  the  self-similar  distribution  function  /m(e)  with  a  same-temperature  Maxwellian. 
Curves  A,B,C,D  correspond  to  integral  m  values  of  2, 3, 4  and  5,  in  order  of  greater  tail  depletion. 

Figure  4.  Evolution  of  an  initially-Maxwellian  distribution  under  ohmic  heating  only,  in  a  low 
magnetic  field.  Here  u?c/i/e//  =  1/99  and  kT  =  500  eV.  Shown  is  the  ratio  of  the  distribution 
'ith  a  same-temperature  Maxwellian;  the  curves  show  the  evolution  of  distribution  shape  as  it 
heats,  from  numerical  solution  of  eq.  (39). 

Figure  5.  Evolution  of  an  initially-Maxwellian  distribution  under  ohmic  heating  only,  in  a  high 
magnetic  field.  Here  wc/i/e//  =  99,  and  kT  =  500  eV.  Shown  is  the  ratio  of  the  distribution  with 
a  same-temperature  Maxwellian;  the  curves  show  the  evolution  of  distribution  shape  as  it  heats, 
from  numerical  solution  of  eq.  (39). 

Figure  6.  Evolution  of  an  initially-Maxwellian  distribution  under  inelastic  cooling  only.  Shown  is 
the  ratio  of  the  distribution  with  a  same-temperature  Maxwellian;  the  curves  show  the  evolution  of 
distribution  shape  as  it  cools,  from  numerical  solution  of  eq.  (39). 

Figure  7.  Ratio  of  magnitudes  of  the  ohmic  heating  and  inelastic  cooling  terms,  as  a  function  of 
energy.  This  plot  gives  only  the  dependence  of  e;  the  absissa  is  in  units  of  Abb /A*,  as  given  in 
eq.  (56). 

Figure  8.  Electrical  resistivity  (sec)  as  a  function  of  turbulence  level  Veff/Vti-  Solid  curve  is  for  a 
Maxwellian,  dashed  curve  is  for  the  m  —  5  depleted-tail  distribution.  A  strong  (1  MG)  magnetic 
field  emphasizes  the  difference  between  curves. 

Figure  9.  Power  density,  in  W/cms,  generated  by  each  selenium  inelastic  collisional  transition 
(listed  in  Table  5),  at  kT  —  500  eV.  The  total  radiative  power  density  was  3  x  1015  W/cm~5  at 
this  temperature. 
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Figure  10.  Power  density,  in  W/cm3,  generated  by  each  selenium  inelastic  collisional  transition 
(listed  in  Table  5),  at  kT  =  1500  eV.  The  total  radiative  power  density  was  8  x  101S  W/cm_s  at 
this  temperature. 

Figure  1 1.  Numerical  results  for  the  evolution  of  the  distribution  function  in  time  under  energy- 
equilibrium  conditions.  Shown  is  the  ratio  of  the  distribution  function  with  a  same-temperature 
Maxwellian.  Time  between  successive  curves  A,B,C. . .  is  ree-  This  is  low  magnetic  field,  low 
temperature  case  1:2?  =  0.93  MG,  kT  =  500  eV.  Other  parameters  are  shown  in  Table  4. 

Figure  12.  Numerical  results  for  the  evolution  of  the  distribution  function  in  time  under  energy- 
equilibrium  conditions.  Shown  is  the  ratio  of  the  distribution  function  with  a  same-temperature 
Maxwellian.  Time  between  successive  curves  A,B,C. . .  is  ree.  This  is  high  magnetic  field,  low 
temperature  case  2:  B  =  29.26  MG,  kT  =  500  eV.  Other  parameters  are  shown  in  Table  4. 

Figure  13.  Numerical  results  for  the  evolution  of  the  distribution  function  in  time  under  energy- 
equilibrium  conditions.  Shown  is  the  ratio  of  the  distribution  function  with  a  same-temperature 
Maxwellian.  Time  between  successive  curves  A,B,C. . .  is  ree.  This  is  low  magnetic  field,  high 
temperature  case  3:  B  =  0.34  MG,  kT  =  1500  eV.  Other  parameters  are  shown  in  Table  4. 

Figure  14.  Numerical  results  for  the  evolution  of  the  distribution  function  in  time  under  energy- 
equilibrium  conditions.  Shown  is  the  ratio  of  the  distribution  function  with  a  same-temperature 
Maxwellian.  Time  between  successive  curves  A,B,C. . .  is  ree.  This  is  high  magnetic  field,  high 
temperature  case  4:  B  =  10.54  MG,  kT  =  1500  eV.  Other  parameters  are  shown  in  Table  4. 
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IV.  PRS  NOZZLE  DESIGN  AND  GAS  PUFF  SIMULATION 

(a)  PI  Simulations  and  Nozzle  Designs 

During  this  past  year,  considerable  effort  was  concentrated  on  developing  a 
better  understanding  of  the  2D  dynamics  of  gas  puff  implosions  associated  with  dif¬ 
ferent  nozzle  designs.  In  particular,  the  role  of  tilting  nozzles  inward  was  examined 
in  order  to  minimize  the  so  called  “zippering”  effect  in  which  the  implosion  first 
reaches  the  axis  near  the  nozzle  throat  and  then  progresses  up  the  axis.  This  type 
of  2-D  implosion  results  in  a  lower  peak  power  and  a  longer  radiation  pulse  width 
than  would  be  expected  from  a  uniform  1-D  implosion.  However,  if  the  nozzle  is 
tilted  inwards,  the  peak  power  of  K-shell  emission  is  higher  and  the  pulse  width  is 
shortened.  The  density  on  the  central  axis  at  the  time  of  assembly  also  appears 
more  uniform.  It  would  appear  that  the  total  yield  would  then  be  approximately 
independent  of  tilt  angle  and  that  the  tilted  nozzle  design  would  be  desirable  purely 
from  the  standpoint  of  higher  peak  power  during  the  pulse.  This  is  not  the  case, 
however.  In  our  simulations,  and  in  experimental  results  from  Physics  International, 
the  tilted  nozzles  also  produced  higher  yields  in  addition  to  the  higher  powers.  This 
is  because  the  final  assembly  on  axis  occurs  over  the  full  length  of  the  pinch  in 
the  tilted  nozzles  with  the  result  that  the  average  density  is  substantially  higher. 
Evidently  the  quality  of  the  tilted  nozzle  implosions,  in  terms  of  density  uniformity, 
is  such  that  the  resultant  large  increase  in  peak  power  also  leads  to  greater  yields 
despite  the  fact  that  the  FWHM  time  of  the  pulse  is  shorter. 


Nozzle  Width 

Yield  (kJ) 

Yield  (kJ) 

Yield  (kJ) 

(cm) 

(0  deg.  tilt) 

(-7.5  deg.  tilt) 

(-10  deg.  tilt) 

0.8 

1.4 

3.0 

4.6 

0.4 

2.0 

3.7 

7.1 

Table  I.  Double  Eagle  Argon  K-shell  yields  from 
different  nozzle  designs  (50  fx gm/cm) 
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Table  I  shows  the  Argon  K-shell  yield  values  from  six  Double  Eagle  nozzle  de¬ 
sign  simulations.  It  can  be  seen  that  there  is  about  a  4:1  ratio  in  the  argon  K-shell 
yield  between  the  best  case  (tilted  nozzle,  0.4  cm  width)  to  the  worst  case  (untilted 
nozzle,  0.8  cm  width).  This  is  roughly  the  ratio  seen  in  the  PI  experiments.  The 
simulation  results  have  shown  that  for  a  straight  nozzle,  the  bulk  of  the  radiation 
comes  from  the  dense  region  which  moves  up  the  z-axis  during  zippering.  In  con¬ 
trast,  when  the  nozzle  is  tilted  the  radiative  power  is  spread  more  uniformly  over 
the  plasma  volume.  In  addition,  the  average  bulk  density  over  the  entire  plasma  fill 
is  twice  as  high  in  the  best  yield  case  versus  the  worst  case  (1  x  1019  cm-3  versus  5  x 
1018  cm-3  shown  in  Table  I).  An  explanation  for  the  lower  densities  in  the  untilted 
nozzle  designs  is  as  follows.  During  a  “zippering”  implosion,  the  large  amount  of 
axial  kinetic  energy  which  is  generated  results  in  preheated  mass  being  pushed  up 
the  axis.  As  a  result,  as  the  implosion  progresses  the  remaining  mass  implodes  onto 
this  hot  plasma  which  limits  its  compression  and  final  density. 

Figures  1  and  2  show  density  fill  plots  for  the  worst  case  (0.8  cm  nozzle  width, 
0°  tilt)  and  best  case  (0.4  cm  nozzle  width,  —10°  tilt)  nozzles.  Figure  1  (a)  and  Fig. 
2  (a)  show  the  initial  density  profile  while  Fig.  1  (b)  and  Fig.  2  (b)  show  the  density 
near  final  assembly  at  the  time  of  peak  power.  The  importance  of  these  results  is 
that  in  a  straight  nozzle  the  radiative  output  is  localized  and  persists  during  the 
so-called  “zippering”  time.  This  leads  to  the  relatively  broader  radiation  pulses 
observed  in  implosions  using  these  nozzles.  In  a  more  uniform  implosion,  with  a 
mass  loading  initiated  by  flow  from  a  tilted  nozzle,  both  higher  peak  power  and 
higher  yields  are  observed. 

Historically,  one  of  the  troubling  problems  in  Z-pinch  simulations  is  the  fact 
that  the  radiation  pulse  width  seen  in  the  simulations  tends  to  be  of  much  higher 
power  and  shorter  duration  than  that  observed  in  experiments.  Despite  this  prob¬ 
lem,  however,  there  is  still  very  good  agreement  with  K-  and  L-  shell  yields.  Several 
explanations  have  been  given  for  this.  These  include  2-D  effects  in  which  instabili¬ 
ties  and  non-uniform  implosions  lead  to  different  times  of  collapse  for  different  parts 
of  the  plasma.  Another  theory  has  been  that  enhanced  transport  coefficients  in  the 
experiment  tend  to  broaden  the  current  and  thermal  profiles  and  hence,  soften  the 
implosion. 
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In  order  to  address  this  issue,  a  series  of  simulations  was  conducted  which 
compared  2-D  results,  obtained  with  tilted  and  untilted  nozzles,  against  1-D  results. 
For  a  series  of  simulations  of  the  PI  Double  Eagle  device  over  a  range  of  linear  mas6 
densities  and  enhanced  resistivities  we  found  that  that  the  yields  are  nearly  the  same 
in  2-D  calculations  as  in  1-D  calculations  but  that  the  width  of  the  radiative  power 
pulse  and  the  value  of  the  peak  power  can  be  very  different.  In  fact,  for  a  series  of 
Double  Eagle  implosions,  we  found  that  for  three  mass  densities  (50  /xgm/cm,  100 
/igm/cm,  and  150  /xgm/cm)  the  ratio  of  the  2-D  radiative  power  curve  pulse  width 
to  the  1-D  pulse  width  was  almost  the  same  with  a  value  close  to  2.5.  A  standard 
straight  nc  zle  (i.e.  typical  of  the  design  used  in  the  experiment)  with  a  10  degree 
gas  expansion,  and  a  3  cm  length  was  used.  We  then  went  back  and  repeated  the 
1-D  calculations  but  increased  the  resistivity  with  a  constant  multiplier,  17*.  We 
chose  rj*  factors  of  1,  10,  and  100.  Increasing  this  value  resulted  in  a  broadened 
pulse  width  as  well  as  some  increase  in  the  yield.  When  17*  =  100  and  with  the 
highest  mass  loading,  the  1-D  pulse  width  came  to  within  85%  of  the  2-D  value.  It 
must  be  emphasized,  however,  that  the  mechanisms  for  broadening  the  pulse  widths 
(at  least  for  the  case  of  a  straight  nozzle)  are  different  in  the  two  cases.  The  2-D 
broadening  is  due  to  “zippering”  or,  motion  along  the  axis,  while  the  broadening 
in  the  1-D  simulations  is  due  to  a  penetrating  magnetic  field  which  softens  the 
implosion.  Figure  3  summarizes  these  results  in  graphical  form. 

The  PI  simulation  results  show  the  same  qualitative  trend  as  the  experimental 
results.  However,  there  wa3  a  systematic  error  in  the  yields  obtained  from  the 
simulations  as  compared  with  the  experiment  in  that  the  simulation  yields  were 
consistently  less  than  were  observed  experimentally.  There  may  be  several  reasons 
for  this.  First,  these  simulations  are  designed  primarily  to  test  hypotheses  about 
implosion  uniformity.  Since  the  initial  mass  distribution  formed  from  a  nozzle  can 
fill  nearly  the  entire  computational  domain,  the  orthogonal  mesh  is  fixed  in  time. 
Hence,  resolution  for  the  core  plasma  is  reduced  over  what  could  be  achieved  with  a 
1-D  Lagrangian  or  Eulerian  with  a  moving  mesh  calculation.  In  these  simulations, 
there  were  approximately  100  radial  zones.  Therefore,  at  the  final  point  of  assembly, 
the  resolution  in  the  K-shell  emitting  zones  is  reduced  to  just  a  few  zones.  This 
necessarily  reduces  the  peak  density  which  can  be  achieved  and  could  lead  to  a 
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systematic  error.  A  second  possibility  is  that  the  pulse  width  for  the  experiment  is 
still  much  larger  than  those  which  are  seen  in  the  simulations.  The  experimental 
widths  are  approximately  double  that  of  the  simulations  although  the  trend  of 
shorter  pulses  in  tilted  vs.  untilted  nozzles  holds.  It  appears  that  there  are  still  other 
physical  processes,  besides  2-D  effects,  which  will  have  to  be  invoked  to  explain  the 
broader  pulses.  This  is  an  ongoing  and  important  topic  for  future  research.  Despite 
these  problems,  in  can  be  said  that  the  modeling  has  resulted  in  good  agreement 
with  the  PI  experiments  since:  (a)  the  ratio  of  the  yield  of  untilted  and  tilted  nozzles 
was  nearly  the  same  (4:1  from  best  case  to  worst  case)  and;  (b)  the  trend  in  the 
reduction  of  pulse  widths  is  the  same  in  the  simulations  as  in  the  experiment.  The 
model  also  confirms  the  experimental  hypothesis  that  jetting  occurs  along  the  axis 
during  a  “zippered  implosion”.  These  results  have  been  compiled  into  a  paper  which 
has  been  accepted  for  publication  in  Physics  of  Fluids  (to  appear  in  the  March,  1993 
issue). 

(b)  Saturn  Simulations 

With  the  success  of  the  interaction  with  PI  and  the  qualitative  agreement 
between  simulation  and  experiment,  it  was  felt  that  the  2-D  code  was  sufficiently 
benchmarked  and  could  be  used  to  look  at  implosion  problems  on  other  pulse  power 
machines.  The  goal  has  been  to  advance  the  simulation  capability  through  bench¬ 
marking  with  existing  machines  while  at  the  same  time  contributing  to  the  under¬ 
standing  and  improvement  of  ongoing  experiments. 

With  this  in  mind,  our  next  step  was  to  examine  gas  puff  implosion  dynamics  on 
the  Sandia  National  Laboratory  Saturn  pulse  power  device.  This  was  a  particularly 
advantageous  choice  because  the  nozzle  design  used  in  Saturn  experiments  was 
similar  to  that  used  in  the  Double  Eagle  experiments. 

In  these  simulations,  we  used  a  circuit  model  to  drive  the  implosion.  This 
model  consisted  of  the  characteristic  voltage  waveform  of  the  SATURN  generator 
coupled  to  a  circuit  which  consisted  of  a  machine  inductance  of  9.75  nh,  0.17  ohms 
of  machine  resistance,  and  return  current  bars  set  so  that  the  initial  load  inductance 
is  2.5  nh.  The  simulated  gas  puffs  were  argon. 

A  modified  ”thin”  radiation  lookup  table  routine  was  used  for  both  equation 


of  state  properties  and  radiation  power.  A  ”thin”  radiation  model  implies  that  ail 
radiation  energy  is  immediately  lost  from  the  emitting  computational  zone  of  the 
plasma  and  is  not  absorbed  in  other  parts  of  the  pinch  plasma.  In  the  initial  phase 
of  performing  these  simulations,  for  lack  of  opacity,  it  was  found  that  an  overly  large 
amount  of  radiated  power  was  coming  out  in  the  L-shell.  In  an  attempt  to  correct 
this  problem,  a  scheme  was  developed  for  modeling  the  effective  opacity  within 
each  computational  zone.  This  allows  us  to  continue  to  use  the  lookup  tables  but 
requires  an  additional  parameter  when  making  the  table,  namely  the  effective  length 
over  which  to  calculate  the  optical  depth.  Once  the  length  is  chosen,  the  effect  of 
opacity  within  each  individual  zone  and  for  each  table  entry  is  calculated  with 
lowered  the  Einstein  A’s  obtained  by  dividing  by  the  optical  depth  (i.e.,  A  =  A/ r). 
Several  tests  showed  that  this  can  dramatically  effect  the  L-shell  line  radiation  by 
decreasing  it  by  factors  of  two  or  more.  However,  the  K-shell  is  relatively  unaffected. 
We  now  feel  that  this  “local”  approximation  allows  us  to  more  accurately  model  2- 
D  Z-pinch  implosions.  In  addition  to  radiative  power,  the  equation-of-state  lookup 
tables  include  internal  energy  and  density  based  values  for  electron  temperature 
and  charge  state.  These  values  are  then  used  to  determine  the  electron  pressure. 

Two  sets  of  Saturn  argon  simulations  were  conducted  with  the  2-D  code 
PRISM.  The  first  set  of  simulations  were  conducted  by  changing  the  nozzle  de¬ 
sign  and  adding  progressively  more  tilt  angle.  We  used  the  PI  nozzle  design  that 
has  been  discussed  in  the  previous  section  (0.4  cm  nozzle  width)  except  that  the 
diode  length  was  taken  to  be  2  cm.  The  tilt  angles  covered  a  range  from  0°  to 
—12.5°.  This  series  of  simulations  thus  covered  designs  from  straight  nozzles  to 
nozzles  which  slightly  overcompensate  for  the  10°  expansion  of  the  puff  gas  as  it 
exits  the  nozzle  throat.  A  mass  loading  of  450  figm/cm  was  chosen,  since  as  dis¬ 
cussed  below  (see  Table  III),  this  is  near  the  mass  for  which  Saturn  achieves  the 
maximum  K-shell  yield.  The  results  are  not  expected  to  change  dramatically  with 
differing  mass  loadings. 

The  results  of  this  series  is  shown  in  Table  II.  It  can  be  seen  that  the  peak 
power  increases  by  a  factor  of  two  from  the  0°  to  the  —10°  case.  The  FWHM  (i.e., 
the  length  of  time  to  rise  from  one  half  peak  to  peak  and  then  fall  again  to  one 
half  peak)  also  falls  by  a  factor  of  two.  In  terms  of  these  two  measures,  the  tilted 


nozzle  design  is  obviously  superior  to  those  of  the  untilted  nozzle.  The  power  and 
pulse  width  values  in  Table  II  also  show  that  even  5°  of  tilt  substantially  improves 
the  implosion  quality  as  evidenced  by  the  higher  power  and  shorter  pulse  width. 
However,  the  K-shell  yield  is  roughly  30%  higl  for  the  tilted  case.  This  is  below  the 
yield  improvement  seen  in  the  simulations  of  the  PI  Double  Eagle  device  discussed 
above.  One  conjecture  for  this  discrepancy  between  the  two  sets  of  simulations  is 
the  different  nature  of  the  two  machines.  Due  to  the  low  machine  inductance  of 
Saturn,  the  current  peaks  well  before  the  time  of  implosion.  Therefore,  the  plasma 
is  not  so  much  driven  into  the  final  assembly  by  the  magnetic  field  pressure  but 
rather  implodes  under  its  own  inertia.  These  implosions  are  sometimes  referred  to 
as  “kinetic  energy”  implosions.  Double  Eagle,  by  contrast,  has  a  higher  inductance 
and  the  current  peaks  very  close  to  the  time  of  final  assembly.  In  this  case  the 
zippering  motion  can  be  driven  with  much  higher  gradients.  Consequently,  the 
emitting  region  can  be  more  localized.  The  result,  for  Double  Eagle,  is  that  there 
is  not  a  simple  tradeoff  between  peak  power  and  pulse  width.  With  Saturn,  the 
simulation  results  indicate  that  this  one-for-one  tradeoff  may  be  closer  to  what  is 
seen  in  the  0°  and  the  —10°  results  in  Table  II.  This  is  a  topic  which  needs  to  be 
investigated  further  in  the  future  in  order  to  understand  its  importance  in  DECADE 
and  fixture  machines. 


Tilt  Angle 
(degrees) 

Power 

(TW) 

Yield 

(kJ) 

FWHM 

(nsec) 

0 

1.7 

27 

13.5 

-5 

3.1 

31 

7 

-7.5 

3.2 

31 

6 

-10 

3.6 

33 

6 

-12.5 

3.5 

35 

5 

Table  II.  Saturn  Argon  K-shell  yields  as 
a  function  of  nozzle  tilt  angle  (450  /jgm/cm) 

Figure  4  shows  the  current,  input  generator  voltage,  and  radiative  power  for 
the  straight  nozzle  case.  The  same  quantities  are  shown  in  Fig.  5  for  the  case  with 
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—  10°  nozzle  tilt.  In  both  figures  it  can  be  seen  that  the  motional  impedance  of  the 
inwardly  moving  plasma  causes  the  current  to  roll  over  at  about  70  nsec.  The  result 
is  that  the  peak  current  is  between  8  MA  and  9  MA.  However,  the  most  dramatic 
difference  between  the  two  figures  are  6een  in  the  radiative  power  curves.  In  the 
tilted  case,  the  curve  is  symmetric  with  a  higher  peak  power  (roughly  3.6  TW).  The 
radiative  power  curve  in  the  untilted  case,  in  contrast,  is  not  as  symmetric,  broader 
at  its  FWHM,  and  has  slightly  less  than  one  half  of  the  peak  power  of  the  tilted 
nozzle. 

In  addition  to  the  tilted  nozzle  series  discussed  above,  a  series  of  simulations 
were  conducted  with  the  —10°  nozzle  design  and  various  mass  loadings.  In  the 
discussion  above,  it  was  seen  that  the  peak  power,  at  the  expense  of  a  shortened 
FWHM,  dramatically  increases  when  a  tilt  angle  is  used  to  counteract  the  expansion 
of  the  puff  gas  as  it  exits  the  nozzle.  Here,  we  detail  the  yield  with  mass  relationship 
for  the  tilted  nozzles.  The  results  obtained  using  a  1.25  cm  radius  nozzle  with  a  0.4 
cm  exit  width  and  a  10°  inward  tilt  are  given  in  the  table  below.  Again,  the  length 
of  the  puff  along  the  axis  was  taken  to  be  2  cm. 

When  these  simulations  were  made  with  different  values  of  enhanced  resistivity, 
we  found  that  K-shell  yield  remained  fairly  uniform  but  the  total  yield  increased 
(i.e.,  K-shell  plus  L-shell)  substantially.  Since  the  experimental  results  have  in¬ 
dicated  that  there  is  much  more  L-shell  radiation  than  is  typically  seen  in  the 
simulations,  we  chose  to  use  an  enhanced  resistivity  for  these  runs.  Evidently,  a 
large  fraction  of  this  L-shell  radiation  must  come  from  resistive  heating. 


Mass 

(10~®  gm/cm) 

K-shell  Yield 

(kJ) 

Total  Yield 

(U) 

300 

22 

122 

400 

32 

139 

450 

33 

148 

600 

31 

154 

Table  III.  Saturn  Argon  K-shell  yields  as 
a  function  of  mass  for  — 10°  tilt  nozzles 
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Additional  post-processing  was  done  using  the  complete  atomic  physics  package 
and  thin  radiation  model  (as  opposed  to  the  table  lookup  used  during  the  simula¬ 
tion)  at  selected  times  during  the  radiative  power  pulse  using  values  dumped  during 
the  simulation.  This  typically  requires  about  10  post-processing  runs,  based  on  one 
dump  per  nsec,  per  mass  load.  This  was  done  as  an  exercise  to  determine  the  accu¬ 
racy  of  using  the  table  lookup  scheme  which  is  must  less  expensive,  computationally 
speaking,  than  using  the  full  atomic  physics  package.  For  the  450  ftgm/cm  simu¬ 
lation,  this  method  gave  a  K-shell  yield  of  35  KJ  which  is  very  close  to  the  33  KJ 
obtained  from  the  lookup  values.  This  has  added  confidence  that  the  table  lookup 
method  is  validated. 


(c)  Summary 

In  summary,  when  the  zippering  effect  is  eliminated  in  argon  gas  puff  implosions 
by  using  a  tilted  nozzle  design,  the  K-shell  power  can  be  significantly  increased.  The 
K-shell  yield  is  also  increased  although  this  effect  is  more  dramatic  in  the  PI  Double 
Eagle  device  than  on  the  SNL  Saturn  device.  This  appears  to  be  due  to  the  fact  that 
the  plasma  which  participates  in  the  K-shell  radiation  is  more  efficiently  imploded 
in  the  tilted  nozzle  designs  and  reaches  higher  densities. 

The  use  of  tilted  nozzles  in  the  upcoming  DNA  devices  (i.e.,  DECADE  and 
follow  on  machines)  should  provide  for  more  optimized  implosions.  Further  exper¬ 
imentation  on  present  devices  (Double  Eagle  and  Saturn)  and  2-D  simulations  of 
nozzle  designs  will  provide  a  path  for  achieving  these  optimized  implosions. 

In  the  future  there  are  several  areas  which  should  be  explored  with  the  2- 
D  PRISM  code.  First,  the  discrepancy  in  the  enhancement  seen  in  the  Double 
Eagle  experiments  and  simulations  versus  the  Saturn  simulations  is  still  not  fully 
understood  although,  there  are  plausible  explanations.  It  is  important  that  this 
issue  be  resolved  because  the  two  devices  have  different  current  driver  characteristics 
and  this  is  an  important  design  issue  for  PRS  experiments  on  future  machines. 
Second  we  have  used  a  generic  nozzle  design  in  these  studies.  There  are  other 
parameters  of  the  nozzle  flow  which  may  be  important  including,  but  not  limited 
to,  the  profile  in  the  cross-section  of  the  gas.  The  ultimate  study  will  be  to  initialize 
the  code  with  results  input  directly  from  a  nozzle  design  code.  We  have  had  informal 
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preliminary  conversations  with  C.  Deeney  at  P.I.  about  doing  this.  Finally,  this  work 
represents  a  benchmarking  of  the  code  through  close  interaction  with  experiments 
at  P.I.  The  goal  in  the  upcoming  year  will  be  to  examine  optimum  PRS  designs  for 
DECADE  and  follow-on  pulse  power  machines. 
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Figure  Captions 


Figure  1.  Density  fill  plots  showing  the  initial  profile  (a)  and  the  profile  at  as¬ 
sembly  (b)  for  an  argon  gas  puff  implosion.  Untilted  nozzle  design. 

Figure  2.  Same  as  Fig.  1  but  for  an  implosion  initiated  with  a  -10°  tilted  nozzle. 

Figure  3.  Radiation  pulse  width  as  a  function  of  linear  mass  density  for  2-D  nozzle 
implosions  and  1-D  implosions  with  different  multipliers  on  Spitzer  resistivity. 

Figure  4.  Current,  input  generator  voltage,  and  radiative  K-shell  power  for  a 
Saturn  simulation.  This  argon  implosion  was  initiated  with  an  untilted  nozzle. 

Figure  5.  Same  as  Fig.  4  but  for  an  implosion  initiated  with  a  —10°  tilted  nozzle. 
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Nozzle:  2.5  cm  diam.,  0.4  cm  width,  0°  tilt 

Initial  Profile 


Figure  1.  (a) 
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Nozzle:  2.5  cm  diam.,  0.4  cm  width,  0°  tilt 

110  nsec 


Figure  1.  (b) 
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Nozzle:  2.5  cm  diam.,  0.4  cm  width,  -10°  tilt 

Initial  Profile 
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Radiation  Pulse  Width  as  a  Function  of  Mass 
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V.  Coupling  PRS  Loads  to  DECADE  Class, 

Long  Conduction  Time  POS  Inductive  Generators 

I.  Introduction 

The  introduction  of  inductive  energy  stores  (IES)  to  pulse  power  driven  Plasma  Radiation  Sources 
(PRS)  has  made  possible  an  order  of  magnitude  or  so  reduction  in  the  total  size  of  the  pulser  modules 
and  a  potentially  significant  decrease  in  the  risetime  of  the  output  current  pulses.  The  risetime 
decrease  would  make  available  the  added  option  of  smaller  radial  implosions  in  PRS  load  designs. 
NRL  is  involved  in  a  specific  effort  to  understand  power  flow  to  the  PRS  load  in  an  IES  machine  in 
support  of  the  DNA  DECADE  simulator.  As  part  of  that  program  we  are  developing  a  transmission 
line  model  that  resolves  the  Plasma  Opening  Switch  (POS)  module  into  a  series  of  shunts  in  the 
transmission  line  and  determines  the  electrical  properties  of  those  shunts  from  four  basic  theoretical 
models  (described  in  IL,  below).  This  work  is  intended  to  provide  first  a  simplified  but  comprehensive 
model  to  study  power  flow  to  any  PRS  operating  in  an  inductive  generator  environment,  and  later  some 
scaling  guidance  for  the  design  of  PRS  loads  in  DECADE.  In  particular,  useful  design  guidance  will 
require  an  understanding  of  several  issues,  such  as:  the  low  overall  system  efficiency,  the  possibility 
of  POS  restrike  and  other  MTTL  losses,  and  the  expected  peak  voltage  at  PRS  implosion. 

These  issues  are  in  fact  fundamental  to  IES/FRS  systems  and  essentially  unavoidable.  The 
operation  of  the  POS,  an  essential  IES  element,  tends  to  lower  the  general  efficiency  of  the  pulser. 
The  fraction  of  initial  stored  energy  delivered  to  the  load  can  be  low  (~  5%),  if  only  due  to  the  basic 
energy  transfer  constraints  associated  with  inductive  store  systems.  The  interaction  between  a  POS 
and  its  load  depends  upon  the  time  development  of  the  load  impedance.  The  efficiency  of  a  PRS 
implosion  depends  upon  the  quality  of  power  flow  to  the  immediate  vicinity  of  the  pinch.  Therefore, 
a  primary  question  is:  how  can  we  make  the  best  trade-off  between  the  high  voltages  developed  at 
PRS  implosion  and  the  amount  of  front  end  inductance  that  controls  the  probability  of  restrike  in  the 
POS?  If  restrike  occurs,  it  could  curtail  needed  power  flow  to  the  PRS.  If  one  can  speak  of  a  single 
"coupling  problem"  as  such,  it  is  this  PRS  voltage/front  end  impedance  trade-off.  The  experience  to 
date  would  imply  that,  as  the  present  switch  and  load  designs  are  scaled  to  DECADE,  neither  the 
PRS  nor  the  POS  can  be  optimally  designed  in  ignorance  of  one  another. 
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We  have  been  quite  successful  previously  in  using  simple  transmission  line  and  PRS  descriptions 
for  modeling  the  power  flow  to  PRS  loads  in  capacitively  driven  pulse  power  machines,  and  in  inves¬ 
tigating  the  consequences  of  various  load  and  machine  parameters  on  PRS  performance.  Early  EES 
investigations  involved  the  use  of  phenomenological  lumped  circuit  modules  to  mock  up  the  POS  — 
typically  as  a  time  dependent  resistance  with  a  conduction  delay.  The  work  presented  here  constitutes 
a  natural  extension  of  these  previous  studies,  but  different  in  that  it  proposes  several  new  theoretical 
descriptions  for  the  recent  observations  in  the  long  conduction  time  POS  experiments  done  on  Hawk 
at  NRL.  The  scope  of  this  summary  is  a  description  of  the  POS  model,  its  benchmark  with  the  Hawk 
experiments,  and  some  preliminary  observations  on  power  flow  and  POS/PRS  interactions. 

IL  Structure  and  Performance  of  the  POS  Model 

While  motivated  by  the  need  to  carry  out  PRS  load  design  calculations  for  DECADE  that  en¬ 
compass  all  the  relevant  power  flow  issues,  the  new  theory  initiatives  to  be  set  forth  are  rooted  both 
in  recent  observations  from  the  Hawk  POS  device  and  in  some  weaknesses  of  early  theoretical  mod¬ 
els  in  treating;  (i)  the  progression  of  current  during  the  POS  conduction  phase,  and  (ii)  the  cathode 
sheath.  In  some  early  work  the  (apparently  linear)  progression  in  time  of  the  current  channel  along 
the  cathode  could  be  described  as  a  process  of  matching  the  switch  current  to  the  injected  current, 
and  opening  was  viewed  as  a  consequence  of  running  out  of  new  area  for  conduction,  which  forced  ion 
current  to  be  drawn.  This  increased  ion  flow  would  erode  the  plasma  ions  to  open  a  channel.  In  this 
picture  the  local  cathode  area  and  plasma  density  set  an  upper  limit  to  the  current,  no  matter  what 
the  timescale  for  charging  the  storage  inductor.  However,  the  width  of  the  observed  current  channel 
seemed  to  argue  for  a  turbulent  diffusion  process. 

A  linear  current  progression  is  inconsistent  with  respect  to  both  observation  and  theory.  Mag¬ 
netic  probes  generally  show  a  very  broad  current  channel,  and  even  given  recent  observations  limiting 
the  spatial  resolution  of  the  probes,  any  apparent  turbulent  diffusion  indicated  by  the  wide  current 
front  should  p:ovide  a  current  progression  into  the  switch  with  a  speed  determined  by  the  physics  of 
the  anomalous  transport  processes.  If  such  a  diffusion  picture  is  accurate,  then  there  is  every  reason 
to  expect  a  nonlinear  progression  of  the  current  channel  through  the  switch.  No  basis  for  the  very 
special  case  of  a  linear  progression  has  been  established.  Moreover,  a  good  deal  of  evidence  to  date 
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suggests  that  the  actual  progress  of  the  current  channel  is  not  linear,  at  least  at  higher  conduction 
currents.  It  is  certainly  possible  therefore  that  the  switch  is  not  a  "current  limited*  device  at  all.  Cer¬ 
tainly  in  the  long  conduction  time,  high  current  limit,  the  nonlinear  scaling  of  conduction  time  and 
peak  current  is  potentially  describable  by  a  variety  of  snowplow  style  models.  The  Hawk  observations 
demonstrate  clearly  that  bulk  plasma  motion  and  number  density  depletion  are  crucial  processes  in 
the  long  conduction  time  POS. 

A  second  weakness  in  earlier  models  was  the  treatment  of  magnetic  insulation  in  the  cathode 
sheath  as  it  progressed  along  the  switch.  The  use  of  an  unmagnetized  bipolar  Child-Langmuir  cur¬ 
rent  density  over  the  whole  "active*  portion  of  the  cathode  is  inaccurate  and  ignores  the  important 
modifications  to  the  gap  space  charge  density  arising  in  a  magnetized  sheath.  As  shown  below,  only 
an  extraordinary  enhancement  of  electron  space  charge  can  draw  enough  ions  from  the  anode  surface 
of  the  sheath  to  produce  a  significant  level  of  erosion,  and  only  a  complete  description  of  the  ion 
emission  process  that  includes  a  presheath  model  can  establish  whether  or  not  erosion  can  occur. 
The  issue  must  be  addressed  by  a  model  sufficiently  complete  to  account  for  the  evolution  of  electron 
space  charge  over  the  whole  length  of  any  gap  or  depletion  channel  that  arises.  Unfortunately  the 
gap  and  emission  surface  structures  in  the  long  conduction  POS  are  not  well  known  -  the  "physical 
gap"  (inferred  from  density  depletion)  is  quite  a  bit  larger  than  the  "electrical  gap’  (inferred  from 
magnetic  insulation  conditions).  Some  mriel  accounting  for  the  buildup  of  electron  space  charge  in  a 
semi-colliaional  gap  is  needed,  either  to  justify  or  to  ignore  the  idea  of  erosion. 

New  POS  Theory  Elements 

In  addition  to  the  motivating  weaknesses  just  noted  there  are  as  well  some  motivating  obser¬ 
vations  that  argue  for  a  fresh  model  and  arise  from  the  Hawk  experiment.  The  principal  innovation 
in  Hawk  is  the  observation  of  the  axially  integrated  electron  line  density  over  the  entire  conduction 
and  opening  phase  at  various  radial  locations.  That  line  density  decays  with  time  over  a  gap  about 
1  cm  wide,  referred  to  here  as  a  "depletion  uhannel".  If  the  peak  current  prior  to  opening  is  denoted 
bond,  one  finds  that: 
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[  i.]  A  steady  drop  in  line  density  precedes  the  opening  event . 

[  ii.]  Ir.nn<i  scales  sublinearly  with  density,  a  n°'2'\ 

tiii]  The  observed  physical  width  D  of  the  depletion  channel  is  larger  than  the  mm  size  electrical 
width  inferred  from  the  opening  data. 

A  model  that  seeks  to  capture  these  observations  quantitatively  requires  several  characteristic 
structures  and  some  rules  governing  the  relationships  among  them.  The  basic  conception  of  such  a 
model,  for  a  coaxial  FOS,  is  shown  in  Fig.  1. 


Fig.  1  Conduction  Phase  POS  Structures 
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With  physical  structures  shown  to  the  right  of  the  center  line,  and  field  profiles  schematized  to  the 
left,  the  POS  is  pictured  here  in  the  conduction  phase.  Generator  current  circulates  from  anode 
to  cathode  through  a  Hall  current  channel  that  is  compressing  the  switch  plasma  and  propagating 
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toward  the  load.  In  the  depletion  channel  (shown  just  upstream)  one  expects  a  partially  magnetically 
insulated  Child-Langmuir  shunting  current  to  develop,  but  the  primary  conduction  current  is  in  the 
Hall  channel.  The  axial  extent  of  the  Hall  channel  is  determined  primarily  by  turbulent  diffusion 
processes  scattering  the  electron  flow  emitted  from  the  upstream  cathode  (interior  electrode).  The 
shunting  conduction  current  determines  the  line  voltage  by  attenuating  the  radial  electric  field  of  the 
coaxial  TEM  mode. 

The  POS  description  summarized  below  requires  four  theoretical  developments  to  build,  from 
the  concept  above,  a  mathematical  model  that  is  free  of  the  problems  already  mentioned,  and  is  also 
readily  extendable  to  deal  with  more  subtle  issues  later.  This  model  is  comprised  of  (i)  a  shunted 
transmission  line  picture  of  the  POS  electrodynamics,  (ii)  a  Hall  channel  snowplow  picture  of  current 
penetration  into  the  POS,  (iii)  a  warm,  magnetized  emission  theory  to  describe  the  magnetic  insulation 
of  the  depletion  channel,  and  (iv)  an  explosive  loss  of  force  equilibrium  to  initiate  the  opening. 

Transmission  Line  Representation 

The  normally  complex  2-D  electrodynamics  present  in  a  POS  can  be  approximated  by  the  evo¬ 
lution  of  a  TEM  mode  in  a  shunted  transmission  line.  The  shunted  transmission  line  description 
is  possible  because  one  can  simplify  the  electrodynamic  description  by  employing  axial  electric  field 
modes  which  allow  a  self-similar  response  of  the  POS  plasma  to  the  stress  imposed  by  the  generator. 
The  physical  character  of  the  transmission  line  in  the  POS  region  is  determined  by  the  effective  gap 
width  D  and  the  effective  cathode  radius  rr,  which  may  be  located  several  mm  into  the  physical  gap 
as  the  conduction  phase  develops.  The  absorbed  TEM  mode  propagates  down  this  line,  defined  by 
the  low  density  annular  "depletion  channel"  region  just  upstream  of  the  Hall  current  channel.  As 
the  field  energy  propagates  into  the  denser  plasma  of  the  Hall  channel,  it  develops  a  longitudinal 
component,  Ez,  viz. 

E  =  f  EjFM  +  zE‘z 

which  can  be  strong  enough  to  support  whatever  radial  Hall  shunt  current  (oc  Ez  1 )  is  required 
by  the  difference  in  line  currents  that  enter  and  leave  any  annular  slice  of  switch  plasma  in  the  Hall 
channel. 
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The  axial  self-similarity  of  the  plasma  response,  over  such  an  annular  slice  containing  the  radial 
Hall  current,  is  possible  because  the  radial  structure  of  the  axial  electric  field  conforms  to  a  particular 
set  of  TM  boundary  conditions.  It  is  not  obvious  that  such  a  mode  set  should  exist,  but  it  turns  out 
that  any  POS  bulk  plasma  density  profile  has  a  compatible  radial  profile  for  the  axial  electric  field 
which  admits  this  decompostion.  The  radial  Hall  current  density  Jr  can  then  be  described  by  a 
self-similar  form  which  varies  only  in  z,  viz.  ren^r,  2,  t.)  j3*  =  Jr(z,  t)  ■ 

Current  Penetration:  Conduction  =>  Kinematics 

Conversely,  when  (i)  the  axial  current  difference  is  taken  to  be  the  Hall  shunt  current  needed  to 
maintain  current  conservation  and  (ii)  an  axial  extent  for  the  Hall  channel  is  fixed  by  an  appropriate 
turbulent  diffusion  length,  then  the  local  axial  electric  field  can  be  uniquely  determined  over  the 
entire  length  of  the  switch  plasma.  But  now  the  electrodynamics  determines  the  kinematics  of  the 
Hall  channel  and  the  length  of  the  conduction  phase  is  determined  by  the  time  required  for  this 
snow-plow  like  Hall  cnannel  disturbance  to  propagate  the  length  of  the  switch  plasma  and  then  push 
through  the  rear  or  downstream  surface  to  open  tie  circuit.  In  this  process  the  ion  number  density 
is  compressed  axially  but  is  lost  radially,  and  therefore  the  conduction  phase  model  only  tracks  the 
axial  compression  until  an  initial  cell  size  is  squeezed  to  a  fixed  b.  rail  fraction  of  its  original  size. 

Is  the  POS  then  a  current  limited  device?  In  this  hydrodynamic  limit  it  would  not  show  any 
intrinsic  upper  limit  to  I r.„nd-  If  this  conduction  mechanism  iB  excited,  the  plasma  can  accept  an 
essentially  arbitrary  amount  of  current  into  such  a  "Hall  channel'  by  means  of  progressively  smaller 
charge  separations  at  progressively  higher  bulk  plasma  densities.  The  net  effective  stress  on  the  ion 
plasma  component  can  be  calculated  from  (i)  a  partial  force  balance  on  the  electron  plasma  component 
and  (ii)  a  fixed  amount  of  radial  shunt  current  generated  per  unit  axial  length.  The  current  that  can 
be  carried  by  this  process  is  as  large  as  ever  needed  for  the  theory  of  present  or  future  POS  devices. 
The  scale  size  for  the  current  solution  is  given  by  60.359  [r;?nOjj/101R][.2/2]  MA.  In  practice,  since 
the  typical  POS  ion  density  is  on  the  order  of  lO1'1*  and  the  charge  state  Z  is  of  order  one,  this 
normalizing  current  is  about  60  MA.  Since  very  small  fractions  of  this  current  are  encountered  in 
real  POS  devices,  the  resulting  fractional  charge  separations  in  the  Hall  channel  are  very  small. 


91 


Magnetized  Emission:  Gap  Equilibria  ==>  Erosion  Constraints 


The  required  magnetized  emission  picture  considers  the  dynamics  of  “free"  electrons,  viz.  those 
having  a  mean  free  path  on  the  order  of  or  greater  than  the  gap  size.  Unlike  the  vacuum  M1TL 
problem,  the  effective  emission  surface  in  the  POS  depletion  region  is  not  a  dense,  highly  conductive 
electrode  plasma,  but  rather  a  diffuse  plasma  cloud  with  a  sharp  radial  density  gradient.  Electron 
emission  can  be  space  charge  limited  in  the  depletion  region,  but  the  characteristic  length  for  the  decay 
of  any  net  charge  is  of  the  order  of  a  few  Debye  lengths,  or  more,  up  to  the  plasma  skin  depth  (if  the 
conducting  layer  is  microturbulent  or  exhibits  neutral  scattering).  Under  these  conditions  a  detached 
electron  current  would  not  be  expected,  and  the  emitted  population  of  electrons  can  only  exhibit  a 
few  relevant  orbit  classes  that  contribute  to  an  equilibrium  flow  in  the  ambient  magnetic  field  of  the 
depletion  region.  In  a  steady  Vlasov  limit,  there  must  exist  a  detailed  balance  between  electrons 
entering  a  volume  element  on  the  “cathode”  surface  along  a  given  trajectory,  those  leaving  the  surface 
along  that  trajectory,  and  those  electrons  resident  on  the  surface,  establishing  the  local  space  charge. 
The  radial  TEM  mode  electric  field  at  the  surface  then  establishes  a  local  charge  density  and  its  radial 
integral  defines  the  voltage  across  the  line.  The  fraction  of  the  electron  population  surviving  to  enter 
the  “anode"  surface,  rather  than  propagate  axially,  then  determines  the  local  shunt  current  allowed 
in  the  space  charge  limited  gap. 

Using  this  model,  one  may  calculate  the  impedance  of  a  POS  gap  in  the  depletion  region  accu¬ 
rately  over  the  whole  of  DECADE  relevant  parameter  regimes.  One  output  of  the  emission  model  is 
a  rule  for  the  space  charge  limited  shunt  impedance  to  be  used  at  any  point  in  a  depletion  channel  or 
sheath  when  it  is  partially  or  heavily  magnetically  insulated.  It  turns  out  that  the  gap  impedance, 
relative  to  the  usual  Child-Langmuir  impedance,  can  be  roughly  approximated  as  a  function  of  a 
single  variable,  the  product  of  cyclotron  frequency  and  sheath  free-fall  transit  time,  which  serves  as 
a  measure  of  the  magnetization. 

A  second  output  of  the  model  is  an  estimate  of  the  ion  current  drawn  from  the  cathode  side  of 
the  gap  as  a  function  of  the  local  magnetization.  As  shown  in  the  table  below,  as  the  magnetization 
increases,  the  normalized  electron  flux  drops,  but  the  relative  ion  current  enhancement  increases 
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dramatically.  The  final  (or  absolute)  ion  flux,  scaling  as  the  product  of  these  parameters,  is  not  much 
enhanced  for  any  degree  of  magnetization.  Since  the  largest  enhancement  factor  is  so  weak  the  ion 
flux  may  never  become  sufficiently  large  to  produce  erosion. 

Table  I.  Ion  Fluxes  and  Magnetic  Insulation 


Magnetization 

Electron  Flux 

Ion  Factor 

Final  Ion  Flux 

0.0 

0.81313 

0.96976 

0.78850 

0.25 

0.15198 

4.2370 

0.64394 

0.50 

0.15005 

4.5547 

0.68343 

0.75 

0.18714 

4.4787 

0.83814 

0.85 

0.22396 

4.3375 

0.97143 

1.0 

0.15551 

7.9983 

1.2438 

1.15 

0.6930 

1.25 

bhu i 

TTFTTTT 

0.4940 

Gap  Equilibria,  the  Role  of  Erosion,  and  Force  Balance 

In  Fig.  2  a  sheath  structure,  in  the  neighborhood  of  an  ’’anode’’  plasma  surface  bounding  the 
depletion  region,  has  been  sketched.  The  spatial  variable  £  denotes  a  length  scaled  by  the  gap  width 
T>.  The  ion  velocity  profile  V(£)  has  a  scale  speed  V„  determined  by  the  voltage  across  the  gap. 
Once  an  ion  flow  is  established  in  the  space  charge  limited  region,  the  density  is  determined  by  the 
space  charge  limited  (spatially  uniform)  ion  flux,  T„  and  the  ion  velocity  profile.  As  indicated,  the 
removal  of  ions  from  this  surface  is  controlled  by  the  relative  strength  of  two  competing  terms  in 
the  continuity  equation.  For  any  ion  velocity  profile  whatever,  fluid  kinematics  sets  a  speed  for  the 
apparent  movement  of  any  density  level  on  any  spatial  domain 


v  =  v.(£B)[i  + 


d(  In  n„ 


Clearly,  if  V  is  to  expand  because  the  gap  erodes,  then  the  density  level  must  recede,  viz. 


In  nn 
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Fig.  2  POS  Conditions  Needed  for  Gap  Erosion 
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and  the  density  level  moves  opposite  to  the  Sow!  This  is  a  totally  general  criterion  for  the  erosion  of 
an  emission  surface,  and  from  it  one  can  derive  two  important  results. 

First,  it  is  easy  to  see  that  in  regions  of  space  charge  limited  flow  an  erosion  process  is  impos¬ 
sible  because  for  any  of  the  equilibrium  flows  with  fixed  r„, 

^  In  Vn  =  _x 
d^\nn„  ~ 

Second,  it  is  clearly  the  presheath  region  that  requires  an  accurate  theoretical  treatment  to  assess 
the  speed  of  any  erosion  process.  In  the  presheath  the  ion  flow  makes  a  smooth  transition  from  the 
steady  interior  flow  u„  characterizing  the  source  plasma  to  the  free  fall  profile  determined  by  the  gap 
space  charge.  Only  if  the  velocity  gradient  in  the  presheath  solution  is  sufficiently  strong  will  the 
erosion  process  be  triggered;  the  magnitude  of  the  ion  flux  that  finally  leaves  to  enter  the  strictly  space 
charge  limited  region  is  only  a  boundary  condition  for  the  presheath  problem.  Establishing  erosion 
or  computing  an  expected  speed  for  the  erosion  process  must  proceed  from  a  complete  solution  over 
the  presheath  region  -  large  ion  flux  cannot  by  itself  establish  that  erosion  must  occur. 

For  a  experiment  like  Hawk,  a  first  approximation  to  the  presheath  problem  can  be  developed 
using  the  ion  fluxes  that  come  from  the  magnetized  emission  model.  The  scale  factor  called  "Final 
Ion  Flux"  in  the  table  above  enters  as  a  calibration  constant  for  the  emitted  ion  flow.  Balancing  the 
ion  flux  into  the  equilibrium  gap  with  a  thermal  flux  from  the  interior,  and  assuming  the  presheath 
spatial  scale  is  k  Debye  lengths,  the  time  scale  for  presheath  density  depletion  can  be  estimated  by 

Ion  Flux  Out  Thermal  Flux  In 


— dt  Inn* 


$ 

1.96  ~  9.79-10 'y/TJAvi 

743k 

Spatial  Scale 


and,  even  for  the  most  generous  choice  of  k=1  and  typical  Hawk  parameters  fT^slOeV.n.sl.O 
•1016,D=0.25  cm],  this  estimate  implies  that  Vaaj)  ~  5.0  MV  would  be  needed  to  overcome  the  ther¬ 
mal  flux.  For  Hawk  parameters  even  the  largest  ion  flux  enhancement  due  to  thermal,  magnetized 
emission  implies  weak  erosion. 
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If  erosion  is  not  a  strong  opening  mechanism,  then  what  additional  physics  will  determine  the 
size  of  the  depletion  region?  Is  the  size  likely  to  be  static,  and  what  conditions  might  force  it  to 
expand? 

Unlike  a  conventional  space  charge  limited  diode,  the  POS  depletion  region  gap  is  defined  by 
fluid  surfaces  that  are  free  to  move.  Thus  paired,  charged  layers  having  finite  dimension,  large  density 
gradients,  and  small  to  moderate  resistivity  are  considered  in  this  model  as  likely  souroes  of  electron 
and  ion  emission  into  the  thinner  material  between  them.  An  opposing  pair  of  such  surfaces  is  then 
subject  to  strong  attractive  forces  if  the  voltage  across  their  gap  is  large  and  strong  repulsive  forces 
if  the  intervening  magnetic  field  within  that  gap  is  large.  As  the  conduction  phase  progresses  the 
axial  compression  and  radial  exclusion  of  mass  by  the  passage  of  the  Hall  channel  can  leave  behind 
a  depletion  region  gap  at  radius  R  in  a  radial  force  balance 

Outward  Magnetic  Push  Inward  Electric  Pull 

M  In2  (i  +  |)  V  =  i  [(hin,Zlin,)2  -  ( JWZW )2]  , 

where  M  is  the  mass  of  an  annular  slice  of  "cathode"  or  "anode*  plasma,  while  and  Igap  are  the 
currents  along  the  plasma  surfaces  and  across  the  upstream  slice  of  the  Hall  channel,  respectively. 

Here  the  arrangement  of  forces  is  such  that,  at  fixed  line  current  and  voltage,  smaller  separations 
than  the  equilibrium  value  will  drive  the  gap  shut  and  larger  separations  will  drive  it  apart.  Cast  in 
terms  of  the  line  and  gap  currents  and  voltages,  the  equilibrium  constraint  becomes 

Geometrical  Electrical 

Rbz  (l  4-  7?)  .  ,  *£\  _  506  I ggj,  — .  . 

D>  V  &)~  y/V^Iiin*-  V' 

Typical  solutions  for  D,  for  the  currents  and  voltages  inferred  at  the  end  of  the  depleted  gap  region 
during  POS  simulations,  are  on  the  order  of  3  — ►  5  mm.  During  conduction  this  force  equilibrium 
can  be  established,  but  such  equilibria  are  unstable.  The  loss  of  the  shunt  current  Ig(ip  as  the  Hall 
snowplow  pushes  off  the  downstream  surface  of  the  POS  leaves  the  system  out  of  balance.  When 
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this  equilibrium  fails,  the  net  radial  forces  can  be  quite  strong  and  readlily  open  the  gap  on  a  10  ns 
timescale. 

The  onset  of  the  opening  phase,  the  gap  size,  and  the  opening  rate  of  the  POS  are  thus  determined 
by  the  constraints  imposed  by  the  transmission  line  current  and  voltage  profiles  on  an  electrostatic 
and  magnetostatic  equilibrium  between  the  surfaces  forming  the  effective  cathode  and  anode  in  the 
depletion  channel.  The  opening  is  caused  by  the  rupture  of  the  Hall  channel  as  it  breaks  free  and 
moves  downstream  of  the  POS  leaving  a  depletion  channel  with  an  equilibrium  width  already 
large  enough  to  become  magnetically  insulated  by  the  typical  currents  and  voltages.  As  an  effect 
enhancing  the  opening,  when  the  Hall  channel  breaks  free  of  the  bulk  POS  plasma,  the  failur  e  nf  the 
force  equilibrium  pushes  the  gap  open  further.  Erosion,  if  it  occurs  at  all  at  the  higher  densities  of 
Hawk-like  switches,  is  likely  to  be  only  a  second  order  effect  serving  to  strengthen  the  opening  rather 
than  cause  it. 

Conduction  Phase  Benchmarks  of  the  POS  Theory 

As  a  fundamental  benchmark  for  the  model  proposed  here,  a  sequence  of  POS  conduction  phase 
solutions  with  a  short  circuit  load  was  computed  with  varying  background  plasma  densities.  In  this 
case  the  length  of  the  switch  was  held  fixed,  as  was  the  radius  erf  the  depletion  channel.  The  peak 
conduction  current  in  the  model  POS  as  a  function  of  the  average  plasma  density  that  uniformly  fills 
the  gap  as  the  Marx  bank  begins  its  discharge  was  calculated.  Note  that  the  conduction  time,  and 
the  peak  conduction  current  are  fixed  predictions  in  the  model  once  the  density  has  been  chosen  from 
experiment;  the  results  are  shown  in  Fig.  3. 

In  the  Hawk  experiment  the  measured  densities  are  on  the  order  of  those  used  here  (»  2.5  1075 
ions/cm3)  and  a  Hawk  simulation  was  computed  with  the  model  for  three  cases.  In  the  first,  the 
charging  voltage  or  the  Marx  is  the  same  as  in  a  typical  Hawk  shot,  Rs640  kV.  One  density  calibration 
point  was  checked  and  the  predicted  conduction  time  was  ss  1/rsec.  In  the  second,  the  charging  voltage 
on  the  Marx  is  the  same  as  the  rated  peak  voltage  for  Hawk,  «720  kV.  In  the  third,  the  Marx  voltage 
has  been  doubled  to  admit  more  current  into  the  inductive  store  and  thus  gain  access  to  a  larger 
range  of  densities  while  keeping  the  conduction  time  within  the  quarter  cycle  time  of  the  bank. 
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Peak  Conduction  Current  [kA] 


Fig.  3.  Conduction  Current  Scaling 
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The  conduction  current  in  the  simulation  line  is  somewhat  larger  than  in  the  experiment,  perhaps  due 
to  the  systematic  uncertainty  in  the  density  correspondence  between  the  model  and  the  experiment 
The  simulation  and  the  Hawk  data  scale  with  density  more  slowly  than  the  r?0  25  rule,  and  a  faster 
charging  rate  for  the  inductive  store  shifts  the  model  scaling  curve  toward  a  larger  conduction  current 
scaling  like  The  kinematics  of  the  depletion  channel  may  not  be  perfectly  represented  by  the 

existing  model,  but  any  more  detailed  treatment  based  on  the  same  fundamentals  will  probably  show 
a  similar  trend.  To  the  extent  the  experiment  moves  away  from  this  slower  density  scaling  and  source 
voltage  dependence,  the  corrections  to  the  current  penetration  model  should  be  sought  in  terms  of 
more  detailed  microturbulence  resistivity  models  and  diffusion  theory. 
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m.  POS/PRS  Coupling  Behavior  -  Preliminary  Findings 

Among  the  many  possible  paths  to  an  understanding  of  the  coupling  of  stored  electrical  energy 
to  radiation  through  an  inductive  energy  store  and  a  PRS,  the  traditional  choice  has  focused  on 
kinetic  energy  delivery  to  the  imploding  annular  plasma  and  compression  to  a  given  fraction  of  the 
initial  radius.  Here  the  assessments  will  be  made  in  a  similar  way  but  with  some  minor  changes 
and  additional  considerations  arising  from  the  POS  model  being  employed.  In  particular,  while  the 
IES  (being  essentially  a  current  source)  is  a  very  natural  PRS  driver,  the  risetime  of  that  current 
into  the  load  will  depend  upon  the  intervening  inductance  and  the  associated  possibility  of  restrike 
in  the  POS  gap.  For  the  PRS  this  means  that  the  usual  similarity  of  energy  transfer  for  fixed  values 
of  initial  mass  and  radius  in  the  combination  mnR?n  may  be  modified  by  POS  restrike  effects  as  R„ 
changes  shift  the  front  end  inductance.  For  the  same  reason  the  choice  to  measure  kinetic  energy 
delivery  at  a  fixed  fraction  of  the  initial  radius  is  also  modified  here.  The  "gas  bag”  load  model  used 
here  undergoes  a  transition  from  an  annulus  to  a  solid  cylinder  at  a  prescribed  lower  radius  rather 
than  a  prescribed  fraction  of  the  initial  radius.  Subsequent  to  this  assembly  the  load  usually  implodes 
further  and  the  final  stagnation  radius  is  fair  measure  of  the  change  in  inductance  achieved.  The 
peak  load  kinetic  energy  is  recorded  shortly  after  the  time  of  assembly,  and  is  then  reduced  to  zero  at 
peak  compression.  By  selecting  two  fixed  radial  measures  rather  than  imploding  to  a  fixed  fraction 
of  the  initial  radius,  the  absolute  values  of  initial  and  final  radius  play  properly  independent  roles  in 
defining  the  downstream  inductance  seen  by  the  POS. 

The  PRS/POS  in  a  Hawk-like  Line 

With  peak  currents  on  the  order  of  0.75  MA  Hawk  is  at  best  a  weak  PRS  driver,  but  it  could  be 
equipped  with  a  gas  jet  front  end  and  used  as  a  research  facility.  As  a  first  step  in  examining  the 
coupling  problem  that  was  defined  in  the  Introduction,  consider  the  Hawk  pulser  fit  with  a  PRS  load 
characterized  by  an  outer  (return  current)  radius  of  3.0  cm,  a  pinch  length  of  2.0  cm,  and  a  initial 
load  radius  of  1.5  cm.  The  initial  Marx  charge  is  taken  to  be  720kV,  the  POS  conduction  time  is  set  to 
about  1150  ns  by  adjusting  the  switch  plasma  density.  Masses  smaller  than  ~  37.5  //gm/cm  can  be 
imploded  in  such  a  device,  with  the  sum  of  the  POS  conduction  and  the  implosion  time  less  than  the 
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quarter  cycle  time  of  about  1.25  //s.  Since  the  load  is  Ar  the  X-ray  yields  are  a  small  fraction  of  the 
pinch  energy,  even  scaling  the  emission  by  n2r,  and  this  provides  a  useful  test  case  precisely  because 
the  energy  is  fairly  well  bottled  up.  The  initial  mass  is  varied  from  15.0  — »  35.0  //gm/cm  and  various 
measures  of  energy  coupling  and  efficiency  are  developed.  It  is  found  that  the  overall  efficiency  of 
energy  into  the  load  rises  monotonically  with  load  mass  for  all  cases  the  pulser  is  able  to  implode 
within  the  quarter-cycle  time.  There  is  a  small  peak  in  the  kinetic  energy  coupling  to  the  load  as  the 
mass  is  varied. 

Table  II.  PRS  Coupling  in  a  Hawk  line 


Mass 

//gm/cm 

E  Kin 
[kJ] 

E  Jnd 
[kJ] 

teV] 

npk 

[1019cm~3] 

Y  K 
[J] 

Y  r 
[J] 

15.0 

1.82 

3.60 

1,280 

1.51 

1.32 

1.98 

17.5 

1.79 

4.04 

1,140 

1.68 

1.75 

2.86 

20.0 

1.82 

3.63 

1,040 

1.79 

2.14 

3.81 

22.5 

1.83 

3.97 

960 

2.04 

2.66 

5.24 

25.0 

1.81 

3.82 

899 

2.03 

3.01 

6.45 

27.5 

1.86 

3.71 

838 

2.26 

3.47 

8.24 

30.0 

1.84 

4.00 

798 

2.40 

3.97 

10  1 

32.5 

1.85 

3.70 

755 

2.39 

4.18 

11.8 

35.0 

1.86 

3.74 

715 

2.61 

4.60 

14.2 

The  electrical  measures  of  energy  delivery  aie  seen  to  be  roughly  constant  for  all  the  load  masses; 
but  the  load  velocities,  at  the  beginning  of  assembly  and  stagnation,  were  certainly  not  constant. 
What  is  seen  here  is  the  familiar  inductance  rule  in  action,  viz.  the  energy  to  the  implosion  varies 
directly  with  the  change  in  load  inductance  over  the  rundown.  The  variation  in  radiative  yield  then 
arises  because  each  distinct  mass  is  locked  onto  a  distinct  stagnation  path  by  the  thermodynamic 
constraints  of  the  implosion.  Because  the  radiation  yield  is  so  small,  these  implosion  paths  follow 
essentially  adiabatic  trajectories  in  the  pressure  and  density.  In  the  Hawk  electrical  environment  the 
implosion  dynamics  are  fixed  by  the  stagnation  adiabat  which  any  given  load  mass  is  placed  upon, 
while  the  energy  delivered  to  the  front  end  is  essentially  constant . 
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The  PRS/POS  in  a  DECADE  Line 


Consider  a  generic  DECADE  pulaer  fit  with  a  PRS  load  characterized  by  an  outer  (return  current) 
radius  of  3.125  cm,  a  pinch  length  of  2.0  cm,  and  a  initial  load  radius  of  3.0  cm.  The  initial  Marx 
charge  is  taken  to  be  750kV,  the  POS  conduction  time  is  set  to  about  390  ns  by  adjusting  the  switch 
plasma  density  to  6.25  -101S  ions/cm3.  Each  module  accumulates  about  2MA  in  the  POS  before 
the  opening  event  is  triggered  by  the  emergence  of  the  Hall  channel  from  the  downstream  plasma 
boundary.  The  risetime  of  current  into  the  load  is  now  slower,  RJ50-55  ns,  and  the  initial  voltage 
on  the  load  generally  spikes  to  about  1700  kV  during  the  opening  phase.  The  peak  POS  voltage  is 
essentially  the  same  as  the  peak  load  voltage,  and  a  second  (weaker)  spike  in  both  occurs  as  the 
implosion  stagnates.  The  second  spike  is  generally  not  sufficient  to  spoil  the  magnetic  insulation  of 
the  POS  gap,  but  much  more  detail  about  the  front  end  must  be  incorporated  to  explore  the  question 
definitively. 

Masses  smaller  than  «  4.0  mgm/cm  can  be  imploded  in  such  a  device,  with  the  sum  of  the 
transfer  capacitor  charging  time,  the  POS  conduction  and  the  implosion  time  less  than  the  quarter 
cycle  time.  Recent  1-D  and  2-D  implosion  calculations  at  DECADE  current  levels  have  indicated  the 
persistence  of  a  rather  high  level  of  preheating  and  a  consequent  hollowness  to  the  implosions.  A 
gas  bag  model  can  only  simulate  this  through  the  addition  of  viscous  dissipation  if  the  resistivity 
is  to  remain  classical,  and  the  timing  of  the  enhanced  heating  is  to  correspond  with  the  earliest 
decelerations  in  the  stagnation  process.  As  a  consequence  the  model  load  is  equipped  here  with  an 
enhanced  ion  relaxation  time  which  provides  a  sufficiently  strong  viscous  preheat  to  the  implosion 
calculation  to  obtain  an  early  average  temperature  similar  to  that  in  the  1-D  and  2-D  MHD  models. 
The  relaxation  time  is  chosen  to  obtain  a  viscous  heating  rate  competitive  with,  but  not  exceeding, 
the  compressions!  heating  rate  that  sets  the  models  overall  stagnation  timescale.  The  load  material 
is  Ar  and  the  X-ray  yields  are  now  a  more  sizable  fraction  of  the  pinch  energy. 

If  the  initial  mass  is  varied  from  4.0  — +  0.40  mgm/cm  then  the  various  measures  of  energy 
coupling  and  efficiency  discussed  above  can  be  summarized  in  the  table  shown. 
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Table  HI.  PRS  Coupling  in  a  DECADE  Line 


Mass 

mgm/cra 

Eft'in 

[kJ] 

E/w 

[kJ] 

AL 

[nH] 

T  r* 
teV] 

KSH 

KM 

Y  r 

IkJl 

4.0 

440.0 

816.0 

12.8 

146.0 

l.n 

0.00 

178.0 

3.0 

379.0 

1,080 

19.8 

845.0 

28.7 

26.5 

1,600 

2.0 

316.0 

926.0 

20.4 

925.0 

24.9 

38.2 

1,320 

1.0 

237.0 

680.0 

20.3 

735.0 

11.9 

50.4 

877.0 

0.8 

226.0 

589.0 

18.6 

1,190 

4.20 

114.0 

741.0 

0.7 

213.0 

515.0 

17.2 

1,510 

1.78 

138.0 

368.0 

0.6 

193.0 

446.0 

16.6 

1,870 

1.12 

131.0 

204.0 

0.5 

178.0 

404.0 

16.3 

2,450 

0.804 

93.0 

119.0 

0.4 

162.0 

365.0 

16.1 

3,220 

0.597 

51.7 

61.1 

The  electrical  measures  of  the  energy  delivered  to  the  front  end  are  seen  to  peak  with  the  total  yield, 
and,  while  the  driving  current  at  the  time  of  stagnation  was  essentially  constant,  the  peak  kinetic 
energy  delivered  drops  off  less  rapidly  than  the  mass,  indicating  progressively  higher  implosion  speeds 
as  the  mass  falls.  The  K-shell  yield  peaks  at  a  mass  slightly  lower  than  that  which  optimizes  the  total 
yield.  The  peak  in  K-shell  yield  with  lower  mass  occurs  because  the  very  light  loads  get  hotter  and 
turn  around  faster  at  lower  density.  The  lower  yields  at  large  mass  are  simply  a  matter  of  insufficient 
energy  transfer. 

IV.  POS/PRS  Coupling  Behavior  -  Conclusions 

The  results  detailed  above  confirm,  once  again,  that  a  DECADE  class  machine  can  couple  energy 
to  a  PRS  load  and  do  so  in  a  relatively  efficient  manner.  From  these  calculations  one  would  expect 
between  0.5  — - ►  1.5  MJ  to  be  coupled  forward  of  the  POS  modules  into  the  front  end.  More  important 
than  the  power  flow  details  in  this  preliminary  study  however  is  the  result  that  a  POS  model  anchored 
in  a  few  simple  hypotheses  can  be  successful  at  DECADE  power  levels  and  properly  responsive  to 
the  electrical  environment  generated  by  the  PRS  and  any  front  end  Mil'L  loss  processes.  Several 
important  results  should  be  noted  from  the  POS  model: 

i.  Self-similar  Hall  current  channels  can  exist  for  any  bulk  plasma  density  profile. 
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ii.  These  Hall  current  channel  can  carry  very  large  currents,  implying  the  POS  is  not  a  current 
limited  device. 

iii.  The  model  captures  the  conduction  time,  peak  conduction  current,  opening  time,  and  the  ’ate 
time  impedance. 

iv.  The  model  exhibits  I rnnd  scaling  with  density  and  Marx  voltage  that  is  compatible  with  experi¬ 
mental  data. 

v.  The  concepts  that  work  to  explain  the  conduction  phase  have  important  implications  for  the 
opening  phase  and  the  viability  of  any  erosion  hypothesis.  In  particular,  any  erosion  speed 
calculation  must  account  for  the  presheath  structures;  and,  the  emission  theory  appropriate  to 
the  dense  Hawk  POS  conduction  phase  predicts  weak  ion  flux  enhancement. 

vi.  A  causal  opening  mechanism  has  been  identified  that  is  distinct  from  erosion.  In  particular,  gap 
equilibria  can  exist  which  predict  the  inferred  electrical  gap  size,  and  the  loss  of  force  balance  in 
these  equilibria  at  the  end  of  the  conduction  phase  predicts  gap  accelerations  and  opening  times 
compatible  with  experiment. 

The  POS  model  can  be  viewed  as  the  union  of  three  distinct,  but  intimately  related  parts.  Each 
element:  magnetized  emission  theory.  Hall  channel  penetration,  and  gap  dynamics,  can  and  should 
be  refined  much  further,  but  the  need  to  do  so  will  be  based  upon  the  performance  relative  to  existing 
and  further  POS  experiments.  However  these  elements  evolve,  the  results  can  now  be  installed  in  an 
integrated  numerical  framework  that  works  rather  well  as  a  model  of  the  entire  DECADE  machine. 

The  coupling  problem  posed  at  the  beginning  of  this  work  has  also  been  addressed  in  the  context 
of  a  particular  generic  machine  design,  with  the  result  that  POS  restrike  appears  unlikely  to  be  a 
major  concern.  While  some  POS  restrike  does  occur  it  generally  is  too  little  and  too  late  to  have 
much  effect  on  the  energy  delivered  to  the  load.  Naturally  this  is  a  very  preliminary  assessment 
and  it  must  be  refined  further  with  detailed  machine  design  parameters  and  careful!  attention  to 
the  experimental  results  as  they  are  available.  In  addition  to  the  painstaking  collection  of  machine 
details  and  further  POS  comparisons,  future  efforts  will  also  include  a  second  modeling  alternative 
for  the  DECADE  line,  viz.  the  combination  of  modules  to  a  single  parallel  line  with  an  appropriate 
POS  segment. 
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VI.  L-Shell  Model  Development  and  Scaling 


The  development  of  x-ray  diagnostics  in  the  L-shell  requires  an  extensive 
atomic  data  base.  These  data  include  energy  levels,  collision  strengths, 
photoionization  cross  sections,  and  dielectronic  recombination  branching 
ratios.  When  accurate  data  is  obtained,  it  can  allow  us  to  develop  scaling 
relations  to  obtain  more  reliable  and  more  extensive  data  files.  This  data 
processing  is  also  necessary  for  scaling  PRS  load  dynamics.  We  are  continuing 
the  process  of  slowly  upgrading  and  evaluating  our  L-shell  atomic  data  base, 
since  these  data  in  our  current  data  base  were  obtained  only  in  an  approximate 
way  and  are  not  very  accurate  or  complete.  We  have  developed  much  improved 
calculations  of  dielectronic  recombination  (DR)  rate  coefficients  over  the 
last  few  years.  A  dataset  was  completed  last  year  of  DR  ratf>'  for 
recombination  from  oxygen-like  ions  for  Ar,  Ti,  Pe,  and  Se.  These  DR 
calculations  involve  detailed  calculations  of  energy  levels  and  autoionization 
as  well  as  radiative  rates.  We  are  in  the  process  of  demonstrating  their  Z- 
scalability.  We  are  trying  to  develop  scaling  relations  for  the  singly  as 
well  as  doubly  excited  energy  levels,  autoionization  rates,  radiative  rates 
and  state  specific  DR  branching  ratios  following  our  scaling  of  dielectronic 
recombination  data  in  fluorine-like  ions  (Phys.  Rev.  42,  2640(1990)). 

By  carrying  out  the  process  of  scaling  using  three  and  four-parameter  fits 
to  the  oxygen-like  DR  data,  one  not  only  efficiently  extends  the  data  base  to 
all  elements  in  the  fourth  row  of  the  periodic  table,  one  also  obtains  a  check 
on  the  accuracy  of  the  data  of  the  Ar,  Ti,  Fe,  and  Se  databases  that  were 
calculated  in  detail.  While  most  of  the  data  can  be  fitted  with  extreme 
accuracy,  there  are  a  few  anomalies.  These  are  shown  in  Figs.  1,  2,  and  3 


using  a  few  of  the  states  that  make  important  contributions  to  the  DR  process. 

3  2 

Fig.  1  shows  how  the  energies  scale  for  the  doubly  excited  states:  2p  3s  , 
4  2  4 

2p  3d  ,  and  2p  3snl(4<n<8).  The  accuracy  of  the  latter  4  data  points  is  seen 
to  be  less  than  that  of  the  first  8  data  points.  The  least  squares  fit 
presumably  improves  upon  the  accuracy  of  the  latter  4  points.  In  Fig.  2  we 
see  that  the  total  autoionization  rates  for  the  states  2pJ3s*  and  2p*3dz  scale 
very  accurately.  The  data  for  the  autoionization  rates  of  the  2s2p^nl(5<n<10) 
state  do  not  scale  as  accurately.  Similarly,  Fig.  3  shows  that  the  total 
radiative  rates  for  the  states  Zp^s^  and  2p *3dz  scale  very  accurately,  but 
that  the  radiative  rates  for  the  2s2p^nl(5<n<10)  state  show  small  errors.  Our 
next  effort  was  thus  to  track  down  sources  of  the  few  erroneous  data, 
especially  if  human  errors  in  bookkeeping  such  a  large  number  of  input  and 
output  data  had  been  made.  Ve  are  going  through  a  thorough  checking  of  our 
database  for  oxygen-like  DR  recombination  data.  A  few  mistakes  have  been 
found  and  corrected  and  our  error  investigation  is  nearing  completion. 

We  have  also  recently  completed  DR  rate  calculations  for  oxygen-like 
molybdenum  (Mo)  in  order  to  determine  the  accuracy  of  extrapolating  the 
scaling  relations.  Once  we  have  finished  the  Z-scaling  relations  using  the  0- 
to  F-like  Ar,  Ti,  Fe  and  Se  data,  we  will  be  able  to  test  them  by  comparing 
the  Z-scaled  data  with  the  explicitly  calculated  Mo  data.  We  have  also  just 
begun  to  calculate  N-  to  0-like  DR  data  to  continue  to  fill  out  the  L-shell 
dataset. 

During  the  development  of  a  detailed  atomic  model  for  neon-like  selenium 
that  utilizes  the  above  DR  data,  we  discovered  that  our  in-house  capabilities 
for  calculating  the  necessary  collision  strengths  and  photoionization  cross 
sections  for  the  model  were  greatly  inadequate.  The  Branch  collision  strength 
capabilities  had  been  developed  for  K-shell  modeling  and  were  unreliable  for 
the  L-shell.  The  neon-like  selenium  model  was  completed  only  by  tediously 
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utilizing  data  that  was  available  in  published  x-ray  laser  literature. 
However,  the  reliability  of  our  DR  rate  calculations  must  be  matched  by  a 
similar  high  accuracy  and  consistency  for  our  collision  strength  and 
photoionization  data.  We  thus  investigated  the  possibilities  of  obtaining  or 
accessing  collision  codes  such  as  the  one  put  together  by  Zhang  and  Sampson  at 
Penn  State  University  or  even  a  modified  version  of  Cowan's  code  which  can  be 
used  to  accurately  calculate  the  needed  large  number  of  L-shell  collision 
strengths. 

We  have  now  acquired  atomic  codes,  developed  by  the  researchers  at  Los 
Alamos  Laboratory,  to  calculate  accurate  collision  strengths  and 
photoionization  cross  sections  (and  in  turn  radiative  recombination  rates). 
Using  these  codes,  we  should  be  able  to  generate  the  large  datasets  that  are 
routinely  needed  for  L-shell  ionization  balance  calculations.  These  codes 
include  the  structure  codes  of  R.  D.  Covan  and  are  available  on  the  Los  Alamos 
Cray.  We  have  obtained  printed  documents  for  properly  using  the  codes.  They 
include  1)  CATS,  a  modified  version  of  Cowan's  Atomic  Structure  Code,  2)  ACE, 
Another  Collisional  Excitation  code,  3)  LINE,  a  code  for  computing  atomic  LINE 
spectra,  and  4)  TAPS,  a  Theoretical  Atomic  Physics  code  which  is  used  to 
handle  and  and  display  the  data  generated  by  the  three  previous  codes.  The 
purpose  of  the  atomic  physics  code  development  effort  at  Los  Alamos  was  to 
bridge  the  gap  between  simple  models  and  sparse  sets  of  more  accurate  data  by 
providing  a  collection  of  computer  programs  for  conveniently  calculating  the 
quality  atomic  data  needed  for  ionization  balance  calculations  and 
diagnostics.  These  codes  are  easy  to  run,  and  their  required  input  has  been 
minimized.  The  fact  that  they  generate  atomic  physics  data  in  the  greatest 
detail  possible  as  well  as  perform  averages  when  such  detail  is  not  necessary 
makes  them  very  useful.  For  example,  the  programs  store  data  based  on  fine 
structure  levels,  but  if  quantities  based  on  multiplet  terms  or  configurations 
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are  desired,  then  the  fine  structure  data  is  averaged.  We  are  beginning  to 
use  these  codes  by  learning  hov  to  set  up  the  atomic  structure  files, 
input/output  commands,  code  action  commands,  data  file  management  commands  and 
the  commands  needed  to  calculate  multiplet  strengths  and  excitation  cross 
sections.  We  are  also  exploring  several  physics  issues  such  as  vhat  kinds  of 
wave  functions  or  what  approximations  need  to  be  used  in  these  codes  for  our 
particular  applications.  We  can  now  use  the  Los  Alamos  codes  to  calculate 
collision  strengths. 

In  order  to  test  the  Los  Alamos  collisional  codes  we  calculated  a  few 
neon-like  selenium  and  iron  collision  strengths.  The  results  obtained  agreed 
with  other  published  and  presumably  more  accurate  collision  strengths  except 
at  very  low  energies.  We  therefore  investigated  the  use  of  ACE  (as  mentioned 
before),  which  uses  a  distorted  wave  approximation  and  gives  more  accurate 
collision  strengths  than  CATS.  However  this  code  takes  a  large  amount  of 
computer  time.  The  researches  at  Los  Alamos  Laboratory  are  still  in  the 
process  of  optimizing  it  by  adding  options  for  calculating  more  atomic 
parameters.  We  are  presently  communicating  with  them  and  as  more  information 
becomes  available  we  will  incorporate  it  accordingly  in  our  calculations. 
Since  we  need  to  generate  large  numbers  of  collision  strengths,  the  ACE  code 
may  not  be  very  economical  and  practical.  Thus,  at  present,  we  are  using  the 
PWB  approximation  for  most  of  the  collision  strength  calculations  and  using 
ACE  for  the  most  important  diagnostic  transitions,  especially  at  low  energies. 

We  have  used  the  plane  wave  Born  (PWB)  method  of  Cowan  to  calculate  the 
collision  strengths  for  N-,  0-  and  F-like  selenium.  Fig. 4  shows  the  energy 
level  diagram  for  the  F-like  system  that  we  are  currently  using  for  our 
database  development.  The  other  L-shell  ionization  stages  have  a  similar 
structure:  only  the  degeneracy  of  the  states  changes  from  one  ionization  stage 
to  another.  The  collisional  excitations  A,  B,  and  C  (labeled  in  Fig. 4)  from 
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the  ground  state  to  the  3s,  3p  and  3d  excited  states  respectively  are 
especially  important.  In  Figs.  5  and  6  we  show  the  calculated  collision 
strengths  (solid  lines)  obtained  from  CATS  for  these  transitions  for  F-like 
and  O-like  selenium  and  compare  them  with  those  that  had  been  previously 
calculated  (dotted  lines)  from  our  adapted  in-house  codes.  Fig.  5  shows  that 
the  collision  strengths  for  the  dipole  transitions  (excitation  from  the  ground 
state  to  3s  and  3d  states)  in  the  CATS  calculations  are  close  to  those 
calculated  previously  except  at  very  low  energies.  Agreement  ij  within  a 
factor  of  two  for  the  monopole  (ground  state  to  the  3p  excited  state) 
transition.  However,  there  is  no  agreement  between  the  present  and  previous 
calculations  of  collision  strengths  for  the  same  transitions  in  O-like 
selenium  as  can  be  seen  from  Fig.  6.  Since  the  CATS  collision  strengths  are 
obtained  in  a  more  consistent  fashion  and  the  collision  strengths  for  the  same 
transition  (for  F-like  and  O-like  systems)  show  similar  behavior  we  believe 
that  our  present  results  are  much  more  accurate  and  should  be  scalable.  The 
previous  collision  strength  data  for  O-like  selenium,  on  the  other  hand,  were 
especially  inaccurate  and  unscalable. 

It  is  known  that  collision  strengths  decrease  as  the  nuclear  charge  Z 
increases.  We  are  in  the  process  of  obtaining  collision  strengths  from  CATS 
for  the  same  ionization  stages  and  for  the  same  transitions  for  iron  as  we 
have  done  for  selenium.  Once  these  data  are  obtained,  we  will  be  able  to 
verify  the  Z  dependence  of  collision  strengths  as  described  in  our  1990  final 
report  and  eventually  scale  the  data  to  obtain  collision  strengths  initially 
for  any  Z  between  iron  and  selenium  and  possibly  beyond. 

Presently  our  L-shell  data  base  does  not  include  any  radiative 
recombination  (RR)  rates  to  excited  states.  In  order  to  have  an  upgraded  and 
complete  model  we  have  to  eventually  include  such  rates.  There  is  also 
available  an  ionization  code  caxled  GIPPER  on  the  Los  Alamos  Cray.  However, 
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this  code  is  still  being  developed  for  the  calculation  of  more  improved 
collisional  ionization  rates,  and  there  is  no  available  written  document  on 
it.  We  can,  however,  use  it  to  calculate  photoionization  cross  sections; 
thus,  we  intend  to  learn  to  use  GIPPER  in  the  future  to  calculate 
photoionization  cross  sections  and  then,  using  detailed  balance,  calculate  the 
radiative  recombination  rate  coefficients  needed  to  compute  the  continuum 
emissions  of  importance  in  the  L-shell. 
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Fig.  6.  Collision  strengths  of  the  3s,  3p  and  3d  excited  states 
of  O-like  selenium  for  excitations  from  the  ground  state  of  the 
ion.  The  solid  lines  represent  the  present  calculations  and  the 
dotted  lines  represent  previous  calculations. 


